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AMBIGUITY -RESISTANT  THREE-  AND  FOUR-CHANNEL  INTERFEROMETERS 


1.0  INTRODUCTION 

Electronics  support  measures  (ESM)  systems  used  in  military  applications  such  as 
reconnaissance  and  threat  reaction  often  require  estimates  of  the  angle  of  arrival  (AOA) 
of  radio-frequency  signals  from  distant  emitters.  Accurate  real-time  AOAs— to  1°  rms 
resolution  and  accuracy  or  better— facilitate  efficient  emitter  sorting  routines  aa  well  as 
establish  emitter  locations  through  the  processing  of  successive  relative  bearing  measure- 
ments. For  ESM  use,  increasing  emphasis  is  being  placed  on  the  implementation  of  AOA 
techniques  that  can  provide  azimuth  coverage  over  wide  fields  of  view,  approaching  ±60° 
with  respect  to  the  direction-finding  (DF)  system  boresight.  Multiple-element  phase-only 
interferometers  processing  electrical  phase  differences  between  the  signals  received  at 
spaced  apertures  turn  out  to  be  excellent  for  achieving  accurate  AOAs  for  certain  ESM 
requirements.  These  techniques  appear  particularly  applicable  to  implementing  high- 
accuracy  DF  from  airborne  platforms  in  the  frequency  bands  of  interest  above  1 GHz. 

The  three  major  goals  of  this  report  are 

• To  clarify  some  past  misconceptions  concerning  the  theory  of  multielement 
interferometers. 

• To  expound  a general  theory  of  interferometers  which  may  stimulate  re- 
search into  the  less- well-qualified  aspects  of  phase-only  Interferometers,  such  as  the  con- 
straints imposed  by  operation  in  severe  multipath. 

• To  provide  the  ESM  system  designer  with  exact,  readily  applied  techniques 
for  obtaining  the  lowest  probability  of  ambiguity  for  a given  overall  array  length  in  three- 
and  four-element  arrays. 

This  report  treats  the  phase-only  interferometer  exclusively.  Interferometers  that 
process  relative  amplitude  information  as  well  as  phase  difference  information  from  multi- 
ple apertures  - techniques  prevalent  in  radio  astronomy— are  beyond  the  scope  of  this 
report.  Attention  will  further  be  constrained  to  line  arrays.  Within  those  apparently 
severe  restrictions  of  scope,  there  are  many  areas  of  applications  for  interferometers. 

The  reason  for  fixing  attention  exclusively  on  the  phase-only  Interferometer  is  that 
in  ESM  systems,  the  designer  often  must  maximize  the  instantaneous  (nonscanned)  field 
of  view.  He  is  faced  with  the  choice  of  implementing  directivo-gain  antenna/receiver 
channels,  relatively  nondirective  antennas  driving  phase-only  channels,  or  combinations 
of  these  approaches.  For  those  requirements  in  which  the  loss  in  system  detection  sensi- 
tivity from  nondirective  apertures  is  acceptable,  interferometer  techniques  that  inherently 
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achieve  a large  variation  in  output  parameter  for  a small  change  in  input  hearing  angle 
(large  gradient)  with  comparatively  few  channels  are  quite  suitable. 

The  simplest  conceivable  interferometer  uses  two  essentially  omnidirectional  antennas 
and  must  employ  supplementary  techniques  to  discriminate  against  energy  arriving  from 
the  rear  hemisphere  vs  energy  arriving  from  the  front.  In  this  elemental  interferometer, 
the  maximum  allowable  aperture  spacing  without  ambiguous  indications  of  AOA  is  one- 
half  wavelength  at  the  operating  frequency.  Arrays  of  n elements  (a  > 3)  '.'ovide  (n  - 1) 
phase  differences  that  can  be  used  to  reBolve  these  ambiguities,  t:  *hough  the  spadngs 

between  one  or  more  pairs  of  adjacent  channels  exceeds  one-half  w..voiength.  Multielement 
arrays  must  be  used  in  order  to  achieve  both  good  angular  resolution  and  low  probability 
of  ambiguity  for  realistic  values  of  receiving  channel  phase  errors. 


The  amount  of  published  material  on  radio  direction  finding  is  incredible.  Travers 
and  Hixon  [1]  have  assembled  abstracts  on  the  DF  literature  (including  interferometers) 
covering  the  period  1899-1965.  More  recently,  Barton  [2-4]  has  included  entries  on 
interferometers  in  his  index  (and  supplements)  on  material  published  in  the  I.E.E.E.  Tram • 
actions  on  Aerospace  and  Electronic  Systems  and  its  predecessor  publications.  The  em- 
phasis of  recent  Russian  work  available  in  translation  [6-7]  is  oriented  more  toward  “fre- 
quency" interferometers  (multifrequency  ranging  schemes)  than  toward  addressing  the 
“spatial”  interferometer  synthesis  problem.  Apparently,  no  work  clarifying  past  incon- 
sistencies in  three-element  arrays,  and  oxtending  the  theory  to  optimum  arrays  of  four 
elements,  has  appeared  to  date. 


Limiting  the  scope  of  this  report  to  a consideration  of  only  three-  and  four-element 
systems  will  not  unduly  restrict  design  freedom.  These  arrays  are  suitable  for  many  re- 
quirements, us  the  following  example  will  Bhow. 

The  example  pertains  to  locating  surface  emitters  from  an  airborne  platform  and  is 
offered  to  illustrate  some  of  the  tradeoffs  between  using  arrays  with  either  three  or  four 
elements.  Figure  1-1  is  a plan  view  of  the  geometry  typical  of  an  airborne  collection  sys- 
tem taking  a series  of  relative  bearings  on  a distant  emitter.  There  are  many  treatments 
of  the  factors  affecting  location  system  performance  [8-11]  for  the  problem  depicted  in 
the  illustration.  The  example  shown  is  Butterly’s  [8]  “asymmetrical  15°  sector.”  The 
emitter  range,  normal  to  the  assumed  straight-line  flight  track,  is  100  n.mi.  (186.2  km). 

A flat-earth  approximation  is  made,  and  errors  because  of  altitude  are  neglected.  The 
collection  platform  takes  16  successive  bearings  (forced  to  be  equally  spaced  in  Butterly’s 
analysis,  for  convenience,  with  no  great  loss  in  generality).  The  first  DF  cut  is  made  at 
an  angle  of  46°  right  of  the  array  boresight;  the  last  is  takei  at  a bearing  of  30°.  But- 
terly's  analysis  shows  that  to  restrict  the  area  of  uncertain*  of  emitter  location  to  n 
n.mi.2  (3.43kt  km2)  to  a probability  of  0.96  requires  accuracies  of  0.21°  rms  on  the  bear- 
ings over  this  15°  runge  in  angle. 


The  spatial  accuracy  in  degrees  rms,  oyu,  of  an  interferometer  is  related  to  the  elec- 
trical phase  error  in  the  largest-spaced  pair  of  antennas  by  the  well-known  expression 


- °<t> 
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Flu.  1-1—  Problem  geometry  for  emitter  location  example 


= electrical  phase  error,  degrees  rms, 

(}  = 27 r/X  = phase  constant  in  the  medium  of  propagation, 
d - channel-pair  spacing,  longest-spaced  channels. 

If  oa,  = 12°  rms  (zero-mean)  and  if  9a  is  in  the  vicinity  of  37.B°  (mean  of  48* _ and 
30°)  to  achieve  an  accuracy  in  AOA  of  0.21°  rms  (zero-mean)  requires  an  overall  base- 
line length  d = 11.5  wavelengths  = 23  half-wavelengths.  (Note:  In  developing  the  theory 
of  low  probability-of-ambiguity  arrays,  it  is  convenient  to  work  with  spacing * 2 

in  integer  hulf-wavelengths.)  The  optimum  f/tree-element  array  employs  a channel  1 to  l 
spacing ^jf  12  hdf-wavelcngths,  and  a channel  2 to  3 spacing  of  H half-wavelengths  For 
aPreadily  obtained  = 12°  rms,  the  probability  of  ambiguity  of  this  array  is  45%.  T is 
is  clearly  unacceptable  performance. 

The  optimum  four- clement  array,  synthesized  according  to  the  principles  presented 
in  Secs  4.0  and  5.0  of  this  report,  has  channel-pair  element  spacings  of  6,  4,  and  id 
half-wavelengths  between  adjacent  channels.  For  the  same  channel-pair  phase  error,  12 
rms,  the  probability  of  ambiguity  is  0.076%-a  600:1  improvement  over  the  three- 

element  array. 

It  is  obvious  that  in  dense  signal  environments,  initial  estimates  of  emitter  location 
(and  later  refinements  of  these  esrimatos)  can  be  accomplished  much  more  rapidly  using 


ROBERT  L.  GOODWIN 

the  more  reliable  estimates  from  the  four-element  array.  It  may  also  be  observed  that  to 
reduce  the  probability  of  ambiguity  of  the  three-element  array  to  0.26%  requires  a re- 
duction in  channel-pair  phase  error  to  approximately  3.0°  ms.  At  the  present  state  of 
the  art  in  microwave  component  technology,  it  is  doubtful  whether  such  channel-pair 
phase  tracking  can  be  obtained  over  more  than  6%  to  10%  bandwidths,  even  with  auto- 
matic calibration. 

The  above  example  clearly  shows  the  advantage  of  implementing  a four-channel  sys- 
tem, if  the  obtainable  performance  can  be  justified  against  the  need  for  (a)  the  additional 
channel  and  (b)  the  additional  processing.  There  are,  of  course,  less  stringent  require- 
ments that  can  be  addressed  very  competently  by  a three-element  array. 

For  example,  suppose  that  it  is  possible  to  achieve  system  channel-pair  phase  track- 
ing to  10°  rms,  and  that  a baseline  length  of  d/X  = 4 (li  = 8 half-wavelengths)  is  available, 
At  boresight,  the  angular  accuracy  of  such  a system  is  0.398°  rms.  Over  the  symmetric 
sector  defined  by  a starting  bearing  of  +15°,  and  an  ending  bearing  of  -16°,  a direction- 
finding system  need  take  only  21  bearing  cuts,  spaced  at  1.5°  increments,  to  attain  a 
0.96  probability  of  determining  the  location  of  an  emitter  to  within  a n n.mi.2  (3.437T  km2) 
area  of  uncertainty  at  a cross-track  range  of  102.66  n.mi.  (190,1  km).  This  particular 
system  has  a probability  of  ambiguity  of  1.01%,  implying  a somewhat  longer  initial 
processing  interval,  on  the  average,  before  “outliers”  in  the  data  could  be  identified  and 
discarded  (in  contrast  to  the  0.076%  pa  of  the  four-element  system). 

For  many  applications,  the  obtainable  performance  in  three-element  systems  is  quite 
satisfactory,  provided  that  (a)  the  location  geometry  is  favorable,  (b)  a sufficient  number 
of  bearing  cuts  can  be  taken,  and  (c)  the  Bystem  channel-pair  phase  tracking  is  good. 

Figure  1-2  is  one  configuration  of  a four-element  interferometer.  The  four  functional 
elements  shown  are  (a)  relatively  nondirective  antennas  for  wide  spatial  coverage,  (b)  phase- 
tracked  receiver  channels  incorporating  hard  limiting  of  channel  signal  levels  to  remove 
amplitude  fluctuations,  (c)  phase  comparators,  and  (d)  an  ambiguity-elimination  and  angle- 
processing circuit.  As  shown  in  the  diagram,  the  channel  at  the  far  left  is  the  phase 
reference.  Other  four-element  array  configurations  are  possible.  Sec.  4.0  will  show  that 
this  particular  array  configuration,  called  “cascaded  end-phase,”  is  the  canonical  configura- 
tion for  a four-element  array.  Other  array  configurations  can  equal,  but  not  exceed,  its 
tolerance  to  angular  ambiguities.  This  configuration  is  optimal  because  of  its  relatively 
simple  processing  compured  with  other  array  configurations. 

For  the  configuration  shown  in  Fig.  1-2,  the  electrical  phase  differences  between  a 
signal  in  the  reference  channel  and  signals  in  the  other  channels  are 

'Plj  = Pd\,j  sin  0a  , (1-2) 

where  is  the  physical  spacing  between  the  phase  centers  of  antennas  1 and  j in  the 
linear  array  (j  = 2,  3,  or  4).  Thus,  three  phase  differences —the  necessary  and  sufficient 
number  to  extract  all  the  AOA-dependent  electrical  phase  information  the  array  can 
provide  —are  made  available  to  the  processing  circuits. 

Suppose  dj  4 - 4 di(3  - 16  dj  2 - 8 wavelengths.  Normalized  to  half-wavelengths, 
the  spacing  integers  (see  Fig.  1-2)  are  p = 1,  p + q - 4,  and  p + q + r = 16;  these 
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Pig.  1-2— Block  diagram  of  four-elermmt  phase  Interferometer 


define  q = 3 und  r - 12.  According  to  Eq.  (1-2),  as  0a  increases,  the  electrical  phase 
angles  associated  with  channel-puir  1 and  4 and  channel-pair  1 and  3 will  eventually  ex- 
ceed |±7r|  radians,  However,  since  each  member  of  the  ensemble  of  phase  comparator 
outputs  can  only  be  known  modulo  27T,  ambiguities— multiple-candidate  AOAs— will  be 
manifested  in  the  phase  comparator  outputs.  On  the  other  hand,  the  phase  comparator 
output  associated  with  channel  pair  1 and  2 can  never  exceed  |±tt|  radians.  This  fact 
can  be  exploited  to  yield  an  ambiguity  resolution  process  that  will  enable  an  estimate  of 
AOA  associated  with  the  correct  mean  AOA  to  be  recovered. 

There  are  l!  - 1 ambiguous  0a  in  an  arruy  whose  overall  length  is  2 half-wavelengths. 
The  probability  that  one  of  these  umbiguous  AOAs  will  be  computed  from  the  set  of 
modulo  2it  phase  angles  applied  to  the  processing  circuits  is  inversely  related  to  the  array 
spacinga  d\j  or  more  appropriately,  to  the  spacing  integers  p,  q,  and  r,  when  four-element 
arrays  are  under  consideration. 


The  Interferometer  synthesis  problem  that  is  the  main  scope  of  this  report  is  to  de- 
termine p and  q in  three-element  arrays;  and  p,  q and  r in  four-element  arrays  so  that  an 
acceptable  balance  is  achieved  between  the  mutually  incompatible  objectives  of  maxi- 
mizing the  accuracy  of  the  AOA  estimate  (increasing  the  overall  array  length)  and  mini- 
mizing the  probability  of  ambiguity  (minimizing  the  overall  array  length)  subject  to  given 
channel-pair  phase  errors. 

The  two  groups  of  readers  to  whom  this  report  is  addressed  are  (a)  electronics  sup- 
port measure  (ESM)  system  designois  whose  direction-finding  requirements  may  be  met 
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by  implementing  arrays  synthesized  according  to  the  principles  given  here,  rnd  (b)  those 
who  may  find  the  general  theory  of  assistance  as  a point  of  departure  in  discovering  the 
principles  of  arrays  of  more  than  four  elements. 

Amy  designs  that  achieve  the  best  obtainable  probability  of  ambiguity  for  given 
overall  array  length  and  channel-pair  phase  errors,  and  are  guaranteed  decodable,  are 
described  here.  Thus,  the  report  may  be  considered  to  be  an  extended  “existence 
theorem”  on  low-ambiguity  arrays.  These  arrays  are  guaranteed  decodable  because  the 
spacing  integers  p,  q,  and  r are  relatively  prime,  but  no  specific  decoding  procedures  are 
given  in  the  report.  The  choice  not  to  include  some  material  on  decoding  was  made 
reluctantly  (in  order  to  stress  ambiguity  aspects  of  array  design),  but  was  done  in  the 
knowledge  that  techniques  based  on  the  Chinese  remainder  theorem  are  well  known. 

Section  2.0  of  this  report  is  a short  review  of  basic  interferometer  theory.  It  pro- 
vides a thorough  exposition  of  just  what  constitutes  an  ambiguity  in  arrays,  and  discusses 
the  ambiguity  constraints  in  arrays  defined  by  two  spacing  integers. 

Section  3.0  begins  with  a discussion  of  the  phase-error  sources  in  multichannel  re- 
ceiving systems,  and  the  magnitudes  of  these  errors  in  current-art  receiver  components. 

Then  an  expression  for  probability  of  ambiguity  based  on  a channel-pair  phase  error 
formulation  is  given  for  three-  and  four-element,  two-integer  set  arrays.  Finally,  the 
optimum  spacinga  for  these  arrays  is  given,  along  with  tabulations  of  pa  for  various  a.iray 
lengths  as  a function  of  channel-puir  phase  error. 

Section  4.0  treats  the  fundamentals  of  four-element,  throo-integer  set  interferometers 
by  first  considering  the  properties  of  various  configurations  of  four  elements.  Resolvable 
and  unresolvable  ambiguities  in  the  two  distinct  three-element  subarrays  that  constitute 
a four-element  array  are  then  discussed.  The  role  of  the  integer  factor  common  to  mem- 
bers of  each  of  the  suburray  ratios  as  this  factor  Influences  overall  array  ambiguity  is 
explored . Then,  explicit  forma  ure  given  for  the  ambiguity  variables  in  each  of  the  two 
subarrays,  and  some  sample  calculations  are  performed  of  probability  of  ambiguity.  Last, 
it  is  shown  that  a particular  form  of  four-element  array,  the  cascaded  end-phase  array,  is 
the  optimum  array  configuration. 

Section  5.0  provides  the  theoretical  basis  for  synthesizing  optimum  arrays  of  various 
longths.  Two  forms  of  array  are  introduced  and  defined  to  achieve  a rapid,  readily  applied 
procedure;  the  “ideul  unrealizable”  array  and  the  “optimum  realizable”  array.  An  approxi- 
mate synthesis  procedure,  ussumtng  independence  of  subarray  ambiguities,  is  given.  For 
most  cases  of  current  practical  concern,  this  approximate  synthesis  procedure  is  very  satis- 
factory. Indeed,  this  procedure  often  provides  several  arrays  of  the  same  overall  length, 
but  with  differing  p,  qt  and  r spacings,  that  achieve  the  same  overall  probability  of 
ambiguity.  However,  by  treating  the  two  subarray  ambiguity  variables  as  members  of  a 
joint  probability  density  function,  with  correlation  between  the  variables,  it  is  possible 
to  derive  an  exact  formulation  for  the  overall  probability  of  ambiguity  in  four-element 
arruys.  The  impetus  for  presenting  this  exact  analysis  is  threefold: 

• The  optimum  array  for  a given  overall  array  length  and  Identical  zero-mean 
channel-puir  phase  errors  can  bo  unequivocally  specified. 
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• The  analysis  can  be  extended  to  arrays  using  antenna  elements  sequenced 
into  a shared  channel  to  achieve  component-usage  economies.  (This  analysis  is  not  given 
in  the  present  report,  but  it  is  a straightforward  extension  of  the  analyses  given  here.) 

• Extensions  of  the  basic  analyses  oriented  toward  other  aspects  of  multi- 
element array  performance  may  be  encouraged. 

A convenient  classification  of  overall  array  length  that  leads  directly  to  a readily  applied 
computer-aided  synthesis  is  then  given.  The  section  closes  with  tabulations  of  probability 
of  ambiguity  for  optimum  arrays  for  a wide  variety  of  array  lengths. 

Extensions  of  the  research  presented  in  this  report  are  in  progress,  and  are  discussed 
briefly  in  a final  section  that  summarizes  the  effort  and  provides  some  concluding  remarks. 


2.0  BASIC  INTERFEROMETER  THEORY 

This  section  presents  the  basic  principles  underlying  multielement  phase  interferom- 
eters. Initial  attention  is  directed  to  the  concept  of  characterizing  arrays  by  two-integer 
sets.  This  generalized  array  concept  emphasizes  that  the  susceptibility  of  interferometers 
to  providing  ambiguous  (actually,  grossly  erroneous)  estimates  of  angle  of  arrival  is  a func- 
tion  of  two  integers  and  two  associated  channel-pair  phase  errors. 

This  method  of  characterizing  an  interferometer  is  compatible  with  either  three- 
element  or  four-element  arrays.  For  the  former,  a portion  of  the  channel-pair  phase 
errors  reside  in  a "common-channel”  phase  error,  whereas  for  the  latter  case,  phase  errors 
in  one  pair  of  channels  are  independent  of  phase  errors  in  the  other  pair  of  channels. 

There  has  been  some  confusion  in  the  literature  concerning  the  apparent  superiority 
of  one  configuration  of  three-element  array  with  respect  to  another  array  configuration 
as  regards  susceptibility  to  ambiguities  for  a given  value  of  channel-pair  phase  error;  as  in 
Kendall  [12]  and  Margerum  [18],  Apparently,  Kendall’s  efforts  in  generalizing  three- 
element.  array  theory  to  embrace  arbitrary  ratios,  i.e.,  arrays  whose  spacings  were  not 
relatively  prime  integers— obscured  the  fact  that  if  one  particular  configuration  of  three- 
element  array  is  ambiguous  (due  to  some  set  of  channel-pair  phase  errors),  then  the  other 
configurations  must  also  be  ambiguous.  The  proof  is  trivial  and  is  given  in  Sec.  2.4. 

Margerum  analyzed  a midphase  three-element  array,  but  neglected  to  consider  the 
effects  of  "common-channel”  phase,  Consequently,  the  ambiguity  constraints  he  derived 
were  actually  those  for  the  four-channel  (independent)  two-integer  set  array.  However, 
Margerum’s  use  of  the  ambiguity-plane  method  of  illustrating  ambiguity  boundary  rela- 
tionships seems  to  predate  its  employment  by  others.  As  will  be  seen  in  Sec.  2.3,  the 
ambiguity  plane  is  a valuable  concept  for  understanding  the  mechanism  underlying  am- 
biguities in  multielement  interferometers. 


2.1  Two-Integer  Set  Interferometer  Fundamentals 

Consider  the  linear  array  shown  in  Fig.  2-1.  where  the  physical  spacing  between 
antennas  1 and  2 is  dL,  and  the  spacing  between  untennus  3 and  4 is  ds.  As  depicted, 
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=$dL  iinO, 
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PRIME 
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Fig,  2-1— One  form  of  two-lr\'qger*ot  interferometer 


spatial  angles  of  arrival  are  defined  with  respect  to  the  normal  to  the  baseline.  With  no 
loss  in  generality  as  far  as  interferometer  ambiguities  are  concerned,  the  sources  of  distant 
signals  are  assumed  to  lie  in  the  horizontal  plane.  It  is  further  assumed  that  these  spatial 
angles  of  arrival  are  restricted  to  |0a|  < 90°  or  that  other  means  have  been  employed  to 
resolve  the  gross  front-toback  ambiguity  inherent  in  phase  interferometers  using  antenna 
elements  with  little  directivity. 

The  electrical  phase  differences  between  the  signals  in  the  two  sets  of  antennas  for 
radio-frequency  (RF)  energy  arriving  from  angle  0O  are 


Elements  1 and  2 


<1 >L  = Ml  sin  °a 


= X dL  9a  , 


(2-la) 


Elements  3 and  4 


*3  = Ms  sin  0a 


IT  dS  sin  0O  . 


(2-lb) 


These  phase  differences  can  be  determined  by  applying  the  antenna  outputs  to  phase- 
tracked  receiver  channels  incorporating  hard  limiters  (to  remove  amplitude  fluctuations) 
and  terminating  in  phase  comparators.  As  is  well  known,  in  an  elemental  two-element 
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phase  interferometer  that  does  not  have  any  available  supplemental  amplitude  informa- 
tion from  the  antennas  for  use  in  resolving  ambiguities,  the  maximum  spacing  between 
antenna  phase  centers  is  restricted  to 

. ff  - (—  rt) 

“max  ~ o- 

Y [sin  (90°)  - sin  (-90°)] 

= (2-2) 


In  interferometers  providing  two  or  more  phase  differences,  it  is  possible  to  resolve 
ambiguities  even  though  one  or  more  spacings  exceed  a half-wavelength  at  the  operating 
frequency. 


Suppose  di  and  d$  in  the  array  of  Fig.  2-1  are  in  the  ratio  of  small  integers,  as 


(2-3) 


where  p and  q are  relatively  prime. 
If 


p complete  cycles  of  27r-rad  phase  change  will  occur  between  the  output  signals  of  an- 
tennas 1 and  2 and  q complete  cycles  of  2?r-rad  phase  change  will  occur  between  the  out- 
put signals  of  antennas  3 and  4 as  0a  varies  between  -90°  and  +90°. 

Electrical  phase  angles  can  be  determined  only  (mod  2ir);  hence,  the  true  phase  dif- 
ferences <1»£,  and  must  be  expressed  as 


d>L  = 'Pl  -pc  (m°d  2tt)  + 2iry 

(2-4a) 

and 

♦ 

4>s  = <2S_pc  (mod  2tt>  + 27tjc  , 

(2-4b) 

where 


subscript  pc  - voltages  available  at  the  output  of  a phase  comparator  defining 
angles, 

x,  y ~ pair  of  integers  which  must  be  determined  so  that  d>£  and  d>5  of 
Eq.  (2-4)  are  equivalent  to  those  of  Eq.  (2-1). 

Determination  of  integers  for  x and  y is  basic  to  any  vernier-resolving  problem  (e.g., 
distance  measurement  applications  iri  which  phase  differences  between  sidetones  displaced 
from  a carrier  are  used  to  resolve  range  ambiguities)  and  is  readily  achieved.  Noting  that 


ROBERT  L.  GOODWIN 


*L  = f *S  . 


and  equating  Eq.  (2-4a)  to  (2-4b)  yield 


ql'fiL- pC  (mod  2it)  + 2iryl  = p^..,  .,*  (mod  2?r)  + 2nx] 


An  equivalent  form  is 


p*s  - wl  + 2*(px  ~ </y)  - 0 . 


where  the  pc  subscripts  and  the  (mod  2it)  notations  have  been  suppressed  for  clarity. 

One  method  of  solving  Eq.  (2-7)  is  to  doubly  iterate  through  the  allowable  Bets  of  x 
and  y until  the  set  giving  equality  is  found.  As  will  be  shown  in  Sec.  3.0,  the  best  strategy 
in  a system  with  phase  errors  is  to  accept  the  x(,  y;-  that  result  in 

\f('PL>'fiS>xbyj)\  < ff  • (2-8) 

2.2  Interferometer  Array  Classification 

Figure  2-2  shows  the  four  array  configurations  capable  of  providing  two  electrical 
phase  differences  for  estimating  spatial  angle  of  arrival.  Only  antennas  and  the  associated 
phase  comparators  are  depicted;  it  is  understood  that  receiving  channels  with  phase  devia- 
tions (as  contrasted  to  ideal  channels)  are  interposed. 

Figure  2-2a,  b,  and  c illustrate  the  three  possible  ways  of  obtaining  two  phase  dif- 
ferences using  three  antennas,  a consequence  of  the  fact  that  the  number  of  combinations 
of  three  elements  taken  two  at  a time  = 31/21(3  - 2)1  = 3.  Fig.  2-2d  shows  an  array  using 
four  elements. 

The  arrays  may  be  classified  according  to  the  channel  used  for  phase  reference.  Thus, 
in  Fig.  2-2u,  channel  1 is  the  reference;  the  an'ay  in  called  End-phase  Left.  Similarly,  in 
Fig.  2-2c,  channel  3 is  the  reference;  the  configuration  is  called  End-phase  Right.  Finally, 
in  Fig.  2-2  b,  channel  2 is  tho  reference;  this  array  is  called  Midphase. 

Obviously,  for  the  array  configuration  in  Fig.  2-2d,  i»  unaffected  by  phase  errors 
in  channels  1 and  2.  Hence,  this  four-channel  interferometer  is  called  Independent. 

Suppose  that  for  each  of  the  three-element  interferometers  of  Fig.  2-2,  the  element- 
to-eloment  spacings  (normalized  to  half  wavelengths  at  the  operating  frequency)  are  p be- 
tween elements  1 and  2,  and  <j  between  elements  2 and  3.  Then,  the  interferometer 
ratios  in  terms  of  electrical  angles  and  spacing  integers  are 

„ , , , (0  f,+s  P + 7 

End-phase  Left  &L  - — - — — — 


1 1 • ' , ■ ...  i ‘.il ‘ - 
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P + q 


(A)  END-PHASE  LEET 
,jO  *1 » S P + q 


(B!  MIDPHASE 

- f 

♦s  q 


tC)  END-PHASE  RIGHT 
jj)  *L»s  P + R 

»s  q 


ID)  INDEPENDENT 
#1  P 


u -*—  n 


q 


Ft*.  2-2— Simplified  block  diagrams  of  the  four  possible 
two-integer  set  interforometurs 


Midphase  - ~ 

ic  P + (7 

End-phase  Right  !RK  = — jp— • ” — q' — 


where 


2 7T 


“ ~ (<*/,  4 ds)&\nOa, 


H = X dL  sin  0o, 
*s  * y dcj  sin  0a . 


Fot  the  four-element  interferometer,  the  analogous  quantities  are 
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Independent 

<H  */,  P 
1 “ d ' 

A complete  classification  should  include  consideration  of  phase-error  parameters  so 
as  to  account  for  unequal  channel  electrical  lengths.  Table  2-1  list*-  the  channel-pair  phase 
error  parameters  appropriate  to  the  four-interferometer  array  configurations,  in  accord 
with  the  notation  of  Fig.  2-2.  Expression  of  deviations  from  some  nominal  channel  phase 
length  in  terms  of  channel-pair  errors  will  be  useful  in  establishing  ambiguity  constraints 
for  the  various  interferometer  configurations.  In  following  sections,  a channel-pair  phase 
error  formulation  will  be  useful  in  deriving  analytic  expressions  for  probability  of  ambi- 
guity in  three-  and  four-element  arrays. 


2.3  Ambiguity  Constraints  for  Two-Integer  Set  Interferometers 

The  constraints  on  channel-pair  phase  errors  that  preclude  erroneous  determination 
of  the  spatial  AOA  (because  of  ambiguities)  will  now  be  derived  for  two-integer  set 
interferometers. 


Consider  un  interferometer  in  which  dL  and  ds  are  in  the  ratio  m:n.  The  integers 
m and  n are  used  to  emphasize  that  a generalized  interferometer,  rather  than  one  of  the 
specific  four  configurations  introduced  earlier,  is  being  discussed.  The  maximum  allow- 
able spacings  for  the  rutio  $ = min  are,  of  course,  dL  ~ mX/2,  and  ds  - nX/2.  Otherwise, 
more  than  one  0a  within  the  -90°  to  +90°  field  of  view  (FOV)  will  produce  a specific  (mod 
2n)  <I>£,  <I>s  set.  If  m - 3,  n - 2,  dfj  - (3/2)X,  ds  =>  X,  then 


*L-  max 

= 1 

- max 

“I 

• ~-X  = ±3tt  - ±mir 

• X - ±2tt  = ±mtt|  . 


If  it  were  possible  to  measure  <I>^  and  *I>S  unambiguously  over  their  ranges,  a phnse- 
plune  plot  of  »l>£  vs  would  appear  as  shown  in  Fig.  2-3a,  in  a form  due  to  Marge  rum 
[13].  A more  convenient  representation  is  Fig.  2-3b,  which  is  centered  on  0a  - ±90° 
rather  than  on  0a  ~ 0°.  The  expressions  for  <l>£  and  are  thus  modified  to 


‘I’i,  “ PdL  sin  0a 
«l'.v  = Pds  sin  0a 


< 0 < 90", 


(2-9a) 


and 


4’i,  =•  Pd^(2  + sin  0a) 
<l\s  = Pds(2  + sin  0o) 


> -90°  < 0a  < 0°. 


(2-9b) 
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1 to  4 are  deviations  of  channel  f elect 
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(A)  CENTERED  ON  6,-0" 


(B)  CENTERED  ON*,  - ±80° 


Flu,  2-3— I’haso-plane  diagram*  for  m:n  « 3:2 


Sinco  and  'l»Jj  can  be  known  only  to  within  2tt  rad  at  the  outputs  of  phase  com- 
parators, the  phase-plane  plot  must  be  collapsed  to  a square  2ir  rad  on  a side,  as  shown 
in  Fig.  2-4.  In  this  diagram,  the  intersections  of  the  <l>fL  vs  ‘hjj  trajectories  on  the  <pL  axis 
are  spaced  2n/n  (=  7r).  The  <ps  axis  intersections  are  spaced  2ir!m  (=27r/3). 


Fig.  2-4  — Ambiguity  diagram  for  m:n  « 3:2 
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Suppose  0a  --  45°.  Then 

4>i,  - 2ir  X 3/2  X lh/2  , 
and 


y'L  » (mod  2 ff)  « 2tr  (3/2v/2  - 1) . 

Also, 

4>s  = 2ir  xix  l/s/5  j <fis  * 2ff/v/2 . 


From  Eq.  (2-7),  we  have 


+ 2w (rrwe  - ny)  = 0 . 

Substituting  for  ^ and  s5s  and  dividing  by  2ff  yields  a trajectory-establishing  equa- 
tion of  the  form 

3#  + 2 ■=  2y  . (2-10) 

The  only  solution  of  Eq.  (2-10)  subject  to  constraints  0 < x < (n  - 1)  and  0 < y < (m  - 1) 
is  jc  = 0,  y = 1.  (Note:  The  Fig.  2-3b  formulation  of  the  ambiguity  plane  leads  only  to 
positive ’solutions  of  Eq.  (2-10)  and  is  much  less  cumbersome  than  an  approach  based  on 
Fig.  2-3a.)  Thus, 

*‘L  = H>L  + 
and 

+ 2?r(0) 

2ff 

The  spatial  angle  of  arrival  0a  can  be  recovered,  for  example,  from  the  expression 


r i 

=•  sin-1 

a 

1*1.  -max. 

L 27rXi  J 

Errors  in  the  channel-pair  signals  sent  to  the  phase  comparators  will  move  the 
<fls  set  off  the  phase  trajectories  (solid  lines  0,  la,  lb,  or  2)  in  Fig.  2-4.  For  simplicity, 
the  prime  notation  on  4’^,  d’g  and  y?/,,  will  be  dropped  henceforth. 
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If  any  <!>£,  <I>S  3et  is  in  error  because  of  channel  deviations,  then  the  (mod  2tr)  repre- 
sentation at  the  phase-comparator  output  is  also  in  error.  The  actual  phase  differences 
(including  errors)  compared  to  the  phase  differences  made  available  by  the  instrumenta- 
tion are  as  follows: 


Actual  Electrical  Phase  Differences 


2tt 

" X dl  Sin  00  + A<^  “ + 

(2-lla) 

*a„t  " X ds  sin  6°  + A*s  " + A{fis  • 

(2-llb) 

Available  Electrical  Phase  Differences 

<Pr,,v  ■ + A^z.](mod2ir) 

(2-1 2a) 

^sav  “ l+s  + Av?sl  (mod  2ir) . 

(2-l2b) 

The  geometry  of  Fig.  2-4  shows  that  to  equalize  the  possibility  of  incurring  ambi- 
guities over  the  -90°  to  +90°  FOV,  the  dashed  boundaries  separating  the  region  around 
one  trajectory  from  the  region  around  another  should  be  located  parallel  to,  and  equi- 
distant from,  adjacent  phase  trajectories.  The  ambiguity  boundary  directly  above  tiajec- 
tory  ‘O’  intersects  the  axis  at 

aV>x,-a/  “ TP 

(2-13a) 

and  the  boundary  directly  below  trajectory  ‘O’  intersects  the  <ps  axis  at 

&PS-AI  ==  W' 

(2-13b) 

In  attempting  to  generalize  the  conditions  for  obtaining  minimum  tendency  toward 
ambiguities  over  the  field  of  view,  Kendall  [12]  considered  ratios  ft  = min  (or  n/m)  in 
which  m and  n were  not  restricted  to  relatively  prime  integers.  He  employed  a need- 
lessly complex  analytical  formulation,  whereas,  with  just  the  previously  used  geometric 
arguments,  it  is  easy  to  show  that  m and  n must  be  relatively  prime  integers,  as  in  the 
following  example. 

Figure  2-5  shows  a phuse-plane  plot  of  the  trajectories  for  ft  - y/Tl/2  - 3.3166/2, 
subject  to  the  element  spacing  associated  with  m equal  to  a/TI/2  wavelengths.  The  solid 
trajectories  in  the  figure  thus  terminate  with  3a,  at  a Vl  value  of  0.3166  X 27r.  ThiB 
indicates,  of  course,  that  over  a 180u  range  on  0a,  <I»j,  sweeps  through  3.3166  X 2n  rad, 
and  ‘hg  sweeps  through  2ir  rad.  Since  ft  is  irrational,  one  can,  in  theory,  extend  the 
array  length  Indefinitely  without  encountering  a condition  whore  multiple  0a  give  rise  to 
the  same  yL,  <ps  sot.  As  more  trajectories  are  added  to  Fig.  2-5  (implying  increased  array 
length),  the  spacing  between  trajectories,  and  hence,  the  tolerance  to  ambiguities,  decreases. 
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Of  the  additional  trajectories  3b  and  4 (shown  dashed  In  Fig.  2-5),  4 terminates  near 
the  point  “ 2ir,  2n.  This  might  have  been  anticipated  by  noting  that  the  integer- 

set  ratio  nearest  id  = ^/Tl/2  it  Id  = 5:3.  For  the  latter  ratio,  there  would  be  five  equally 
spaced  trajectory  intersections  on  the  horizontal  axis  of  the  phase  plane  and  three  inter- 
sections on  tho  vertical  axis.  Because  -^/TI/2  is  irrational,  however,  there  is  no  way  of 
drawing  ambiguity  boundaries  to  equalize  their  spacing  from  adjacent  trajectories,  theroby 
making  the  tendency  toward  ambiguities  implicitly  Independent  of  0a.  Consequently, 
ratios  formed  by  relatively  prime  integers  are  preferable. 

All  phase  trajectories  and  ambiguity  boundaries  in  phase-plane  plots  intersect  the 
i ps  axis  at  an  angle  7 = tan-1  [(2ff/n)/(2ff/m)l  = tan"1  [(m/n)  = ill . For  any  0a,  it  is  easy 
to  show  that  tho  two  ambiguity  constraints  are 


AiPl  > KAi ps  + Ai£>£_Ai  , 


(2-14a) 


I 


< idAv>S  - ^L- AI> 


(2-14b) 


for  !d  = min,  with  m and  n relatively  prime  Integers. 

If  the  upper  inequality  is  satisfied,  the  ambiguity  boundary  above  a given  trajectory 
has  been  crossed;  an  analogous  statement  holds  for  the  lower  inequality. 

An  expression  equivalent  to  Eq.  (2-14a,  b)  which  makes  use  of  fd  ~ min  and 
A = rt /n  is 


InA^  - mAi/>s|  s*  Jr , 


(2-15) 


where  > Implies  an  ambiguity,  and  < implies  rio  ambiguity. 
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Equation  (2-15)  is  the  ambiguity  constraint  for  any  two-integer  set  interferometer, 
provided  that  the  channel-pair  phase  error  terms  A *pL  and  Aips  are  defined  properly  with 
respect  to  the  particular  array  configuration  being  analyzed. 


2.4  Ambiguity  Constraints  for  Specific  Array  Configurations 

The  ambiguity  constraints  for  the  four  possible  two-integer  set  interferometer  con- 
figurations are  now  obtained  by  making  use  of  the  constraint  just  presented  for  a gen- 
eralized Interferometer. 

Equation  (2-15)  is  expressed  in  terms  of  generalized  array  integers  m and  n and 
channel-pair  error  parameters  A^  and  Atpg.  The  array  ratios  (from  Sec.  2.2)  and  the 
channel-pair  eitor  parameters  (from  Table  2-1)  for  each  of  the  array  configurations  are 

End-phase  Left 


m P 


co  = — 

n p 


Midphase 


<d  = m.  P . 

T ?' 


End-phase  Right 


m 


“ 7T  " ' q 


Independent 


(D  _ J22  _ P , 

K ' ~ n " q ' 


A <PL  = 

Aipx  - 

A(p3 

A< fs  ~ 

Aip!  - 

A<p  2 

Ai fiL  = 

~Aip2  4 

' Aip^ 

> 

£ 

n 

A<p2  " 

Aip3 

A Pi  = 

A<P3  - 

AiPi 

Av>s  “ 

Alps  - 

Aip2 

AiPl  != 

A<pi  - 

Aip2 

AiPs  = 

AiPs  - 

Aip4 

(2-15) 

yields 

End-phase  Left 


|p(Av?i  - Aip3)  - (p  +q)(Acp1  - A<p2)l  £ it, 


(2-16a 


(2-16b 


(2-16c 


(2-16d 


(2-17a 


Midphase 


|q(-Ap2  + Ay?i)  - p(Ay?2~Ap3)|  £ it, 


(2-17b 
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End-phase  Right 

\q(A*p3-  - (p  + q)(  A<p3  - Aip2)l  £ « , (2-17c) 

Independent 

\q(&Vl  ~ A<p2)  - p(Avj3  - Av?4 )i  £ ff  . (2-17d) 

It  is  easily  shown  that  Eq,  (2-1 7a,  2-'17b,  and  2-17c)  are  all  equivalent  to  the  follow- 
ing equation. 

Three-element  Interferometer 
Ambiquity  Constraint 

iqAifi!  - (p  + q)Aif>2  + pAip3|  ^ it.  (2-18) 

Equation  (2-18)  is,  of  course,  a consequence  of  the  fact  that  the  channel-pair  error  param- 
eters in  a three-element  interferometer  are  related  by 

a + c = b , (2-19) 

where 

O = Ay?!  - Av?2,  on  <fiL, 

c - Ap2  ~ AlPg,  on  < PS, 

b = Apx  - Av?3,  on  <pL+s, 

The  results  of  Sec.  2.0  are  the  basis  for  correcting  past  misconceptions  concerning 
the  tendency  toward  ambiguities  for  various  configurations  of  three-element  interferom- 
eters. Kendall  [12]  argued  that  the  end-phase  configuration  is  superior  to  the  midphase 
configuration.  Margerum  [18]  analyzed  only  the  midphase  configuration,  giving  actually 
an  expression  for  the  probability  of  ambiguity  for  the  independent  four-element  con- 
figuration. This  left  the  impression,  perhaps  by  omission,  that  tho  end-phase  array  con- 
figuration is  inferior  to  the  midphase  configuration. 

The  correct  statement  for  the  ambiguity  constraints  in  three-element  interferometers, 
based  on  tho  development  given  here,  is  the  following:  Any  configuration  of  three- 
element  interferometer  defined  by  spacing  integers  p and  q (implying  identical  overall 
array  length,  regardless  of  configuration)  has  the  same  tendency  toward  ambiguity,  re- 
gardless of  the  channel  employed  as  the  phase  reference. 


3.0  PERFORMANCE  OF  TWO-INTEGElt  SET  INTERFEROMETERS 

For  many  applications,  the  performance  obtainable  with  two-integer  set  interferom- 
eters, either  three-  or  four-element,  is  adequate.  For  example,  a three-element  array 
8 half-wavelengths  long  at  the  frequency  of  operation  with  channel-pair  phase  errors  of 
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10  electrical  degrees  rms  can  achieve  a boresight  angular  accuracy  of  0.40  spatial  degrees, 
and  a probability  of  ambiguity  of  1.01%.  Performance  such  as  this  can  satisfy  those 
requirements  which  do  not  need  the  “super”  accuracy  obtainable  with  multielement 
arrays  that  have  longer  overall  baselines.  The  primary  reason  for  presenting  a compre- 
hensive theory  of  two-integer  set  arrays  is  that  they  are  the  basic  subarrays  that  can  be 
cascaded  to  form  longer  arrays  of  four  or  more  elements.  In  Secs.  4 and  5,  examples  of 
the  dramatically  improved  performance  that  can  be  achieved  with  these  arrays  will  be 
given. 


■1 

I 

J 

I 
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This  section  discusses  error  sources,  carefully  distinguishing  between  a “channel” 
error  description  and  a more  convenient  formulation  for  design  and  analysis  purposes— 
the  “channel-pair”  description  of  error  sources.  -1 

A 

A channel-pair  error  correlation  coefficient  is  introduced  and  defined  for  the  end-  t 

phase,  midphase  and  independent  array  configurations.  Although  it  is  not  mandatory  to  I 

use  a channel-pair  error  and  correlation  coefficient  formulation,  this  approach  to  deriva-  • 

tion  of  probability  of  ambiguity  facilitates  the  development  of  an  exact  expression  for  ’ 

it  in  four-element  threo-integer  set  arrays  in  Sec.  6. 


A recurrence  relationship  for  generating  the  optimum  spacings  for  three-element 
arrays  of  arbitrary  length  is  given  next.  It  is  Interesting  that  the  full  tabulation  of  allow- 
able spacings,  I.e.,  nonredundant,  for  arrays  of  various  lengths  haB  a counterpart  in  the 
Farcy  sequences  of  rational  fractions  from  number  theory.  Travers  [14]  points  out,  how- 
ever, that  propagation  anomalies  in  the  HF  band  (3  to  30  MHz)  may  constrain  the  “small- 
spuced  channel-pair"  spacing  to  one  half-wavelength  rather  than  the  optimum  “Bmall-spacod 
channel-pair”  spacing  for  the  overall  array  being  used,  resulting  in  fewer  ambiguities. 

The  section  concludes  with  a tabulation  of  probability  of  ambiguity  vs  array  length 
(with  channel-pair  phase  error  as  a parameter)  for  both  three-element  (end-phase  or  mid- 
phase)  and  four-element  (independent)  two-integer  set  array  configurations.  Historically, 
two-integer  set  interferometer  arrays  seem  to  have  been  Implemented  first  in  the  inde- 
pendent configuration  as  in  Bailey  and  Moller  [16].  Watters,  Rees,  and  Enstrom,  [16] 
have  reported  on  a two-frequency  technique  equivalent  to  tho  independent  configuration. 
Later,  as  the  theory  of  arrays  improved  und  the  component  art  advanced,  the  three- 
element,  two-integer  sot  arrays  became  much  moro  prevalent.  A recent  example  of  a 
commercially  oriented  three-element  array  is  reported  by  Watanabe,  et  al.  [17] . Many 
current  military  surveillance  systems  also  employ  the  three-element  array. 


3.1  Interferometer  Channel-Pair  Errors 

In  an  actual  interferometer,  the  electrical  length  i/J;  of  u given  channel  muy  be  many 
thousands  of  electrical  degrees.  This  length  may  vary  because  of  changes  in  operating 
frequency,  as  a function  of  temperature,  or  with  input  signal  power  level  (AM-to-PM  con- 
version). The  length  y3(  can  be  represented  (as  in  See.  2.0)  by 


Vi  “ ^nominal  + <V( 


- 


where  Ay?(  is  the  deviation  of  channel 


in  phase  length  from 


(3-1) 


-y- ■ — 


► 
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Individual  channel  errors  A ^ are  accessible  only  by  measurement  to  some  standard. 
However,  channel-pair  errors 

Ay>lj  s <pt  - <pj 

~ (^nominal  + " ('/’nominal  + 

= A (/),•  - Aip;  (3-2) 

are  manifested  in  the  outputs  of  phase  comparator's  connected  between  pairs  of  channels. 

It  is  convenient  to  define  interferometer  performance  in  terms  of  channel-pair  phase 
errors,  rather  than  channel  errors,  for  the  following  reason.  The  individual  components 
that  constitute  u channel  are  usually  more  readily  specified  in  terms  of  their  deviations 
with  respect  to  a corresponding  component  in  another  channel,  rather  than  in  terms  of 
their  deviations  from  nominal  (but  not  conveniently  measured)  electrical  lengths. 

The  variance  of  channel-pair  phase  error  between  (th  and  >th  channels  for  the  feth 
component  in  an  n-eloment  cuscade  is 

ah,k  = t’Mfc}  + E{^U) 

“ 4k  + 4k  (3"3> 

assuming  thut  A« Aipy are  Independent  and  have  aero-mean  deviation  from  'Pnomlnul-fc 
over  applicable  parameters  (RF  frequency,  power  levels,  supply  voltage,  angle  of  arrival, 
etc.). 


The  overall  channel-pair  error  vuriunce,  summing  over  n components,  is  then 

n 

4j  “ E abk  = 4 + 4 • <3-4) 

k- \ 

For  a large  number  of  cascaded  components,  all  of  whose  variances  are  comparuble, 
tho  channel-pair  phase  error  distribution  tends  toward  Gaussian  (central-limit  theorem). 
Statistical  analyses  of  channel-pair  errors  obtained  on  actuul  multichannel  systems  support 
this  contention. 

Figure  3-1  is  u block  diagram  of  one  channel  of  a multichannel  interferometer  show- 
ing the  principal  contributors  to  channel-pair  electrical  phase  error.  The  direct  output 
from  phase  comparators  is  unalog.  That  is,  <1>(  i(0a)  is  defined  by  sin“l(Ujj.#/i»mBX)  and 
cos"1  (Vj  j.clv mux),  where  vtj.,  and  vtj.c  arc  voltages  from  the  sine  and  cosine  phase  de- 
tectors, respectively,  of  the’ phase  comparator  and  umilx  is  the  maximum  phase  detector 
output  voltage. 

In  automatic  systems,  it  is  convenient  to  perform  ambiguity  elimination  and  other 
calculations  digitally.  Hence,  additional  contributions  to  channel-pair  phase  error  are 
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ANTENNA  GROUP 


RECEIVER  GROUP  PHASE-COMPARATOR  GROUP 


Ai^ld  due  to  multivoltage  level  detector  nonlinearities,  and  due  to  the  quantization  of 
electrical  phase  over  2tt  rad.  Detailed  consideration  of  the  effects  of  AipLD  and  A«Pq  on 
probability  of  ambiguity  is  beyond  the  scope  of  this  report. 

Table  3-1  lists  typical  values  of  the  major  contributors  to  channel-pair  phase  error. 

It  is  assumed  that  the  individual  channel-pair  errors  are  uniiormly  distributed  between 
the  limits  shown.  Phase  errors  due  to  finite  signal-to-noise  ratio  (SNR)  are  not  specifically 
included,  but  can  be  readily  added  with  the  aid  of  the  relation  = 180°/7r>/SNR,  where 
SNR  is  in  termB  of  power.  The  associated  standard  deviation  of  channel-pair  phase  error 
is  I A<Px-nux  \l>/$  where  Ay?x_max  is  the  peak  value  of  phase  error  for  component-pair  X. 
Two  values  of  phase  “noise”  due  to  quantizing  are  given.  These  represent  upper  and 
lower  bounds  on  the  degree  of  quantizing  typioully  employed  in  interferometer  systems 
of  the  type  described  in  Sec.  1.0. 

Table  3-1  shows  that  the  phase  error  due  to  the  quantizer  group  is  more  than  half 
the  overall  phase  error  on  a root-sum-square  basis  for  4-bit  quantizing.  On  the  other 
hand,  quantizing-phase  noise  is  negligible  for  7-bit  encoding.  In  general,  higher-speed  de- 
coding (and  possibly  less  complexity  in  the  ambiguity  algorithms)  can  be  achieved  with 
low  degrees  of  quantizing.  Conversely,  systems  employing  higher  degrees  of  quantizing 
perform  closer  to  the  theoretical  probability  of  ambiguity  for  an  analog  processing  system. 

Radio-frequency  calibration  can  be  employed  in  interferometer  systems  to  reduce  the 
magnitude  of  channel-pair  piiase  errors.  Alternatively,  by  the  use  of  calibration  techniques, 
low-quality  components  (poor  phase  tracking)  can  provide  performance  comparable  to  that 
achieved  in  uncalibrated  systems  using  high-quality  components.  Calibration  signals  are 
usually  introduced  into  the  channels  directly  behind  the  antennas. 


* 
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Table  3-1  —One-Sigma  Channel-Pair  Phase  Error* 
No  Calibration 


Group 


Antenna  Group: 
Antenna 

Transmission  line 


Limits  on  A j+  Standard  Deviation 


RF  filter 


Receiver  Group: 

RF  mixer  and  LO 
power  divider 

IF  preamplifier 

IF  filter 


Phase  Comparator  Group : 
IF  limiter 
Phase  comparator 


Quantizer  Group: 
Level  detector 
Quantizer,  7 bit 
Quantizer  (4  bit) 


Overall  channel-pair  phase  error  = 


9.848°  rms 
(11.769°  rms) 


•"Strong-signal"  conditions.  ^Uniform  distribution  assumed. 
Note:  Present  technology,  wideband  systems. 
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The  use  of  calibration  signals  whose  frequencies  are  equal  to  signals  of  interest  makes 
possible,  in  principle,  the  removal  of  channel-pair  phase  errors  due  to  RF  and  IF  filters, 
RF  mixer  and  local  oscillator  power  divider,  IF  preamplifier,  and  phase  comparator  (see 
Fig.  34).  Since  the  power  level  of  an  incoming  signal  within  tho  system  dynamic  range 
is  arbitrary  relative  to  a fixed-level  calibration  signal,  only  a portion  of  the  channel-pair 
phase  error  from  the  IF  limiters  can  be  removed.  AIbo,  the  signal  angle  of  arrival  is 
arbitrary  with  respect  to  the  fixed  angle  of  arrival  synthetically  introduced  into  the  sys- 
tem during  the  calibration  mode.  Thus,  quantizing  group  errors  must  be  accounted  for 
twice:  First,  during  the  signal  reception  mode;  and  Becond,  during  the  calibration  mode. 

The  list  of  error  sources  and  their  magnitudes,  typical  of  a system  employing  calibra- 
tion, is  given  in  Table  3-2.  The  independence  of  these  errors  during  the  two  modes  of 
operation  is  assumed. 


Table  3-2— One-Sigma  Channel-Pair  Phase  Error* 
With  RF  Calibration 


Component 

Limits  on  At p,-  jt 

Standard  Deviation 

ANALOG  PROCESSING 

Antennu 

O 

00 

-H 

°ANT  _ 4-619° 

RF  calibration  network 

= ±3° 

oCAL  = 1.732° 

IF  limiter 

I! 

H“ 

05 

c 

°LIM  - 1-732° 

Root-Sum-Square  - 6.228° 

DIGITAL  PROCESSING 

Antenna 

= ±8° 

0ANT  ” 4.619 

RF  calibration  network 

» +3° 

°CAL  = 1-732° 

IF  limiter 

= ±3° 

°LIM  - 1-732° 

Level  detector  (signal) 

» ±1° 

0LD-SIG  = 0.677° 

Level  detector  (calibration) 

= ±1° 

°LD-CAL  = 0.677° 

Quantizer,  7-bit  (signal) 

- ±1.406° 

ffQ  SlG  ” 0.812° 

Quantizer,  7-bit  (calibration) 

= ±1,406° 

°Q-CAL  = 0-812° 

Root-Sum-Square  » 5.416° 

♦“Strong-signal"  conditions. 

'Uniform  distribution  assumed. 

Note:  Present  technology,  wideband  systems. 
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3.2  Probability  of  Ambiquity  for  Two-Integer  Set  Interferometers 

The  ambiguity  constraints  for  the  three-element  interferometer  configurations  of 
Sec.  2.4  are 


End-phase  Left 

|p(Asf>3  - AyJ3)  - (p  + q)(Ay5,  - A^2)|  £ rr 

where 

Ai/?i  - Av?3  3 Ay>L , Aip!  - Ayj2  h Ai ps 

,R  _ jn  _ Pj_Q , 

" T p 

Midphase 

lq(-Aip2  + Ay^)  - p(Ay?2  - Ay?3)l  ^ rr 

where 

~Ay?2  + A^t  3 A V?L,  Ay>2  - Aipg  3 A(ps 


(3-Ba) 


(3-5b) 


io  _ m _ P 
“ T~~q' 

End-phase  Right 


l<lOV3  - A v>i)  - (p  + <z)(Ay>3  - A(p2)|  ^ ?r  (3-Bc) 

where 

3 AipL  , Ay’s  - Atp2  = 

P +_Jf 

~~q 

it  was  shown  that  Eq.  (2-17a,  b,  and  c),  repeated  above  for  convenience, 
a single  ambiguity  constraint 

I q Ay5|  -•  (p+q)Ay>2  + pAip3|  \ n . (3-5d) 

In  this  section,  an  expression  for  the  probability  of  ambiguity  of  a generalized  two- 
integer  set  interferometer  will  be  derived,  starting  from  Eq.  (3-5a,  b,  and  c)  rather  than 
from  Eq.  (3-5d),  as  could  be  readily  accomplished.  The  motive  for  this  indirect  deriva- 
tion is  to  focus  attention  on  channel-pair  errors  and  channel-pair  error  correlation  coef- 
ficients. A formulation  of  the  probability  of  ambiguity  in  terms  of  these  parameters 


Ay>3  - A <p-L 


m 

n 


In  Sec.  2.4, 
are  equivalent  to 
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which  simplifies  the  derivation  of  an  exact  expression  for  the  probability  of  ambiguity  in 
four-element  interferometers  will  be  seen  later  (in  Sec.  B). 


Expressions  for  the  channel-pair  error  correlation  coefficient  will  be  required.  For 
two  zero-mean  random  variables  fl  and  18,  the  correlation  coefficient  between  them  is  de- 
fined [18]  by 


P 3 


Em) 


°A°B 


(3-6) 


where  J5[x]  * expected  value  of  x. 

From  Eq.  (3-Ba),  for  the  end-phase  left  configuration,  is  derived 


Pe 


_ - Ay32)? 


x/o?  + ojj  y°?  + °2 


We 

- +0.B , all  o,  equal. 


(3-7) 


Measurements  on  the  Joint  statistics  of  aL  and  as  iure  necessary  to  define  pe\ separate 
measurements  of  oL  and  as  in  Eq.  (8-7)  are  insufficient.  For  array  design  purposes,  the 
assumption  that 


of 


f « of  « of  - 4 of  ° 4-  I’l 


is  a reasonable  one. 


It  is  readily  shown  that  pe  = +0.6,  all  ai  equal,  for  the  end-phase  right  configuration 
us  well. 


From  Eq.  (3-6b),  for  the  midphase  configuration,  is  derived 

_ g{(-A<p2  + A^t)(A^g  - Ay?3)} 


Pm  s 


y£f  + s/o§  + o ] 


-at 


°L°S 

-0.6,  all  equal. 


(3-8) 
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I'he  variance  of  the  left  wide  of  Eq.  (3-oa),  the  end-phase  left  configuration,  is 
°e  ~ P2(°l  + °1)  " 2 p(p  + - Av?3)(Av?i  - A« p2)} 

+(p  + q)2(of +o|).  (3-9) 

Substituting  p = n,  p + q = m,  of  + o§  = of,,  and  of  + o § = o$  yields 

02  ■ m2o§  - 2mn£{(Av5x  - AipaKA^  - Aip2)}  + n2ofi . (3-10) 

Since,  from  Eq.  (3-7) 

£{(Av?!  - Av?3)(A^!  - A^2)}  = or?  = Pe°L°S  . 

Eq,  (3-10)  becomes,  for  the  end-phase  configuration, 

o2  ■ m2o§  - 2pemnaLas  + n2o f,  , (3-11) 

where 

m = p + q,  Ol  = channel-pair  phase  error  standard  deviation,  channels  1 and  3 
n = P,  os  - channel-pair  phase  error  standard  deviation,  channels  1 and  2 

Pe  ~ aV°L °s 

= +0.6,  all  Oj  equal. 

The  variance  of  the  loft  side  of  Eq.  (3-6b),  the  midphase  configuration,  is 

o2  = q2(o|  + of)  -•  2pqfi{(-A»p2  + A^!)(A^2  ~ Av>a)>  + P2(o$+o§).  (3-12) 
Substituting  p - m,  q = n,  of  + o$  = of,,  and  o§  + o2  = o§  yields 

o^,  = m2oj;  --  2mnE{(-Atp2  + A<pj )(A<p2  " A<p2)}  + n2oL  . (3-13) 

Since,  from  Eq.  (3-8),  we  have 

E{(-A «p2  + Aipi)(Aip2  - A*s)}  = -o|  = pmoLo6.  > 

Eq.  (3-13)  becomes,  for  the  midphase  configuration, 

o'fn  =*  m2o|  - 2pmmnoLos  + nVf,  , (3-14) 

where 
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m ~ P>  °h  c channel-pair  phase  error  standard  deviation,  channels  1 and  2 
n = q,  os  = channel-pair  phase  error  standard  deviation,  channels  2 and  3 


Pm  ~ 


- -0.5,  all  equal. 

Equations  (3-11)  and  (3-14)  have  the  same  functional  form;  hence  a single  relation 
suffices  provided  m,  n,  olj(  os,  and  p(pe  or  pm)  are  properly  defined.  In  fact,  if  we 
refer  to  Eq.  (2-17d),  Sec.  2.4,  for  the  four-element  independent  interferometer,  it  is  ap- 
parent that 


£{(Av>a  - A<p2)(A*p3  - Av?4)} 

P(  , — 0 


(3-15a) 


" o?  + =■  o§  + of 


(3-  15b) 


Thus,  the  variance  of  the  ambiguity  variable  is 


- m2ojj  - 2pmnoLos  + n'Vjf,  , 


(3-16) 


where 


p » +0.5,  end-phase  ") 

= 0 , independent > all  ot  equal. 

-•  - 0.5,  midphase 

An  alternate  expression  for  </$  is 


= aSn*PRl(m5«)2  - 2pmXnA  + (nA)2], 


(3-17) 


where 

°CH-PR  ~ 11  nominal,  or  design,  valuo  for  channel-pair  rms  phase  error 
A s or,/aCH.PR 
2 tfs/aCH-PR- 

For  design  purposes,  ull  ot  am  usuully  specified  equal,  implying  thut  = os  - ^CH-PR  " 
yj’l  Of.  Thus  X = A = unity,  und  a simplified  form  of  o$  is 

al  rj  o&h-PrI™2  ~ 2pmn  + n2]  . (3-18) 
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The  probability  of  correct  resolution,  that  is,  the  probability  that  the  left  sides  of 
Eq.  (3-5a,  b,  or  c)  will  not  exceed  rad,  is 


(3-19) 


Setting  t'  = y'2- 1,  dt!  = y/Z  dt  and  changing  the  limits  of  integration  gives 

2 v/r/v/J 


2 r 

= "'(Tftj) 


exp  (-t2) dt 


(3-20) 


Thus,  the  probability  of  ambiguity  (the  probability  that  the  left  sides  of  Eqs.  (3-6a,  b, 
or  c)  will  exceed  ±ir  rad)  is 


where 


Pa  = 1 “ Pc 

- 


(3-21) 


oty  = °CH-PR[(mi)2  - 2pm2nA  + (nA)2]1/2  , 

m - integer  associated  with  large  spacing  (m  > n), 
n <=  integer  associated  with  small  spacing, 

,;L<  °S  ~ st«ndard  deviations  of  channel-pair  phase  errors  for  large  and  small  spacings, 
respectively, 

p - £{S?S}  /o^Of;,  correlation  coefficient  between  large  and  small  channel-pair 
phase  errors, 

och-pr  “ design,  or  nominal,  value  for  channel-pair  phase  error, 

2 = O^/Oqh-PRi  A = OfJO  CH-PR’ 


An  asymptotic  expression  for  erfc  (f)  is  [19] 


erfc  (f) 


(T H 
ty/v 


1 + 1-3 


l-3«6 


2f2  (2<2)2  (2<2)3 


+ 4 , , 


(3-22) 


If  the  first  two  terms  above  are  utilized,  pa  is  bounded  by 
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(3-23) 


3.3  Optimum  Amy  Ratios  in  Three-Element  Interferometers 

The  argument  of  the  complementary  error  function  defining  probability  of  ambiguity 
Eq,  (3-21),  with  the  simplified  form  of  o^,  is 


7T 

0CH-PR  s/mi  ~ 2 pmn  + n2 


For  three-element  arrays  of  length  8 = p + q (in  half-wavelengths  at  the  operating  fre- 
quency), legardless  of  which  element  is  phase  reference,  the  quadratic  form  within  the 
radical  in  the  denominator  of  the  expression  above  can  be  written  as 

P2  + pq  + q2  , p > q 


where 

p - integer  related  to  the  channel  1— channel  2 baseline, 
q = integer  related  to  the  channel  2— channel  3 baseline. 

To  determine  which  sets  of  integers  p,  q are  optimum— in  the  sense  of  providing 
arrayN  having  the  lowest  probability  of  ambiguity  for  given  8— subject  to  the  constraint 
p + q - constant  = 8,  it  suffices  to  determine 

[p2  + P(S  - p)  + (8  - p)2]  = 0 , 


which  reduces  to 

JL  [p2  - 8p  + 82)  = 0. 

For  no  restrictions  on  p and  q,  the  solution  iB,  of  course,  p = q - 8/2.  To  form  an  allow- 
able interferometer  ratio  = p/q  or  q/p,  however,  the  two  parameters  must  be  relatively 
prime  integers. 


Thus  p “ q (=1)  only  for 

8 — 2 ; p - 1, 

= 3 ; =2, 

- 4;  = 3 , 


8=2.  The  optimum  p,  q for  8 =*  2,  3,  4 are 

q = 1 
= 1 

= 1 (p  = q = 2 is  not  allowable  since  it  is 
a degenerate  form  of  p = q = 1). 


For  8=5  through  10,  the  optimum  sets  of  p,  q are 


I 

i 


1 


■1 


1 s 
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r 


5; 

P = 

3, 

q = 

2 

6 ; 

= 

6, 

= 

1 

7; 

= 

4, 

- 

3 

8; 

= 

6, 

= 

3 

9; 

ES 

5, 

= 

4 

10; 

- 

7, 

= 

3 

Functions  generated  by  induction  for  p and  q,  given  6 = p + <j,  are 


P/,-1  53  21, 

fl/,-1  = 2/  - 1 ; 

B/,-1  “ " 1 

P/,0  “ 2f  + 1, 

ql  0 « 2i  - 1 ; 

fi/,0  3 41 

Pi,  l 3 2i  + 1, 

<J/,i  2i ; 

6i(1  = 4i  + 1 

Pi, 2 - ?J  + 3, 

7/, 2 = 2/  - 1 ; 

«/,2  =41  + 2. 

(3-24) 

Table  3-3  lints  the  optimum  p,  q and  associated  end-phase  and  midphase  array  ratios 
vs  overall  array  length  6(2  < 6 *S  26). 

It  can  be  readily  deduced  from  the  argument  of  the  probability  of  ambiguity  func- 
tion for  independent  arrays, 

jr 

11  ■■  11  ~"~~T  * 

°CH-PRVmZ  + n2 


that  the  optimum  p,  q vs  urray  length  are  just 


m = p = 6 
n ~ g - 1 


independent . 


Table  3-4  lists  all  possible  p,  q sets  (including  the  optimum  sot  defined  above)  for 
three-element  arrays  of  length  6 = 2 to  26.  It  is  interesting  to  note  that  the  ratio  entries 
for  p fixed  have  their  counterparts  in  a table  of  Farey  sequences  (arrays  of  rational  frac- 
tions between  0 and  1)  in  number  theory  os  shown  by  Niven  and  Zuckerman  (20]. 


3.4  Tabulated  Probability  of  Ambiguity  for  Three-  and  Four-Element 
Two-Integer  Set  Arrays 

Based  on  Eq,  (3-21)  from  Sec.  3.2,  pa  for  three-element  arrays  are  listed  in  Table  3-6 
for  Och-PR  = 6°,  10°,  12°,  16°  and  20°.  Table  3-6  provides  the  same  information  for 
four-clement  (independent)  arrays. 
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Table  3-5— Probability  of  Ambiguity  vm  Channel-pair  Phase  Error  for 
Threo-element  Interferometers 


Array 

Length, 

K 


Probability  of  Ambiguity  for  Given  och.pr 


6:1  , 6:5 


10°  I y 12°  o « 15' 


IB 


■ 

wUkkI'- 

8:3 

, 8:5 

8:1 

, 8:7 

9:4 

, 9:5 

9:2 

, 9:7 

9:1 

, 9:8 

10:3 

, 10:7 

10:1 

, 10:9 

11:5 

, 11:6 

11:4 

, 11:7 

11:3 

, 11:8 

11:2 

, 11:9 

11:1 

, 11:10 

12:5 

, 12:7 

12:1 

, 12:11 

<lE-8 

<lE-8 

<1E*7 

<3E-7 

<2E-6 

4E-6 

1.1E-5 

2.6E-5 

5.1E-5 


1.61E-4 

1.89E-4 

2.57E-4 

3.89E-4 

6.33E-4 


- 

<lE-7 

<6E-6 

<6E-7 

3.18E-5 

8.74E-4 

3.63E-6 

5.79E-4 

6.91E-3 

8.67E-5 

1.06E-3 

t 

S.83E-3 

1.23E-3 

7.06E-3 

0.0311 

3.09E-3 

0.0137 

0.0486 

3.95E-3 

0.0163 

0.0547 

6.05E-3 

0.0222 

0.0673 

0.0101 

0.0321 

0.0171 

0.0469 

0.112 

0.0212 

IjXj » * 9 

0.124 

0.0279 

0.143 

0.0351 

Rff  r*m 

0.160 

0.0429 

0.0916 

0.177 

0.0592 

0.116 

0.208 

0.0592 

0.116 

0.208 

0.0620 

0.120 

0.213 

0.0676 

0.128 

0.223 

0.0761 

0.139 

0.0876 

0.154 

0.256 

0.0847 

0.151 

0.250 

0.119 

0.193 

0.298 

*•*  optimum  ratio  for  givttn  V. 
t-  = pa  < IE-10. 


NRL  REPORT  8005 


Table  3-6— Probability  of  Ambiguity  vs  Channel-pair  Phase  Error  for 
Two-integer  Set  Four-element  Interferometers 


Maximum 
Spacing,  2 


Probability  of  Ambiguity  for  Given  oCp 

-PR 

o=6° 

a = 12° 

B9 

o = 20° 

<lE-7 

2E-6 

1.48E-4 

4.43E-3 

- 

1.3E-5 

2.76E-4 

3.61E-3 

0.0290 

- 

4.16E-4 

3.26E-3 

0.0186 

0.0776 

C1E-8 

3.08E-3 

0,0137 

0.0485 

0.139 

C1E-6 

0.0'l09 

0.0339 

0.0897 

0.203 

8E-6 

0.0266 

0.0628 

0.137 

0.264 

7.0E-6 

0.0468 

0.0976 

0.186 

0.320 

3.41E-4 

0.0733 

0.136 

0.232 

1.12E-3 

0.103 

0.174 

0.277 

0.416 

2.79E-3 

0.136 

0.213 

0.319 

0.455 

=■  pu  < IE-10, 


Aa  mentioned  previously  in  Sec.  3.2,  »CH-PR  *n  the  range  10°  to  12°  is  typical  of 
present  technology  for  wideband  microwave  systems— -2  to  3 octaves  of  frequency  coverage- 
with  lower  and  upper  absolute  limits  of  600  MHz  and  18  GHz.  aCH.PR  on  the  order  of 
6°  typifies  the  channel-pair  errors  in  the  same  category  of  systems  utilizing  RF  calibra- 
tion. The  probability  of  ambiguity  for  aCiI.pR  = 16°  and  20°  has  boon  included  as  an 
admittedly  crude  estimate  of  the  performance  that  might  be  experienced  with  an  inter- 
ferometer operating  in  a severe  multipath  situation, 

A pa  criterion  of  0.01  to  1 percent  maximum  ambiguity  is  often  set  by  designers  of 
systems  to  be  used  for  location  by  triongulution  with  multiple  DF  cuts.  As  Fig.  3-2 
shows,  with  oCH_pK  in  the  range  10°  to  12°,  the  maximum  allowable  array  length  is  be- 
tween 8 to  6 half-wavelengths. 

The  angular  accuracies  of  systems  with  baselines  this  short  may  be  unacceptable. 
This  is  the  basic  reason  for  synthesizing  and  employing  urrays  with  additional  elements. 

In  Sec.  5 It  will  be  shown  that  the  optimum  four-element  array  of  length  2 = 26  half- 
wavelengths has  pa  = 0.116  percent  for  Och-PR  ~ 12°,  mid  pa  - 1.17  percent  for 
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4.0  FUNDAMENTALS  OF  FOUIl-ELFMENT  THREE-INTEGER 

SET  INTERFEROMETERS 

The  theory  of  three-element  interferometers  in  fairly  well  known  (disregarding  some 
minor  inconsistencies  in  the  work  referenced  in  Sec.  3.0).  There  does  not  appear  to  be 
any  corresponding  body  of  theory  available  for  four-element  arrays,  especially  with  re- 
gard to  the  question  concerning  the  existence  of  an  optimum  four-element  array  con- 
figuration. With'an  understanding  of  the  material  presented  in  this  section,  the  ESM 
designer  will  be  prepared  for  the  exposition  of  the  array  synthesis  techniques  of  Soc.  5.0, 
and  he  will  be  able  to  apply  them  without  difficulty. 

A four-element  interferometer  can  be  characterized  by  4 antenna/receiver  channels, 
4-1  = 3 available  channel-pair  phasa  differences  (necessary  and  sufficient  for  ambiguity 
resolution  of  the  overall  array),  and  4 -2-2  subarrayB  formed  by  channels  considered 
three  at  a time.  Section  4.0  introduces  several  four-element  array  configurations  (some 
of  which  are  not  optimum)  which  have  been  widely  used  in  the  past.  These  configura- 
tions are  compared  mainly  on  the  basis  of  how  efficiently  the  signal  processing  and  math- 
ematical operations  needed  in  ambiguity  resolution  can  be  mechanized. 

A discussion  of  resolvable  and  unresolvable  ambiguities  in  the  two  three-element  sub- 
arrays that  constitute  a four-element  array  emphasizes  that  there  is  no  need  for  either  of 
the  two  subarrays  to  be  nonrodundant.  That  is,  there  is  no  requirement  that  either  of 
the  subarrays  be  capable  of  unambiguous  operation  as  a distinct  three-element  array —as 
long  as  the  three  spacing  integers  defining  the  overall  array  are  relatively  prime. 

Each  subarray  is  characterized  by  a subarray  ratio  One  member  of  the 

subarray  rutlo  (either  or  ) contains  a factor  common  to  one  member  of  sub- 
array ratio  Ut2  (either  m2  or  n2),  dependent  on  the  overall  array  configuration.  This 
highest  common  factor  a is  introduced,  and  the  manner  in  which  the  ambiguity  tolerance 
of  suburray  2 is  increased  because  of  this  factor  a is  fully  explored. 

The  section  closes  with  a development  on  canonical  array  configurations  in  four- 
element  arrays.  It  is  shown  that  the  cascaded  end-phase  configuration  of  four-element 
interferometers  is  optimum  from  considerations  of  efficiency  of  hardware  usage  and 
probability  of  ambiguity.  This  is  believed  to  be  a new  result,  and  it  was  obtained  with- 
out recourse  to  the  analytical  artifice  of  a fictitious  off-axis  fifth  channel  used  by 
Hanson  [211. 


4.1  Three-Integer  Set  Interferometers 

Interferometers  whose  performance  could  be  specified  as  a function  of  the  two 
spacing  integers  p and  q were  classified  in  Sec.  2.2  according  to  which  channel  was  des- 
ignated the  phase  reference.  For  arrays  formed  with  the  number  of  channels  n - 4 and 
n = 3,  classification  can  be  done  on  the  basis  of  the  channel  (or  channels)  used  as  phase 
reference.  The  number  of  possible  array  configurations  increases  rapidly  with  n.  It  will 
be  shown  later  (in  Sec.  4.4)  that  in  four-element  interferometers,  the  probability  of 
ambiguity  is  dependent  on  the  array  configuration  as  well  as  on  array  spacing  integers 
p,  q,  and  r—in  contrast  to  three-element  interferometers. 
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Figure  4-1  illustrates  several  configurations  possible  in  four-element  arrays.  Param- 
eters for  these  arrays  of  interest  to  a system  designer  are  given  in  Table  4-1, 

One  array  comparison  parameter  is  the  number  of  phase  comparisons  utilized  in 
resolving  the  overall  array  ambiguities.  The  number  of  these  ambiguities  is8  = p + q + r, 
for  0a  within  -90°  to  +90°;  and  p,  q,  and  r represent  the  element  spacings  in  half-wave- 
lengths  at  the  operating  frequency.  The  minimum  possible  number  of  theso  comparisons 
is  obviously  (n  - 1),  since  the  use  of  less  than  (n  - 1)  phase  comparisons  implies  that  the 
phase  information  from  one  or  more  channels  has  been  discarded.  The  cascaded  mid- 
phase array  requires  four  comparisons  in  resolving  p + q + r = fi  ambiguities.  All  of  the 
other  array  configurations  require  (n  - 1)  = 3 comparisons. 


CASCADED  END-PHASE 


*t,5  # 


IBI  CASCADED  MIDPHASE 

a#'  Hi  ao  . HI 

q &lln\  r 


HYBRID  MIDPHASE 
P « q + r 

T ' q ' 

♦ u ^2.4 


1 K7 


ADDITIVE  MIDPHASE 


y ' y ’ y 


’1.2 


MULTIPLE  END-PHASE 
p+q  __  p+q+r 


•Thla  ratio  uaad  altar  p,  q inaolvod  by  (9t , 


Fig,  4-1 -- Various  configuration*  for  four-olemcnt  nrrayH 


= P/9 
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Note:  For  the  ‘ infinite-resolution”  phase-measurement  systems  considered  here,  imple- 
mentation of  other  comparisons  in  addition  to  some  minimum  set  is  easily  accom- 
plished. Thus,  the  phase  difference  <l>1)3  is  equal  (formally)  to  + ^2,3>  in' 
eluding  any  errors.  In  quantized-phase  systems,  the  direct  representation  of  <t>13 
may  differ  from  <1>1  3 obtained  by  summing  quantized  representations  of  d>12  and 
d>2  3.  Hence,  depending  on  the  degree  of  quantizing  implemented,  the  system  de- 
signer might  be  advised  to  obtain  <b1>3  directly  by  use  of  the  additional  phase 
comparator. 

A second  array -configuration  comparison  parameter  is  the  number  of  phase  refer- 
ences needed  to  Implement  a particular  array.  The  number  of  high-level  input  signals 
that  must  be  provided  is  equal  to  the  number  of  phase  references.  Minimizing  this 
number  is  desirable,  since  the  analog  phase  comparators  in  widespread  use  in  present 
wideband  receiving  systems  usually  require  a drive  level  on  one  of  the  input  ports  to  be 
at  least  10  dB  stronger  than  the  other  to  facilitate  accurate  recovery  of  sin  \p  and  cos  cp, 
the  quadrature  components  defining  phase  difference. 

Another  array-configuration  comparison  parameter  is  whether  or  not  the  8 = p + q + r 
ambiguities  associated  with  the  overall  array  length  are  resolved  directly  us  a consequence 
of  the  ambiguity-resolution  process.  If  not,  additional  operations  and  calculations  will  be 
needed.  Table  4-1  shows,  for  n - 4,  that  only  the  cascaded  end-phase  and  multiple  end- 
phase  array  configurations  meet  this  criterion. 

A final  array-configuration  comparison  parameter  is  the  subarray  spacing  length  com- 
mon to  each  of  the  two  ratios  characterizing  the  arrays  shown  in  Fig.  4-1.  The  subarray 
spacing  length  common  to  the  two  ratios  influences  the  order  of  the  ratio.  An  example 
is  a comparison  between  the  midphase  configuration  and  the  cascaded  end-phase  con- 
figuration for  p ~ 3,  q =»  1,  and  r = 7 half-wavelengths. 


Midphase 


CD  J?  „ _1 

2m  “ r ~ 7 


Common  spacing:  q = .1  , 


Previous  analyses  into  the  theory  of  multiple  element  ambiguity  resolution  have  proceeded 
on  the  assumption  that  composite  arrays  having  low  probability  of  ambiguity  could  only 
be  achieved  by  combining  subarrays  of  low  order.  Thus,  inordinate  inteicst  has  been 
placed  on  the  hybrid  midphase  and  udditive  midphaso  configurations.  It  will  be  shown 
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later  in  this  section  that  the  cascaded  end-phase  and  cascaded  midphase  configurations 
exhibit  tho  lowest  probability  of  ambiguity  for  given  p,  q,  and  r because  they  have  the 
largest  common  spacing  p + q,  as  the  arrays  are  configured  in  Fig,  4-1,  It  will  aho  be 
shown  that  if  the  subarrays  are  redefined  properly,  all  array  configurations  in  Fig.  4-1 
axe  equivalent  except  the  multiple  ond-phase  configuration.  This  is  a new  result  in  the 
theory  of  four-element  Interferometers. 

In  summary,  the  four-element  cascaded  end-phase  configuration 

• Employs  the  minimum  possible  number  of  phase  comparisons— 3 

• Requires  the  smallest  number  of  phase  references— 1 

• Resolves  the  overall  length  p + q + r = 8 directly  as  a consequence  of  the 
ambiguity  resolution  process 

• Exhibits  the  largest  spacing  common  to  the  two  subarrays  constituting  tho 
overall  urray— p + q. 


4.2  Resolvable  and  Unrosolvable  Ambiguities  in  Subarrays  of  Three  Elements 

In  this  section,  tho  distinction  between  resolvable  and  unresolvable  ambiguities  In 
subarrays  will  be  defined.  The  material  will  be  useful  in  understanding  the  discussion  on 
ambiguity  constraints  in  the  section  following, 

For  a four-element  cascaded  end-phase  array,  os  shown  in  Table  4-1,  the  two  sub- 
arrays are  defined  by 


/ 


and 


to  p + g ...  mi 
t*1  p ~ nl 


(41a) 


_ P + Q + r „ ™ 2 

~ p + q "no 


(4-lb) 


Analogous  to  the  requirements  for  an  end-phase  three-element  interferometer,  each  of  tho 
numbers  and  nt  (i  ~ 1,  2)  associated  with  the  two  subarrays  must  meet  the  criterion 
for  a realizable  interferometer.  Thut  is,  ittj  und  must  be  relatively  prime  integers. 

There  is  u corresponding  restriction  that  the  spacing  integers  p,  p + <7,  and  p + </  + r be 
relatively  prime,  ns  will  be  shown  in  tho  next  section  on  ambiguity  constraints. 

As  an  example,  consider  tho  cascaded  end-phase  array  defined  by  p = 3,  q = 3,  and 
r ~ 10  half-wuvelengths  at  tho  operating  frequency.  The  subarray  parameters  ure 


Subarray  l 


P + Q 3 6 ,»  m\  2 

p 3 3 ’ u’  nl  1 ’ 
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Subarray  2 

p + q + r 3 + 3 + 10  16.  q m2  _ 8 

p + q 3 + 3 “ T * 2*  “ T* 

The  individual  three-element  interferometers  are  respectively  “three-times  ambiguous" 

[(p  + q)/m1  = p/nl  =--  3],  and  “two-times  ambiguous”  f(p  + q + r)/m2  = (p  + q)ln2  “ 21. 
The  overall  array,  as  well  as  the  individual  subarrays,  can  be  resolved,  however,  because 
the  four-element  array  was  synthesized  from  two  realizable  ratios.  An  equivalent  way  of 
stating  this  is  to  say  that  a four-element  array  is  resolvable  provided  at  least  two  of  the 
three  spacing  integers  are  relatively  prime.  Thus,  the  array  formed  by  p = 3,  q a 3,  r ■ 9 
is  not  resolvable,  whereas  the  array  formed  by  removing  the  common  factor  3,  p'  * 1, 
q » 1,  r'  = 3,  resolvable. 

In  the  above  example  involving  subarray  1 with  p ■ q ■ 3 half-wavelengths,  theie  are 
p + q m 6 angles  of  arrival  over  a -90°  to  +90°  field  of  view  that  will  give  rise  to  the 
sumo  phase  code  V>i,3.,,vi»llable  = ^1.3  (mod  2ff)-  This  is>  of  cour8e.  » restatement  of  the 
principle  that  a three-element  array  of  length  8 ■ p + q (in  half-wavelengths)  cannot  ex- 
ceed mx  A/2  wavelengths;  otherwise  there  will  exist  ambiguities  beyond  the  capability  of 
the  channel-pair  spaced  p half-wavelengths  to  resolve. 

Subarrays  which  have  no  unresolvable  ambiguities  over  the  field  of  view  are  called 
“unambiguous.”  Subarrays  which  exhibit  one  or  more  unresolvable  ambiguities  over  the 
field  of  view  are  therefore  called  “ambiguous.” 

Table  4-2  lists  all  possible  four-element  arrays  of  length  i<=p*q+r<=  16  half- 
wavelengths,  with  the  associated  subarray  ratios  and  iR2#  for  the  cascaded  end-phase 
configuration  (subject  to  an  example  constraint  of  p =*  3).  The  purpose  of  the  listing  is 
to  illustrate  the  conditions  unambiguous  and  ambiguous. 

The  number  of  unreBolvabTe  ambiguities  for  each  subarray  are 

subarray  1 A j = (p  + q)  - (4-2a) 

subarray  2 A2  =>  (p  + q + r)  - m2  . (4-2b) 

Of  the  entries  in  Table  4-2,  only  4 out  of  12  have  both  subarrays  1 and  2 unambig- 
uous. It  is  stressod,  however,  that  the  usual  impetus  for  implementing  arrays  of  four  or 
more  elements  is  to  achieve  the  higher  angulur  resolution  implied  by  larger  overall  spacings 
(with  an  acceptably  low  probability  of  ambiguity).  Thus,  the  array  designer  is  normally 
indifferent  to  employing  angular  estimates  from  any  spacings  except  the  overall  array 
spacing  8.  An  exception  to  this  statement  occurs  if  all  of  the  array  spacings  are  used  to 
form  an  estimate  of  the  angle  of  arrival,  e.g.,  as  in  perhaps  maximum-likelihood  processing 
of  the  electrical  phases  from  the  appropriate  channel-pairs. 

4.3  Ambiguity  Constraints  and  Probability  of  Ambiguity  for  the  Four-element 

Cascaded  End-phase  Array 

Explicit  relations  for  the  ambiguity  constraints  in  the  cascaded  end-phase  configura- 
tion of  a four-element  array  will  be  derived  in  this  section.  The  derivations  for  the  other 
array  configurations  given  in  Sec.  4,1  are  similar. 
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Table  4-2— Parameters  of  Four-element  Arrays  of  Length  2 = p + y + r=16 
(subject  to  the  constraint  p = 3) 


Spacing  Integers 
p - q - r 


*1  = (P  + q)‘P  *2  ■ (P  + q + r):(p  + <?) 


3 - 

2 - 

11 

3 - 

4 - 

9 

3 - 

8 - 

6 

3 - 

10  - 

3 

3 - 

1 - 

12 

3 - 

5 - 

8 

3 - 

7 - 

6 

3 - 

11  - 

2 

3 - 

6 - 

7 

3 - 

12  - 

1 

3 - 

3 - 

10 

3 - 

9 - 

4 

Number  of  Unre- 
solvablo  Ambiguities, 
At  in  31* 

mmmm 

u 

U 

u 

U 

u 

U 

u 

U 

u 

12 

u 

14 

u 

8 

u 

8 

6 

U 

10 

4 

u 

8 

8 

12 

6 
7 

18:11 

16:13 


16:  9 
16:15 


1.  U ■*  unambiguous  (no  unr«solvable  ambiguities). 

2.  >ti,  number  of  unreiolvable  ambiguities  » (p  *q)~  mj,  In  !Ri. 

3.  A a,  number  of  unreiolvable  nmblKultlea  m (p+q*r)  - m2,  In 


It  will  be  convenient  to  introduce  expressions  for  the  array  spacings  in  termB  of  the 
subarray  ratios  and  !R2(,. 

Consider  a four-element  cascaded  end-phase  array  (see  Fig.  4-la).  The  two  subarray 
ratios,  with  the  subscript  e dropped  for  brevity,  are 


K a ^LL  - „ JLU. 

1 = nx  " nj  p 


(4-3u) 


fc  I2.2  ^?2_  _ p + q + r 

2 = "2 ’ " cm'»  " P + 7 


(4-3b) 


with  a an  integer  introduced  for  generality,  to  account  for  possible  common  factors  in 
mj  and  n2  (common  factors  in  and  m2  are  irrelevant).  Another  ratio  Ro»  the  product 
of  1RX  and  !R2,  may  bo  defined  as 
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Identifying  like  terms  in  the  denominators  of  the  two  forms  for  1R0  provides 


, «1«2 
p = ~ 

Substituting  this  value  for  p into  Eq.  (4-3a)  gives 


2 . 

P + Q _ _o 

p “ n,na  “ »i 

a 


and  solving  for  q. 


q = 


£2 

cv 


('»l  ~ «i)  • 


Adding  Eq.  (4-4)  to  Eq.  (4-6)  gives 


P 


+ q - 


5 


Substituting  this  value  for  (p  + q)  into  Eq.  (4-3b)  gives 


p *•  q + r 
~~P  * Q 


»«l«2  , 

a £ . ^2_ 

“ »»1»2  ' n2 


and  solving  for  r, 


mi 

r “ “ (m2~n 2). 

Summarizing,  wo  have 

"l"a  ') 

P = — 

n2 

q « --  (mt  - nO  > ii  = p + q + r 
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mi  p__+_q . 
nt  p ’ 


A m2  P + q + r 
^2  = 7^  = -p-*~<T' 


a = factor  common  to  mi  and  n2- 

With  the  aid  of  subarray  ambiguity  diagrams,  the  ambiguity  constraints  for  a four- 
element  cascaded  end-phase  array  will  now  be  derived. 

Figure  4-2  shows  the  subarray  ambiguity  diagrams  for  the  array  ratios  >=  mllnl  a 4:3; 
$2  = /n2./«2  “ H;4.  Since  = n2  “ 4,  the  common  factor  a in  Eq.  (4-8)  is  4,  and  p,  q, 
and  r are  given  by 


n ..  8 X_4  _ o 
P “ A 3 


q rn  ± (4  - 3)  - 1 


r - j (11  - 4)  - 7 


y il  ~ p + q + r 
4 X 11 


11  . 


The  umbiguity  diagram  for  subnrray  2,  unresolved  by  subarruy  1,  has  umbiguity  Inter- 
cepts of 


A^l,4- 


AI 


7T 

n2 


<4B“) , 


AlPl,3-Al 


7 r 

m2 


- (16.36°) . 


Similarly,  the  ambiguity  diugrum  for  subarruy  1 has  ambiguity  intercepts  of 


<Vi  ,:i-ai 


JL 

«i 


It  will  bo  noted  that  the  trajectories  for  subarruy  1 aro  lubelled  with  unprimed  numerals 
below  the  trajectories  of  slope  4:3,  whereas  the  primed  numerals  above  the  trajectories 
show  parumetrlcally  the  equivalent  trajectory  for  subarray  2.  This  dual  labelling  reflects 
the  fact  that  us  l)a  vurles  from  - 90°  to  +90°,  the  channel  pair  spucod  (p  + q)  half- 
wavelongths  manifests  4X  2ir  rad  phase  change,  whereas  the  channel-pair  spaced  (p  + q + r) 
half-wavelengths  manifests  11  X 2rr  rad  phase  chunge. 

The  spacing  integers  p - 3,  q - t,  and  r-  7 for  this  example  have  been  chosen  so 
that  both  subarruy  1 und  suburruy  2 are  unambiguous  in  the  sense  of  the  development  of 
Sec.  4.2.  That  is,  the  appropriate  processing  of  ouch  of  the  suburruy  phuso-compurator 
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outputs  can  provide  a single  (though  inexact)  estimate  of  0a  for  sufficiently  small  phase 
errors, 

The  channel-pair  errors  for  the  two  subarrays  are 
Subarray  1 


Subamty  2 


Ayfy  = A^  - Asp2  = A^i  a 

Atp^j  - A<0i  - A 1^3  = Ai01i3  . 

Ai0Sa  « A<0t  - Ai03  = A«0m 

Ai, pLi  * Ai0t  - A 104  «=  Ai01j4  , 


(4»9a) 


(4-9b) 


where  Si%  L(  •*  small  und  large  spacings  for  subarray  i,  For  zero-mean  channel  phase 
errors,  the  aisociuted  standard  deviations  of  channel-pair  phase  error  are 


Subarroy  1 


% ■ fo?  + ogl1/a  « o1(2 

% “ t")  + «3l1/2  = °1,3  • 


(4-10a) 


and 


Subarmy  2 


% * t°l  + = °1,3 

°LU  - f°l  + 0a]1/2  = °M  • 


(4-10b) 


For  all  0\j  « «oh  PR>  j ~ 2 to  4,  the  ambiguity  vuriublo  in  Kq.  (3-18)  for  ouch  of  the 
two  suburraya  is  given  by 


% 3 "CH-PRt”1?  ~ 2 PflBI/Bj  + nf] 


1/2 


r3  ®ch-prI mf  - rn in,  + nf)m  , 
since  pe  = +0.5  tor  end-phase  arrays,  Consequently, 


(4-11) 
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a)f> j - (7ch-pr  \A2  “ 4 * 3 + 3Z 
= WCH-PR 


for  subarray  1 , 


a<t>z  ~ °CH-PR  \/ll2  " 11  • 4 + 4- 
= ^CH-PR  \/93 


for  subarruy  2 , 


Then, 


Suppose  that  the  rms  electrical  phase  error  between  channel  pairs  (Jch-pr  = 12.5°. 


o01  = 12.6°  v/13  - 45.07°, 
0^2  = 12.5°  s/93  « 120,56°  , 


The  integration  limits  in  Eq.  (3-19)  uro 

1 firt0 

Subarray  1:  ±~—  « ±3.994  , 

i 

1 RO° 

Subun'ay  2:  ±—--±1.493. 

°01l 

When  Eq,  (3-21)  1h  used,  the  probability  of  ambiguity  of  the  individual  subarrays  is 


Suburray  1:  p.  - orfo  ( — 


/ 180°  ' 

W^iu 


- 6.497  X 10“  r> , 


Subarruy  2:  pa  - erfc  f •— 4^--— | = 0,1364 

\V  2 OpJ 


The  pa  for  subarruy  2 is  over  three  orders  of  magnitude  larger  than  for  subarruy  1,  This 
dramatically  illustrates  the  desirability  of  utilizing  tho  information  from  chunnols  1-2-3 
(subarray  1)  to  improve  tho  ambiguity  porformanco  of  suburray  2,  composed  of  chunnols 
1-3-4. 

The  ambiguity  diugruins  for  the  two  suburruys  illustrate  how  this  muy  bo  accomplished. 
Suppose  thut  0a  - 0°  and  thut  A*>i,2  * 36.818°,  « 49.091°,  and  A*>M  ---  46°,  (Note: 

These  error  sots  huve  boon  chosen  to  place  tho  error  sots  ( AyJt  y,  A^j  3)  os  well  as  (A^j  a , 
Ay?  1,4)  directly  on  ambigulty-plune  trajectories,  but  this  assumption  is  not  essential  to  the 
following  development.)  Now,  tho  error  set  ( Ay?t  2»  ^t,a)  lies  011  trajectory  0 for  sub- 
array 1 (lower  diugrum  In  Fig.  4-2),  whereus  the  error  set  (A^j  3,  A^1  4)  lies  on  trujectory 
3'  for  suburray  2 (upper  diagram  in  Fig.  4-2). 


-«** « 1 1 ■ * . »*< 
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If  the  phase-comparator  output  from  only  the  larger  baseline  in  each  of  the  two  sub- 
arrays was  used  separately  to  form  an  estimate  of  the  angle  of  arrival,  the  results  would 
bo  as  follows 

Subarray  1: 

0ai  = sin'1 


Subarray  2: 

60s  ■ sin"1 


3 X 360°  + 45° 


11 

. 2 


= 34.62° 


X 360° 


0 X 360°  + 49.091° 


4 


= 3.91° 


X 360° 


Thus,  although  the  channel-pair  phase  errors  are  comparable  in  magnitude,  Bincc  trajectory 
3'  for  suburray  2 is  not  the  truo  trajectory,  using  the  phase  information  from  subarray  2 
alone  provides  a grossly  incorrect  angle-of-arrival  estimate. 

Given  that  the  error  sot  (Aipj  2,  for  subarray  1 lies  within  the  boundaries 

around  the  trajectory  associated  with  the ’true  angle  of  arrivul,  the  following  items  may 
be  noted 


• The  true  trajectory  for  suburruy  1 is  0, 

• The  true  trajectory  for  suburruy  2 is  O', 

• The  contiguous  subarray  2 trajectories  (plotted  on  the  trajectories  for  suburruy  1) 
uro  l'  for  0(l  (+)  and  10'  for  0U  (-), 

• The  adjacent  ambiguous  trajectories  for  subarray  2 in  Fig.  4-2a  are  3'  and  8', 

which  are  n2  -1  = 3 greater  or  leas  than  (mod  11)  0':(0  + 3-  3;  11  - 3 8), 

• Trajectories  6'  (3  + n2  - 1 ~ 6)  und  5'  (8  • no  + 1 ® B)  lie  midway  between  tra- 
jectories 0'  and  1',  und  ()'  and  10',  respectively. 


The  pa  calculations  made  above  show  thut  it  is  almost  a certainty  that  the  truo  tru- 
joctory  for  suburruy  1 is  0 for  the  given  example,  subject  to  channel-pair  errors  of  12. 5° 
rms,  with  Gaussian  distribution.  The  only  suburruy  2 trajectories  possible,  given  trajec- 
tory 0 on  subarray  1,  arc  O',  1',  and  2u'.  Trajectory  3'  In  Fig.  4-2a  for  suburray  2 is 
udjucont  to  trajectory  O',  whereas  trajectories  1'  and  2'  are  four  trajectories  to  the  right 
and  left,  respectively,  of  trajectory  0.  Hence,  the  true  trajectory  for  suburruy  2 is  0'  * 
u very  high  probability. 


to 


Subarray  1,  in  effect,  increases  the  width  of  the  ambiguity  boundaries  around  the 
phase  trajectories  for  suburruy  2 by  the  fuctor  - 4.  This  is  shown  in  Fig.  4-3,  the 
umbiguity  diagram  for  !R2  m2:n2  11:4,  resolved  by  m-\  In  in ( : n i - 4:3.  Rather 

than  u square  2rr  on  each  side,  the  ambiguity-plane  surfuco  Is  now  2tr  long  on  the  m.,  axis, 


► 
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RESOLVED 

BY: 


THAJECTORv 0 


TRAJECTORY  1 TRAJECTORY  2 TRAJECTORY  3 


4*t.«  AI  " ml  ~ ^VI.S.  AI  u m,  15J- 

(RESOLVED  BY  m, I ■»  4 ~-  " »r  1180°!  (RESOLVED  BY  m,)  -•.<-“-(06.46°) 

Fig.  4-3— Ambiguity  diagram  for  !R2  “ m2:n2  " Hi4  resolved  by  mi  in  D\i  ■ mitni  ■ 4:3 


but  is  increased  to  4 X 2tt  long  on  the  n2  axis,  It  is  readily  seen  that  as  trajectory  2' 
pierces  the  first  27r  ordinate  to  the  right  of  Aipj  3 = 0,  there  is  no  need  to  translate  this 
intersection  point  back  to  the  ordinate  at  Ag>i  3’  = 0.  The  new  ambiguity  boundary  be- 
tween trajectory  0*  and  1'  (in  reality  trajectory  6'  in  the  previous  illustration)  is  two  times 
the  distance  between  trajectories  0'  and  3'  in  Fig.  4-2a. 

In  this  modified  ambiguity  diagram  for  subarray  2,  resolved  by  subarray  1,  the  am- 
biguity intercepts  are 

M.A -A!  “ n»t  » 4 ~ = it  (180°)  , 

^l,3-AI  “ "l  ~ « 4 (65.45°). 

The  new  integration  limits  in  ISq.  (3-19)  for  subarray  2 are 

m1  180°  , 4 x 180° 


Thus,  the  pa  for  subarruy  2,  resolved  by  subarray  1,  is 


fmx  180°\ 

“ x 10  • 


If  independence  between  ambiguities  in  subarray  2 and  subarray  1 is  assumed,  the 
overall  array  pu  is  given  by 
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Pa-overall  ~ ^ ~ d " Pad  d ~ Pad  • 


(4-12) 


For  the  example  above, 


Pa-overall  = 1 ~ d “ 6-497  X 1<T5)(1  - 2.3  X 10-9) 

- 6.49723  X 10- 5 . 

This  is  a negligible  increase  from  the  pa  for  subaxray  1 alone.  However,  the  overall  array 
has  a length  2 = 11  helf-wavelengths,  whereas  subarray  1 has  a length  p + q » 4 half- 
wavelengths,  Thus,  the  accuracy  of  angle  of  arrival  estimates  using  the  whole  array  is 
almost  three  times  (11/4)  better  than  for  subanuy  1 alone. 

The  geometric  arguments  above  illustrate  how  the  ambiguity  boundaries  for  sub- 
array 2 in  a four-element  cascaded  end-phase  interferometer  are  extended  by  a factor  re- 
lated to  the  mx  spacing  integer  in  subarray  1.  This  ambiguity  constraint  relaxation  on 
suburray  2 can  also  be  developed  analytically,  as  shown  below. 

The  ambiguity  constraints  for  the  two  subarrays  (considered  as  isolated  three-eloment 
arrays)  cun  be  written  in  several  forms 


Array  Ratio  Form 


Subairay  1 IniAy^a  - 12|  ^ n 

Subarray  2 l«2^l,4  " £ n , 

Spaaing  Integer  Form 

Subarray  1 |pAipia  - (p  + q) Ay?12|  ^ ir 
Subarray  2 |( p + q)Aip1>4  - (p  + q + r) Api^l  ^ n , 
Spacing  Integer-Channel  Error  Form 

Subarray  1 l-gApj  + (p  + q)  Ap2  - pAip3|  ir 
Subarray  2 1-rAySj  + (p  + <7  + r)A<p3  --  (p  + </) Aip4|  rt 


where 


9 “ -5-  («h  -»i) 

mj 

r = --  (ms-n2) 


a = factor  common  to  m,  and  n2. 


p + q •= 


m 1 m<> 

P + q + r = — = 2 


(4-13a) 

(4-13b) 


(4-1 4a) 
(4-14b) 


(4-16a) 

(4-lBb) 


I 


v 
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If  tho  subarray  ratio-in teger  forms  for  the  spacing  integers  p,  q,  r,  ...  are  used,  and  if  sub- 
array  2 is  resolved  by  subarray  1,  the  ambiguity  constraint  Eq.  (4-1 5b)  becomes 


ftl  i rnifna  wtrto  I 

--%•  (mz-n2)A^i  + “ — Ap3 Aip4|  t ^ • 


m1m2 


™Jt«2 


(4-16) 


Multiplying  both  sides  of  the  inequality  above  by  a/mx  yields 

l~(m2  - n2)Awi  + m2Aip3  - n2Ap4|  ^ off, 

which  simplifies  to 


|n2A^i(4  --  m2A<pli3|  ^ ait 


Ambiguity  constraint  for  subarray  2 (4.17) 

resolved  by  subarray  1 ' 1 


using 


Ap]  --  AiC4  k Aip4>4,  Aipt  - Aip3  ---  Alp!  3 . 


The  ambiguity  constraint  discussion  in  this  section  culminating  in  Eq.  (4-17)  can  be  sum- 
marized for  the  four-element  cascaded  end-phase  interferometer  as  follows 


Four-Element  Cascaded  End-phase  Interferometer 
• Array  Constants 

ml  P + 9 


Ratios : 'H  l 


«1  ~ p 
m2  p + q + r 


”2 


P + <1 


(4-18a) 

(4-18b) 


Spacing 

Integers: 


«1"2 


mi 


p = — - , q ~-  — - (mt  - «! ) , r ~ — (m2  - n2) 
m1m2 


p + q + r « ti- 


er 


(4-18c) 


a -•  factor  common  to  mx  and  n2 
• Ambiguity  Constraints  (each  siubarruy  considered  separately) 

Subarray  1:  InjAipj^  - m t Av?  1 1 a ! <c  it  (4-18d) 

Subarray  2:  |n2AiPi4  --  m2AiPi|3l  <:  tt  (4-18e) 
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I 


i 

i 


• Ambiguity  Variables 

°<t>t  = Och-prI"1?  - 2p1,,m1rt1  + «fl1/2 
°0a  = °oh  PR[ml  - 2pa2m2rt2  + nj]112 


• Integration  Limits 
180° 


Subarray  1:  ± 


Subanay  2:  ± 


«180° 

% 


* Probability  of  Ambiguity 


. / 180°  \ 

“ erfc(-7^H 

pa  = erfc 


W%'0<Pt/ 

«180° 


(4-18f) 

(4-18g) 

(4-18h) 
(4*1 8i) 

(4-1 8j) 
(4-18k) 


The*  relutlons  above  were  derived  analytically  and  hold  for  all  four-element  cascaded 
end-phuse  arrays,  regardless  of  whether  the  individual  subarrays  are  unambiguous  or  am- 
biguous in  the  sense  of  Table  4-2  in  Sec.  4-2. 

Table  4-3  shows  examples  of  cascaded  end-phase  arrays  falling  into  two  classes: 

Class  A ? ~ p + q + r is  expressible  directly  as  m2 , 

Class  B V can  only  be  expressed  us  m j m 2 /a . 

Note:  The  meaning  of  the  terms  in  puronthesoB  will  bo  explained  below. 

It  should  be  apparent  from  Fig.  4-3  that  a is  the  factor  by  which  the  width  of  the 
ambiguity  diagram,  for  subarrny  2 is  increased.  Even  though  for  the  Oluss  B,  !Rj  4:1, 

!R2  “•  3:1  example,  a Is  only  unity,  it  should  not  be  ussumod  that  this  is  un  urray  with 
poor  tolcruncc  to  amblguitios. 

The  ambiguity  diagram  for  this  array  (Fig.  4-4)  shows  why  this  is  so.  Since  (n  + q + r)l 
(p  + q)  3,  the  trajectories  3',  4',  and  5'  overlay  respectively  the  trajectories  O',  1 and  2'. 

In  the  terminology  of  Sec.  4.2,  suburruy  2 is  ambiguous-resolvable,  The  factor  a is  unity 
(the  common  factor  between  m t = 4 and  n2  ~ 1 is  1);  consequently,  even  with  tho  infor- 
mation from  subarray  1,  tho  ambiguity  diagram  for  suburray  2 cannot  be  extended,  as  wan 
the  case  in  the  previous  example.  There  is  no  need  to  do  this,  howover,  becuuse  suburruy  2 


l»V 


* 
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Table  4-3-Examples  of  the  Two  Classes  of  B'our-element 
Cuscadcd  End-phase  Arrays 


Amv 

Subanay  Ratios 

a 

Array  Spaci 

ng  Integ 

ers 

Class 

(Rj  = mjiMx 

!R2  = m2:«2 

P 

Q 

r 

ft 

A 

4:3 

11:4 

4 

3 

1 

7 

11 

A 

3:2 

11:6 

3 

4 

2 

5 

11 

7:2 

2 

3 

1 

10 

14 

B 

4:3 

(14:4) 

(4) 

3:1 

1 

1 

3 

8 

12 

B 

4:1 

(12:4) 

(4) 

L!_ 

Notei: 

1.  a n factor  common  to  mi  and  na. 


2.P 


mi 


~ (ma  - n 2) ; ' l <=  p +■  q + r 


Is  expressed  in  a form  which  already  is  “reduced"  from  an  ambiguous  higher-order  form. 
This  can  be  illustrated  as  follows. 


Suppose  that  subarray  2 is  expressed  us  - m2:n2  “12:4,  such  ^ '^  ^oation 
equul  H (disregarding  any  common  factors  in  m2  und  n2).  Then  a - 4.  Dir  I P 
of  the  spacing  integer  equations  will  provide  the  same  spacing  integers  p,  q,  and  r as  g 
°f  Z Kw  « aiZopriato  to  die  m2,  n2  set  under  consideration  is  used  This  is  lndlwtod 
in  Table  4-3  by  listing,  whero  necetwary,  a second  set  of  parameters  within  puronthesis 

the  Class  B entries. 


Now,  from  Eq.  (4-18g), 


«*,  - •oii-ntl1**  " 1,-4  * 


= aCn-i»R(4\/T),  or  132.29°,  for  oCH  ,PK  = 12-B°  • 

The  integration  limits  for  subarray  2 (with  m2  - 12,  n2  - 4,  cv  = 4).  using  Eq.  (4-181),  are 


± <*1S0°  . t ±6.443  . 


0<pu  ' 12.5°  X 4 X s/T 
But,  these  are  precisely  the  same  integration  limits  that  result  with  = 3:1,  a - 1,  m 
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0 


<Pi 


woh-pr£32  " 


+ 12]1/2  = 33.07°  , 


± a 180° 
°<f>  a 


± 1 X 180° 
33.07° 


±5.443  . 


for  ®cu.pf{  ~ 12.5  i 


It  is  apparent  that  one  form  of  expressing  $2  f°r  Class  B arrays  exhibits  m2  and  n2 
without  common  factors,  and  with  a*  mi.  But  the  second  form  manifests  m2  and  n2 
with  S.  “ m2,  and  with  a common  factor  0;  = m^.  Thus,  if  m2  and  n2  are  expressed  in  a 
form  that  preserves  any  common  factors,  and  if  the  appropriate  modification  is  made  in 
at,  the  distinctions  between  Class  A and  B disappear. 


In  Sec.  5.5,  the  tabulated  results  of  an  array  synthesis  procedure  developed  earlier 
in  that  section  will  be  presented.  It  will  be  seen  that  consistent  expression  of  m2  and  n2 
directly  os  m2  - p + q + r - S,  und  n2  * p + q,  with  usm]  highlights  the  importance 
of  basing  a four-element  array  synthesis  procedure  on  the  subarray  1 integer  m*.  It  is 
also  stressed  that  the  subarray  2 ambiguity  constraint  development  ending  in  Eq.  (4-17) 
holds  no  matter  in  what  form  m2  and  n2  are  expressed,  provided  a is  correctly  defined. 


4.4  Ambiguity  Constraints  and  pa  for  Four-element  Arrays  of  Various  Configurations 

The  theory  developed  in  Sec.  4.8  for  the  ambiguity  constraints  and  pa  associated 
with  the  cascaded  end-phuse  array  will  now  be  extended  to  other  four-element  array  con- 
figurations. Attention  will  be  concentrated  initially  on  the  other  four  configurations  of 
Fig.  4-1  (see  Sec.  4.1).  Later,  the  development  will  be  generalized  to  all  possible  array 
configurations  using  three  phase  differences— the  minimum  number  required  to  resolve 
the  ambiguities  in  a four-element  array. 

Suburray  ratios  IRj  und  for  the  remaining  four  configurations  of  Fig.  4-1  can  be 
oxpressod  in  terms  of  the  ratio  numerators  und  denominators  by  currying  out  analyses 
similar  to  those  appearing  in  Eq.  (4-3)  through  Eq.  (4-8)  of  Sec.  4.3.  The  results  of  such 
analyses  are  given  in  Table  4»4,  which,  for  convenience  of  reference,  also  lists  the  appro- 
priate parameters  for  cascaded  end-phase  urruys. 

Consider  the  hybrid  midphuse  configuration,  us  shown  in  Fig.  4-lc  os  an  example  of 
how  tho  ambiguity  constraints,  oouplod  with  tho  factor  a— the  greatest  factor  common  to 
the  two  subarray  integers  associated  with  tho  common  spacing  -leud  to  the  pu  function. 
Tho  derivations  of  the  ambiguity  constraints  are  similar  to  those  for  tho  cascaded  end- 
phase  configuration  in  Sec.  4.3,  Eqs,  (4-13)  through  (4-17).  Since  the  derivations  ore 
similar  and  are  based  on  the  sumo  kind  of  geometric  arguments,  it  will  be  convenient  to 
express  the  parameters  for  the  hybrid  midphuse  configuration  in  the  samo  format  as  the 
summary  for  the  cascaded  end-phuse  configuration  following  Eq.  (4-17)  in  Sec.  4.3. 


\ 
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Table  4-4— Common  Spacing,  Factor  a,  and  Spacing  Integers  for  Various  Configurations  of  Four-element  Arrays 


Note : Array  configurations  are  illustrated  in  Fig.  4-1 
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Four-element  Hybrid  Midphase  interferometer 
• Array  Constants 

el  P 

Ratios:  — 


(4-19a) 


c2  q + r 

2e  “ h " <? 


Spacing 

Integers: 


p + q + r=  <!  = 


^2  + he2 


f\(e2  ~ fl) 


(4-19b) 


(4-19c) 


a - factor  common  to  f]  and  f2 

• Ambiguity  Constraints  (each  subarray  considered  separately) 

Subarray  1:  ~ ej Av>2,3 1 ^ (4-10d) 

Subarray  2:  1/2^2, 4 “ «2^2,3l  ^ (4-19e) 

• Ambiguity  Variables 

~ ^CH-PR^l  “ 2pmci  f 1 + /"l  1 1/2  (4*19f) 


°yft  “ aCH-PRf<!2  - ^Pe0,if,2  + 1 


(4-19«) 


• Integration  Limits 


Subarray  1:  £ 


(4-19h) 


a . „ , al80° 

Suburray  2:  ± — 

a'Pi 

• Probability  of  Ambiguity 


(4-19i) 


Pa,  “ erfc 


(4-19J) 


Pa,  = erf{; 


a 180° 


(4-lBk) 
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With  the  information  from  Fig.  4-1  and  Table  4-4,  similar  summaries  can  readily  bo 
generated  for  the  remaining  array  configurations— cascade  midphaso,  additive  midphase, 
and  multiple  midphase. 

Calculated  p are  given  for  all  array  configurations  in  Table  4-5,  with  either  p - 6, 
q - 4,  r - 11  orp  =»  4,  q - 6,  r = 11,  subject  to  oCH ,PR  = 12.5°,  thus  lending  insight  into 
the  importance  of  the  factor  a on  overall  arruy  performance,  The  overall  array  length  ? - 
p + q + r ■=>  21  half-wavelengths. 

In  the  pa  calculations,  it  has  been  assumed  that  ambiguities  in  subarray  2 (resolved 
by  subarray  1)  are  independent  of  ambiguities  in  subarray  1.  [Note:  Later,  in  Sec.  5,  an 
exact  expression  involving  integration  of  a bivariate  Gaussian  density  function  will  be 
given  for  four-element  arrays.  It  will  bo  seen  that  for  most  cases  of  practical  interest,  the 
assumption  of  subarray  independence  causes  negligible  error  in  overall  array  p0.) 

Table  4-5  shows  that  all  array  configurations  are  characterized  by  identical  subarray  1 
p0.  This  results,  of  course,  from  the  fact  that  any  configuration  of  three-element  array, 
ond-phaao  left,  midphase,  or  end-phase  right,  hua  the  sumo  pa  for  given  p,  q,  and  aCH.pR, 

The  fuetor  a for  the  five  configurations  varieB  from  a maximum  of  6 to  a minimum 
of  1.  In  addition,  although  it  is  not  shown  explicitly  in  the  table,  the  subarray  2 ambi- 
guity variable  also  varies.  Consequently,  the  pa  of  subarmy  2 ranges  widely.  The  tabula- 
tion of  overall  pa  shows  that  the  multiple  end-phase  configuration  is  cloarly  Inferior  to  the 
others  for  p = 4,  q « 6,  r = 11,  For  those  spacings,  multiple  end-phase  exhibits  p0  two 
orders  of  magnitude  inferior  to  cuscaded  end-phase. 

It  might  be  conjectured  that  hybrid  and  additive  midphaso  arrays  cannot  uchieve  as 
low  an  overall  pa  us  cascaded  end-phase  and  cascaded  midphase  arrays.  This  conjecture 
is  Incorrect,  as  the  following  development  will  show. 

Consider  tire  hybrid  midphase  configuration.  Suppose  that 

1.  the  array  spacings  am  transformed  according  to  the  rules 

r -*■  p'  (^ll) 
q -+  q («4) 

P -*■  r'  (-^6)  , 

but  that 


2.  the  location  of  the  phase  reference  ut  ehunnel  2 romuins  unchanged. 
Now,  by  redefining 


^ le  “ 


r'  + q'  ux  otu\ 


6 + 4 


» 5:2V 


und 
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1*2  w2 


we  can  ousily  ahow  that 


2m  r'  + *'  “ »a  <xv'2 

(■  rr+ 3 • 


a - factor  common  to  ut  and  u2 


Figure  4-B  shows  the  original,  as  well  us  the  transformed,  army. 

It  can  bo  soon,  by  referring  to  'l’ablo  4-B,  that  the  parameters  for  the  transformed 
array  aro  Identical  to  those  of  the  original  cuscaded  midphase  array,  which  in  turn  has  a 
fig  equivalent  to  the  cuscnded  end-phase  urray.  Thus,  by  redefining  the  spacing  integers 
and  the  directions  in  which  tire  subarruys  expand,  it  has  boon  shown  that  hybrid  midphase 
arrays  are  equivalent  to  cascaded  midphuse,  and  ultimately,  to  cascaded  end-phase  arrays. 


IA)  ORIGINAL  CONHICfUKATION 


(I  " 1,  ' at,' 

> M q 

r»*  T'  T Hm? 

a (ACTOR  COMMON  TO 
t,  ANDL, 

11  SPACING  COMMON  TO 
BOTH  3UI1AHMAY  RATIOS 


IUI  THANSI  OHMEO  CONITOUHATION 


SPACING  THANBFOHMATIONS: 

i|-.s=$q'.  p=>'' 


/iq'  u,  «u,' 


*’■"  7+7  • “ S7’ 

a (ACTOR  COMMON  TO 
u,  AND  v? 

r'iq'  SPACING  COMMON  TO 

BOTH  SUBARRAY  RATIOS 


KlU,  4-5-  Oi-IrIiiuI  hiuI  Iraiwfonm’d  four-elomont  hybrid  mldphuMi  array* 
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4.5  Canonical  Configurations  for  Four-element  Arrays 

Canonical  array  configurations  may  be  definod  as  those  for  which 

• The  minimum  number  of  phase  comparisons  arc  implemented. 


• The  overall  pa  is  the  lowest  achievable  for  specific  array  spacing  integers 
P,q,  r,  .... 


In  Sec.  4.5,  intuitive  arguments  showed  that  the  cascaded  end-phase  configuration  was 
optimum.  This  optimality  of  the  cascaded  end-phase  configuration  will  now  be  proved. 

It  will  bo  convenient  to  precede  the  specific  development  for  n * 4 by  some  observations 
on  arrays  using  an  arbitrary  number  of  elements. 


Consider  an  n -element  array  where  (as  shown  in  Fig.  4-1  for  Borne  four-element  con- 
figurations) the  antenna  receiver  channels  are  numbered  from  the  left  as  i = 1,  2,  ... , 
n.  Adopt  the  convention  that  for  any  olectrieul  phase  angle  <l>y^  provided  by  in- 
strumentation associated  with  these  channels,  channel  j Is  denoted  the  phase  reference. 
ThuB, 


7.* 


X <l),k  sin  °a 

O tie 

+ ,*  > ;] 

X (x/t  ~ Xj)  sin  0tt  ; 

1 

- , fc  < j 

► 0a  (+) , 

(4-20) 

where  xh  Is  the  spacing  of  antenna  k from  untennu  1,  » 0. 


Obviously,  by  definition  <l>(i  ••  0.  Further,  <l>y  * - --‘l^  y.  For  n « 2,  3,  or  4,  mutrlx 
representations  for  the  channel-pair  phuso  variables  ure 


n ~ 2: 


1 o 

© 

1 

•v 

"o  p 

t3 

o 

I 

-p  0 

(4-21u) 


n - 3: 


0 *1,2 

+1,3 

0 p 

(p  + </) 

"*1,2 

0 

+2,3 

= <l> 

-P  o 

<? 

*1.3 

-+2.3 

0 

“(p  + q)  <1 

0 

(4-21b) 
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n = 4: 


0 *1,2  *1,4 


(P  + g)  ip  + q + r) H 


0 ^2,3  *2,4 


(q  + r) 


'*1,3  "*2,£ 


-(P  + 9) 


“*l,4  ”*2,4 


-(p+q  + r)  -(q+r)  -r 


(4-21c) 


with  4>  = ^ ^ sin  “ ?r  sin  0a 


Extension  of  this  channel-pair  phase  matrix  to  n > B should  be  obvious,  It  is  em- 
phasized that  these  matrices  list  all  possible  channol-puir  differences,  but  that  usually,  a 
particular  sot  of  (n  - 1)  p*  so  differences  aro  the  only  ones  implemented  in  an  urruy  of  a 
given  configuration. 


By  inspection,  It  can  he  seen  that  the  n2  elements  in  those  mutrices  are  apportioned 


• number  of  comparisons  = n(n~l)/2 

• number  of  comparisons  = «(/»-l)/2 

• number  of  comparisons  «t»f j « n 

total  a H2 


Eliminating  the  phase  differences  which  are  not  functions  of  the  angle  of  urrivul  in 
the  systems  undor  consideration  in  this  report  leuvos  only  n(n  - l)/2  different  comparisons 
(excluding  negatives)  possible  in  an  n-eloment  interferometer.  The  number  of  theso 
comparisons  for  2 < n 7 is 


It  is  obviously  necessary  to  implement  only  (n  - 1)  phase  comparisons  in  un  m-g lenient 
interferometer,  if  these  comparisons  make  use  of  the  phase  information  from  all  n cluuinels. 


The  number  of  different  sets  of  (n  - 1)  comparisons  is  given  by  the  number  of  combi- 


nations of  n^,  taken  (n  - 1)  ut  a time,  or 
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N,„  = 


V (n-DKn^-n  + l)! 
For  2 < n < 7,  is  given  below. 


(4-22) 


2 

3 

4 
6 
6 
7 


& 


1 

3 

20 

210 

3,003 

54,264 


Note:  For  n > 3,  includes  mirror-image 
configurations;  for  strict  consistency  then, 
with  n=  2,  should  be  equal  to  2. 


For  the  matrix  form  of  Eq,  (4-21b)  provides  an  illustration  of  the  physical  counter- 

part of  N^.  The  first  tow  of  this  matrix  defines  an  end-phase  left  array,  the  middle  row 
defines  a midphase  array,  and  the  bottom  row  defines  an  end-phase  right  array. 

For  n - 4,  it  is  difficult  to  visualize  all  20  possible  arrays  with  just  the  matrix  form 
of  Eq.  (4-21c).  Therefore,  the  phase  differences  c'  ifined  by  the  matrix  form  are  used 
(»  - 1)  - 3 at  a time,  and  the  20  possible  array  configurations  are  sketched  as  in  Fig.  4-6. 
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Fig.  4-6— Simplified  ruprosunlntions  of  the  twenty  possible  four-element  iimvy 
confi|>im>Uonn  using  three  phiim*  differences 
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Table  4-6  lists  the  spacings  utilized,  the  array  phase  reference(s),  and  the  spacing  common 
to  the  two  subarrays.  The  arrays  numbered  2,  6,  13,  and  20  can  be  immediately  removed 
from  further  consideration,  since  they  do  not  employ  information  from  channels  4,  3,  2, 
or  1,  respectively. 

Arrays  7 and  12  have  spacings  of  p + q + r common  to  the  two  subarrays.  Thus, 
whichever  subarray  is  denoted  subarray  2 in  each  of  these  four-element,  arrays  has  the 
benefit  of  a large  common  factor  ex  in  increasing  the  value  of  the  argument  of  the  ambi- 
buity  variable.  Unfortunately,  the  array  ratio  for  subarray  1 will  be  so  large  (this  subarray 
ratio  also  involves  p + q + r)  that  the  pa  of  subarray  1 will  be  much  larger  than  that  ob- 
tained with  a cascaded  end-phase  configuration. 

Arrays  using  a single  common  spacing—  p,  q,  or  r— cannot,  of  course,  achieve  the  per- 
formance of  those  using  a double  common  spacing  such  as  p + q or  q + r.  This  was  shown 
in  the  calculations  for  pa  in  the  example  urrays  of  Sec,  4,4.  This  fact  removes  from  fur- 
ther consideration  arrays  3 and  5 (with  common  spacing  p),  9 and  14  (with  q),  and  16 
and  18  (with  r). 

The  remaining  arrays  share  two  attributes:  (a)  the  common  spacings  are  either 
p + q or  q + r,  and  (b)  four  of  the  arrays  are  mirror  images  of  the  other  four.  Table  4-7 
lists  the  comparison  factors  for  these  eight  arrays.  Arrays  4 and  11  can  be  shown  with 
no  difficulty  (by  developments  similar  to  those  resulting  in  Eqs.  (4-18a)  through  (4-18kJ) 
to  be  identical  in  pa  performance  to  array  1.  But  arrays  4 and  11  do  not  resolve  the 
overall  array  length  associated  with  the  integers  p + q + r as  does  array  1.  Hence  array  1 
is  superior  to  arrays  4 and  11. 

Now,  it  has  already  been  shown  in  Sec.  4.4  that  the  transformed  hybrid  midphase 
array  is  equivalent  in  pa  performance  to  the  cascaded  end-phase  array.  Once  again,  how- 
ever, array  1 la  superior  to  its  transformed  counterpart,  array  16,  in  the  sense  that  overall 
length  p + q + r is  resolved  directly. 

Each  of  the  four  arrays  discussed  above  has  its  mirror-image  counterpart  in  the  re- 
mainder of  the  table.  Thus,  it  should  be  obvious  that  with  the  aid  of  the  transformations 
p -*■  r',  q -»  q',  and  r -*•  //,  the  performance  of  second  group  of  arrays  is  identical  to  the 
first  group,  and  that  a cascaded  end-phase  configuration  is  indeed  optimum. 

In  contrast  to  two-element  mid  three-element  arrays,  the  performance  of  arbitrary- 
configuration  four-element  arrays  is  dependent  both  on  the  configuration  and  on  the 
array  expansion  direction. 

The  preceding  development  has  shown  that  of  the  20  possible  (by  definition)  four- 
element  array  configurations,  there  are  four  configurations  that  have  particularly  simple 
realizations : 

• Amy  1— Cascaded  end-phase  (Reference  channel  1) 

• Amy  8— Transformed  hybrid  midphase  (Reference  channel  2) 

• Array  16— Transformed  hybrid  midphase  (Reference  channel  3) 

• Amy  19— Cascaded  end-phase  (Reference  channel  4). 
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Table  4-7— Comparison  Factors  for  Four-element  Arrays  with  Spacings  b + q otq  + r 


Array  Configuration 


Common  Array 


'Mirror-image' 


Spacing  Reference(s)  Number  Number 


Cascaded  end-phase 
Equivalent  to  CE-P 
Equivalent  to  CE-P 

Transformed  hybrid  midphuse 


p + q 
P + q 
p + q 


p + q 


Axis  of  configuration  symmetry 


Transformed  hybrid  midphase 

Cascaded  end-phase 
Equivalent  to  CE-P 
Equivalent  to  CE-P 


q + r 

q + r 
q + r 
q + r 


Note!  Amiy  numbers  ure  defined  In  Fig,  4-6, 


Those  four  array  configurations  are  defined  by  rows  one,  two,  three,  and  four  in  the 
channol-phaso  error  matrix  of  Eq,  (4-21c). 

This  development  showing  that  the  cascaded  end-phase  configuration  is  indeed  opti- 
mum for  n - 4,  is,  as  far  as  the  writer  can  determine,  a new  result.  Although  some  of 
the  materiul  presented  by  Hunson  [21]  for  the  case  of  four  apertures  on  a line  uppears 
relevant,  he  neglects  to  consider  the  effects  of  correlation  betwoen  channel-pair  phase 
errors  in  aubarrays.  Thus,  he  is,  in  effect,  formulating  judgements  on  array-configuration 
suitability  bused  on  “independent,”  i.e.,  four-element  two-integer  set  suburruys.  It.  is 
known  from  Sac.  3.0  that  this  leads  to  incorrect  conclusions  when  one  deals  with  three- 
element  subarrays. 

This  concludes  the  exposition  of  the  fundamental  analyses  associated  with  four- 
element  arrays.  Sec.  6.0  will  address  the  synthesis  of  four-element  urruys  and  will  provide 
techniques  by  which  the  optimum  array  of  any  length  K ~ p + q + r half-wavelengths  cun 
be  readily  synthesized. 


5.0  SYNTHESIS  OF  OPTIMUM  FOUR-ELEMENT  ARRAYS 

This  section  provides  the  theoretical  busis  for  synthesizing  the  optimum-realizable 
four-element  cascaded  end-phase  array  of  arbitrary  length  % > 4 half-wavelengths.  (Noto: 
the  case  S ~ 3 is  trivial;  the  urray  consists  of  four  channels  separated  by  spacings  p = q = 
r - 1 half- wavelength.) 
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It  is  shown  first  that  the  ideal  numerators  to  employ  in  the  subarray  ratios  'K,.  = 
nti/tij  are  just  = >/£  Of  course,  this  condition  can  be  met  only  for  those  cases  in 
which  1!  is  a perfect  square.  Next,  an  ideal-unrealizable  array  (achieving  the  minimum 
possible  pa)  is  defined.  The  only  four-element  array  for  which  this  ideal  condition  can 
be  met  is  the  one  defined  by  p « 1,  q = 1,  and  r = 2;  8 = 4.  The  performance  of  an  array 
of  any  other  length  can  approach,  but  never  equal,  the  performance  of  the  ideal-unrealizable 
array  for  that  length.  Interestingly,  Moffot  [22] , citing  a result  due  to  Arsac,  shows  that 
the  longest  zero-redundancy  array  in  spatial-frequency  terms,  from  the  radio-astronomy 
art,  regardless  of  the  number  of  elements,  cannot  exceed  8-6,  with  n = 4 elements. 

The  concopt,  of  the  ideal-unrealizable  array  is  quite  useful,  nevertheless,  because  It 
leads  to  a direct  approximate  synthesis  procedure.  This  initial  procedure  assumes  inde- 
pendence of  subarray  ambiguities  but  achieves  the  twin  objectives  of  (a)  readily  visualized 
physical  spacing  interpretations,  und  (b)  rapid  convergence  to  a region  of  low  pa  so  that 
an  exact -pa  analysis  on  only  a very  small  set  of  candidate  arrays  need  be  done. 

The  synthesis  procedures  given  in  this  section  lead  to  a specification  of  the  channel 
spneings  for  the  optimum-realizable  (minimum  physically  realizable  probability  of  ambi- 
guity) urray  of  arbitrary  length  under  the  assumption  of  channel-pair  phase  errors  that 
are  Gausslun,  zero-moan  with  equal  standard  deviations.  This  assumption  simplifies  the 
analysis  und  subsequent  synthesis,  and  in  practice  is  not  too  restrictive.  This  is  because 
the  phuac  errors  in  channels  comprised  of  several  components  in  cascade  (see  Sec.  3.1) 
upproach  statistical  regularity. 

For  a procedure  to  be  exact,  it  must  of  course,  take  uccount  of  possible  dependence, 
or  correlation  between  the  ambiguities  in  subarray  1 und  those  in  subarray  2.  This  de- 
pendence is  indicated  by  tho  absolute  magnitude  of  on  arruy-to-array  correlation  coeffi- 
cient. For  equal  chunnel-pair  phase  errors,  the  correlation  coefficient  is  u function  only 
of  tho  subarray  ratio  integers  (or  equivalently,  the  array  spacing  integers).  A knowledge 
of  this  correlation  coefficient,  in  conjunction  with  the  umbiguity  variables  previously  de- 
fined in  Soc.  4.0,  enables  one  to  define  the  complete  bivuriuto  umbiguity  density  function 
(actually,  the  no-ambiguity  density  function).  Integrating  this  function  between  the  ap- 
propriate limits  provides  the  probability  of  “no  ambiguity,”  und  ultimately,  the  probability 
of  umbiguity. 

For  many  casos  of  interest,  the  array  pa  is  only  a weak  function  of  the  urruy-lo-urruy 
correlation  coefficient.  The  value  of  tho  exact  formulation,  including  the  lurray-to-urray 
correlation  coefficient,  is  thut  a four-element,  arruy  of  any  configuration,  subject  to  vurious 
channel-pair  phase  errors,  can  be  precisely  characterized.  For  example,  arrays  in  which 
more  than  ono  antenna  arc  switched  sequentially  into  u common  channel  (to  minimize 
hardware)  can  be  exactly  analyzed  upon  derivation  of  the  appropriate  suburruy-to-suburray 
correlation  coefficient. 

In  the  final  portions  of  this  section,  the  urray  length  is  characterized  for  synthesis 
purposes  as  fulling  into  one  of  three  classes.  Thoso  are  Class  I (8  = perfect  square),  Class  II 
(8  =>  geometric  mean  between  successive  square  arrays),  and  Class  III  (8  = any  length  not 
falling  into  Classes  I and  II).  Examples  of  an  optimum-reulizuble  synthesis  for  each  of 
thoso  array  classes  ure  given. 
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The  section  closes  with  tabulated  exact  pa  (including  the  effect  of  array-to-array  cor- 
relation) for  optimum-realizable,  four-element  cascaded  end-phase  arrays  of  various  lengths 
from  8 = 4 to  100  half-wavelengths,  with  channel-pair  phase  error  as  a parameter.  It  may 
be  noted  that  very  early  multielement  interferometer  configurations  baaed  on  a binary 
progression  of  spacings:  p --  l,p  + <?-2,  p + t/+r=4,  p + q + r + s=8,  etc.,  [23] , 
inherently  have  a very  high  resistance  to  ambiguities  for  the  usual  range  of  channel-palir 
phase  errors  (even  when  relatively  coarse  quantizing  of  electrical  phase  is  done).  Thus, 
based  on  the  phase-tracking  performance  associated  with  presently  available  microwave 
components,  such  a multielement  design  approach  represents  extreme  (and  uneconomical) 
overdesign. 


6.1  Optimum  Numerators  m(  for  the  Subanay  Ratios  (R,  in  n-Element  Arrays 

In  synthesizing  n -element  arrays  of  a given  length  8,  it  is  desirable  to  have  a proce- 
dure for  obtaining  the  particular  set  (or  sets)  of  subarray  ratioR  JRj,  i = 1,  2,  . . . , n - 2 
that  will  provide  the  lowest  overall  pa,  subject  to  given  cCH  pr.  (The  optimum  31 , for 
n a 3 elomont  arruys,  has  already  been  formulated  in  Sec.  3.0.) 


Since  the  urruy  length  8 = mt  m2/a  in  four-element  arrays,  It  is  natural  to  focus  at- 
tention on  the  numerators  m(  in  the  subarray  ratios  !fl(  = m^.  For  four-element  cascaded 
end-phase  arrays,  the  optimum  (in  general,  unrealizable)  and  31 2 are  readily  derived. 

The  generalization  to  n > 6,  i.o.,  three  or  more  subarruy  ratios  !R/,  follows  oaaily. 


Making  use  of  Eq.  (3-19)  and  the  summary  materiul  in  Eqs.  (4-18a)  through  (4-181) 
results  in  the  overall  probability  of  “no-ambigulty,”  i.e,,  probability  of  "correct,"  Pc.0vBr«ll> 
in  u four-element  cascaded  end-phase  array,  of 


Pc-ovorull  “ (l“Poi  W Pu'j) 


r«(rr/«0i 

exp  (- tf/2)  dt[  , 


(5-1) 


where 

= minl  + * 

cti  - 1,  i ~ 1. 

= «,  1-2,  ot  r-  factor  common  to  und  n%  . 

Equation  (6-la)  assumes,  of  course,  that  ambiguities  in  subarray  2 (resolved  by  sub- 
urray  1)  are  independent,  of  ambiguities  in  subarray  1.  In  Sec.  5.4,  an  additional  develop- 
ment will  include  the  effects  of  correlated  suburruy  ambiguities.  It  will  be  seen  that  op- 
timum suburruy  ratios  arc  only  weak  functions  of  suburruy  correlation  coefficients. 
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It  is  convenient  to  normalize  the  n,-  in  Eq.  (5-1)  to  the  associated  m;,  as 


and 


2-|l/2 

s ”1  x (n'\ 

“ miwcH-PR  [i  - sq  + 


302 


.2]1'2. 

a = riy/mi  , 

0 < a < 1 , 

(5-2a) 

6>] 1/2 . 

b « n2lm2  , 

0 < b < 1 . 

(5-2b) 

From  Eq.  (4-18c),  the  array  length  6 is 


2 = p + (j  + r- 


For  2 constant,  and  with  subarray  ratio  optimization  carried  out  with  respect  to  mj  (for 
convenience),  w2  is  given  by 


«S 

m2  " mT 


(6-3) 


Thus,  with  Eq.  (5-3),  the  integration  limits  (±)  in  the  integral  forms  of  pCj  in  Eq. 
(5-la)  become  ’ 


ir 


Subarray  1 


gCH-PK 
wt(l-a  + oa)1/a 


(5-4a) 


and 


Subarruy  2 


CV7T 


air 

°CHPR 


m.,(l  - b + 62) 


1/2 


m 


7T 


1 0 


CH-PH 


2(1  -6  + f»2)1,a 


(B-4b) 


Functions  of  the  form  (1  - x + x2)  z,  0 < x < 1,  huve  a single  minimum  of  y/0.75-- 
0,866,  ut  x - 0.5,  compared  to  en.  point  values  of  1.000  at  x - 0,  1.  Independent  selec- 
tion of  a - and  b n2/m2  is,  in  general,  not  possible  because  m j and  n2  are  both 

related  to  the  same  array  spacing  p *■  q.  But,  to  establish  optimum  realizable  mjt  we  can 
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assume  that  a and  b ~ 0 or  1,  with  negligible  error.  The  integration  limits  of  Eq,  (5-4a,  b) 
are  now  approximated  by 


Subarray  1 


It  _ °CH-PR 


(6-Bo) 


Subarray  2 ~ 


mi^o 

% * * ’ 


(6-6b) 


where 


u s 1! 

0 ffCH-PR 


If  m-i  a both  integration  limits  in  Eq.  (6-6)  are  equal,  and  Pc.0vori»il  “ Pcx  ~ P'ty 

As  in  numerous  related  problems  in  communications  technology  in  which  there  is 
freedom  to  vary  parameters  to  minimize  overall  probability  of  error,  a value  of  \/$  is 
readily  shown  by  assuming  that  m\  = m i-opi’  + ^ml  Qctually  the  optimum  value  and 
then  proving  the  converse. 

Assume  that  A - 6 « 1.  Then,  in  Eq,  (6-5),  the  integration  limits  are 


Suburray  1 J-  - ^ d " «> . 


(6-6a) 


Su'oarray  2 “ ^ + 6) 


(B-6b) 


For  small  changes  In  integration  limits,  euch  pc  in  Eq.  (5-la)  can  bo  expressed  as 


. ( ko\  6,{o  2 f *0 


(6-7u) 


where 


/ /i(j\  bkQ  2 / fco\ 

1>)  ' V*  exp  V>  ’ 


(6-7b) 
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‘l>(x) 


exp 


s/%r 

The  overall  pc  thus  becomes 


dt. 


Pc-ovurull  w 


(5-8) 


The  overall  probability  of  correct  ambiguity  resolution,  pc  .ivmu,  *8  maximum  when 
6 = 0,  proving  that  the  optimum  value  for  mi  is 


ml  -OPT  = V^T,  n ~ 4 elements.  (B-9) 

The  generalization  for  arbitrary  n > .5  is  obviously 

mi-oPT  “ i\  m > 3 elements,  (5-10) 

taking  into  account  that  an  n-element  array  defines  (n  - 1)  clement  spacings,  and  u mini- 
mum set  (for  ambiguity  resolution)  of  (m  2)  subarray  ratios  !Kj  and  associated  m(. 


5.2  Integration  Limits  in  Ideal-Unrealizable  and  Realizable  Arrays 

Determination  of  particular  spacings  for  the  lowest  pa  in  a four-clement  array  of  a 
given  length  K can  bo  accomplished  by  exhaustively  iterating  Kq.  (5-1)  through  all  possible 
nij  and  n(  defined  by  />,  </,  and  r,  and  selecting  the  nonrodundant  array  thut  exhibits  the 
lowest  pu.  This  determination  is  not  necessary  if  the  concept  of  the  ideal-unrealizable 
array  is  utilized. 

Consider  u four-element  array  in  which  the  length  for  generality,  Is  not  a perfect 
square.  ICquution  (5-9)  from  the  previous  section  defines  the  ideal  m j ,in  ~ R*2-id  \Z^> 
which  is  not  an  Integer  if  i!  A I2,  i > 2.  If  uny  considerations  of  multiple  ambiguities 
over  the  field  of  view  are  disregarded  (because  of  common  factors  in  the  mj  and  rt,  sub- 
array ratio  integers),  the  ideul  (ID)  n(  associated  with  m(  are  just 

"•/■ID  „ 

bMD  " 2 “ 2 ' 

Then,  the  integration  limits  (II,j)  in  Kq,  (5-1)  for  pc  for  these  subarruys  are 
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r 


f:. 


j 


V 

*■ 

[I, 

f.. 

f 


ILj.jd  - ILa-m 


ir 1 

°CH-PB  - m,,D„,1D  ♦ ,f.IDl>'2 

•_  J_  1 

/TJ71 


h.JL 


(5-11) 


where 

j,  H_ 

0 " °CHPR 

An  example  of  the  use  of  this  concept  is  an  array  whose  length  is  12.  By  definition, 
'”l-ID  K tfitf-iD  = \/l^  * 3.4641  . . . ; = rt2.jD  --  1.7321  ....  No  p,  q,  and  r exist 

such  that  p + 7 + r = 12  and  (p  +■  q)/p  = 2:1,  (p  + q + r)/(p  + q)  = 2:1,  If  such  on  array 
existed,  it  would  exhibit  the  lowest  pa  of  ail  arrays  of  length  12.  Hence,  the  term  ideal- 
unrealizable  array  seems  appropriate. 

Arruys  which  ure  realizable,  said  which  approach,  but  do  not  necessarily  equal,  the 
performance  of  the  ideal  array  are  called  “optimum-realizable.’' 

Realizable  (RE)  ami'  have  un  integration-limit  (IL)  product  given  by 


1 


where 


( 

{mf  - m;rt(  + rtf  )1/2 


«l  = 1, 

«2  ra  factor  common  to  mj  and  «2, 


*o  ~ 7T~^ — ' w in  Kcl*  (5’1D* 
°CH  PR 


(6-12) 


Kuch  of  the  IL/.HE  can  bo  expressed  us  the  sum  of  ILMD  and  a deviation  from  this 
ideal-unrealizable  integration  limit,  us  follows.  For  we  have 


b.0  — ■ 

1 

(613a) 

(mf  - mt 

«l  +nf): 1/2 

2fe0 

>/3? 

2(mf 

- mjMj  + rtf)1/2_ 

(6-13b) 

2*o 

7M1 

[1  + 8 

il. 

(6  13c) 
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Equating  the  terms  within  the  brackets  in  Eq.  (5-1 3b)  and  (5-13c)  defines 

6 s <374  ■ _ 1 _ 

(mf  - mj/ij  + 


(5-14) 


If  the  form  dj  = 2fco /s/3K  * $i  is  used  in  conjunction  vith  Eqs.  (5-13)  and  (5-14),  it,  is 
easy  to  express  ILj.rb  in  the  form 

,r  2feo  , , r i 2 

IL,“  “755  *•  -?n 


“ ili  id  + <*1  • 


(5-15) 


By  a similar  development,  it  can  be  shown  that 

2/?o  ix  2 

ils»-re  “ ~m  + ko  " ITxPi  7 fj? 

v3®  (mi  - m2n,j.  + nj|)1/2  v3® 


IL2.id  + d2 , 


(5-16) 


The  only  realizable  four-clement  array  for  which  IL/.ke  “ ^Ji  ID  ia  the  one  in  which 
p « 1,  q » 1,  r * 2;  fi  « 4,  Here,  « m2  » 2,  **  rt2  » 1,  a « 1.  Thus, 


d‘  ° (22  - 2 < 1 + 1st  » 3),/Z  73  * 4 


" ''»[73  ' T?]  " ° 


t - 1,2, 


and  hence, 


lln.RK  ==  IL2  RE  ~ •"  il2  ID  “ 


vW " 


(,CH-1'R. 


Other  urrays  approach,  but  do  not  satisfy,  the  relations  ® II^.jD  because  the 

rutio  miini  can  equal  2 only  for  :R  / & 2:1.  The  d(  in  Eqs,  (5-15)  and  (5-16)  are  u measure 
of  how  a realizable  array  differs  Iri  performance  from  un  ideal-unrealizable  array.  Synthesis 
procedures  bused  on  the  concept  of  deviation  from  the  condition  of  ideul-unreulizabillty 
will  bo  given  in  the  next  subsection, 


'-lA  ... 


ROBERT  L.  GOODWIN 


5.3  Synthesis  of  Cascaded  End-phase  Four-element  Arrays 

The  basis  of  the  synthesis  procedures  to  be  discussed  involves  the  subarray  1 numer- 
ator and  manipulations  of  the  two  realisable  intonation  limits,  ILj.re  and  1L2.re* 
The  standard  forms  for  these  two  limits  are 


und 


ilmie  " flo 


1 

(mf  - mxnx  +n?)l/2 


1L2-RE 


tv 


(m§  - m2n2  + nj) 


1/2 


(6-17a) 


(5-17b) 


In  this  form,  the  dependence  of  the  integration  limits  on  array  longth  V.  is  implicitly, 
rather  than  explicitly,  expressed.  An  equivalent  form  for  IL2.WE  >H  more  suitable  for 
array  synthesis. 


Use  of  Eq.  (4-18c)  from  Sec.  4 yields  equivalent  forms  for  m2  and  n2; 


und 


n» 


(p + (/)a 


m, 


When  these  substitutions  uro  made  in  Eq.  (B-17b),  and  some  simplification  is  done,  the 
result  is 


IL2-RE  - 


(B-1H) 


The  channel  1 to  channel  3 spuoing  integer,  p + q , cun  be  equal  only  to  a multiple  of  the 
subarruy  1 numerator  integer  mx  us, 


l>  + q r:  /MM  • 


Furthermore,  p + q cannot  exceed  V - p + </  + r,  Thus,  the  1E/.rk  Ur>  u form  suitable 
for  synthesizing  realizable  urruys)  become 


und 


ll'l-RK  ~ *o 


(m i - + /if) 


2\l/2 


(B-19a) 


- k 


Mil 


" d'2  K(/f»,)  + (jMll)a),/8 


(B-19b) 
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, where  *]  indicates  the  nearest  integer  less  than  x . 


So  far,  nothing  more  has  been  accomplished  than  to  express  the  IL/-H.B  ’n  a slightly 
different  form  to  emphasize  that  subarray  2 is  more  logically  thought  of  as  a function  of 
the  overall  array  length  1!  und  the  subarray  1 numerator  integer  in ^ . If  one  wanted  to 
establish  that  a particular  set  of  mh  nt  is  the  optimum-realizable  for  given  8,  integration 
of  Guussian  error  functions  (as  in  Eq.  (5-1))  seems  necessary.  Those  integrations  can  be 
uvolded  by  the  U6e  of  a simple  test  function  which  Involves  the  expansion  of  the  error 
function  about  a specific  value. 

The  first  step  in  determining  this  test  function  is  to  express  Eq.  (5-19b)  in  the  form 
of  Eq.  (5-16);  that  is,  as  the  sum  of  IL2.RE  und  a term  d2  which  represents  the  deviation 
of  ILo.re  from  IL2„m.  This  is  reudily  done.  The  result,  with  Eq.  (5-15)  for  IL1>RK  re- 
peated below  for  ease  of  comparison,  is 


2fep 

,Ll-RK 

“ y/H  + 

*o 

- IL1-1D 

+ d 

2 Zip 

IIj1!-RE 

= 7S? + 

* 0 

= 1L2-IO 

+ </; 

(m'f  - 


(5-20u) 


»ii 


- mmt)  + Um,)2) 1/2 


(6-20b) 


Equation  (A-9)  of  Appendix  A gives  u throe-term  Taylor’s  Hertes  expansion  of  the 
Gausslun  error  integral  about  a specific  argument,  as 


* (-H » 'i  (f)  - z [2  (!-)  (#f  (•’-)  * t (f)  (fj 


(5-21) 


whor« 


\ - jmXI  iJ 

V sf'h flxl0  c 


xln 
7 < 


Z 


For  x/o  2.000  (i.e.,  the  argument  of  A(-)  is  set  ut  the  2«  value).  A(x  + A/o)  becomes 


- £) ' M2)  + [a  l«) 


toJ 


(5-22) 
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Equation  (5-22)  is  an  approximation  to  the  Gaussian  error  integral  in  the  vicinity  of  an 
argument  set  at  the  2a  value.  Note  that  A( 2)  = 0.9544997.  The  corresponding  compound 
pa  (assuming  independence)  is  1 - A(2)A(2)  = 0.08893  or  8.893^. 

Now,  k0  is  defined  as  ?r/ocl;J. PR.  Thus,  to  force  IL,-.ID  to  assume  the  value  of  2.000 
(strictly  for  simplicity  in  generating  a useful  test  function)  regardless  of  array  length  E, 

Eq.  (5-11)  can  be  set  equal  to  2.000  and  solved  for  the  appropriate  value  of  aCH.pR,  as 


1L, 


/•ID 


2n 

y'SE 

°CHPR 


= 2.000, 


or 


°gh-pr12  = ^ • 

Use  of  this  value  of  o^h-pr  in  Eq.  (5-20a)  and  (5-20b)  yields 
ILl-RE  = 2 + 


(5-23) 


and 


« 2 + 


IL2-RE  ~ 2 + 


= 2 + 


]fl„ 


3E 


m\  - injii]  + n\ 


(5-24a) 


mj  >/3E 


V^2  - Wmi)  + Umj* 

(*)■ 


(6-24b; 


The  approximate  probabilities  that  subarray  1 and  subarray  2 (resolved  by  subarray  1)  are 
unambiguous  are 


pc..[  =*  0.9644997  i 


s/2rr 


,~2 


. . . . 2 , . 3 


and 


pc..,  * 0.9544997  + e~2 


/A?\ 

/a2\5 

» 

' u 

e.\ 

2 M 

- 2 — 

+ - 

~ 

\ o / 

\o  / 

\ 

o 1 

(5-25u) 


(5-26b) 


The  ])rocedure  for  obtaining  an  approximate  array  synthesis,  one  not  taking  into  con- 
sideration subarray-tc-subarray  correlation  of  ambiguities,  can  be  summarized  ac  follows. 
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• Choose  ■>  trial  m . ( to  ./fi„ 

• Choose  an  appropriate  value  of  i j from  Table  3-4  of  See.  3,0. 

• lLq.RE  of  Eq.  (5-24a)  is  thus  defined. 

• Compute  a set  of  trial  values  of  IL2-RE  ’°V  iterating  through  Eq  (5-24b)  with 
J - 1,  2,  . ..,  S/m  1]  where  •]  indicates  the  nearest  integer  less  than  the  ratio  fi/m*. 

• Substitute  the  trial  IL1-RE  and  the  set  of  trial  IL2.RE  into  Eq.  (5-25a)  and 
Eq.  (5-25b)  und  compute  the  overall  yrobability-of-correct  product  PClPC2- 

• The  trial  anay  achieving  the  highest  pc  pt.  product  is  the  best  array  for  the 
particular  trial  mj  that  was  chosen  initially. 

These  steps  cun  be  repeated  for  other  trial  values  of  mlt  and  the  best-performing 
arrays  obtained  for  each  rn^  can  tlien  be  compared  to  determine  the  optimum-realizable 
array.  Ibis  approximate  procedure  does  not  take  into  account  a possible  subarray-to- 
subarray  correlation  of  ambiguities.  Thus,  several  arrays  may  bo  obtained,  all  of  which 
have  the  same  independent  probability  of  ambiguity.  Also,  since  this  approximate  proce- 
dure is  based  on  an  expansion  of  the  normal  probubility  integral  (evaluated  in  the  vicinity 
of  argument  2.0000),  a slightly  different  arruy  may  actually  have  a somewhat  better  pa  at 
a lower  value  of  aClI.PR  (implying  a larger  value  of  the  probability  function  integral  argu- 
ment) than  the  array  found  by  this  procedure.  The  designer  may  have  to  perform  an  exact 
analysis  at  the  specific  vulue  of  och-PR  that  will  be  used  for  a particular  application. 

Both  the  procedure  outlined  above  and  an  exact  analytical  form  for  pa,  are  easily 
programmed.  For  the  convenience  of  the  designer,  these  programs  are  collected  in  Appen- 
dix B.  The  use  of  these  programs  will  be  described  fully  in  Sec.  B.5. 

In  Sec.  5.4,  on  exact  relation— one  which  considers  the  effect  of  subarray-to-subarray 
correlation  of  ambiguities— for  the  pa  in  any  four-element  array  will  be  derived.  This  ex- 
pression will  become  the  basis  for  determining  the  optimum-realizable  array  (from  pa 
considerations)  in  the  cascaded  end-phase  configuration  (from  the  development  of  Sec.  4.5) 
of  any  length. 


5.4  Exact  pa  for  Four-element  Arrays 

In  the  previous  section,  a synthesis  procedure  was  given  for  cascaded  end-phase  four- 
element  arrays  of  arbitrary  length  S > 3 (in  half-wavelengths).  The  procedure  was  approxi- 
mate, as  the  possibility  that  ambiguities  in  suburray  1 might  be  correlated  with  ambiguities 
in  subarray  2 was  ignored.  As  will  be  shown  in  Sec.  5.5,  certain  arrays  having  the  sarno 
length  1!  apparently  have  the  name  probability  of  umbiguity,  even  though  the  element 
spacings  p,  q,  and  r are  different  if  correlation  is  ignored.  In  actuality,  there  is  a slight 
dependence  of  pa  on  array  arrangement  for  fixed  l£. 

This  dependence  is  most  conveniently  expressed  in  terms  of  an  array-to-array  correla- 
tion coefficient.  This  correlation  coefficient  p(Jl  u will  be  shown  to  be  a function  of  the 
array  configuration  (e.g.,  cascaded  end-phase,  hybrid  midphuso,  etc.)  as  well  us  array 
spacings  (e.g.,  actual  p , q,  and  r employed). 
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The  ambiguity  constraints  of  a four-element  cascaded  end-phase  array,  using  modi- 
fied forms  of  Eq.  (4-18d,  e),  can  be  written  as 

Subarray  1 -ir  < [</>i  = n1Aifi13  - < ir , (6-26 


Subarray  2 -it  < [02  - «2^j,4  ” m2^i,sl  < n • (6-26b) 

The  above  two  random  variables  (t.v.)  '!>,  and  'b2  are  Gaussian,  since  they  are  linear  com- 
binations of  various  channel-pair  phase  errors  Ai pX  j ~ At - A = 2,  3,  4)  which  are 
assumed  Gaussian,  zero  mean  according  to  the  development  in  Sec.  3.1.  The  variances  of 
these  ambiguity  variables,  if  Eq.  (3-11)  is  used,  are 

<4  = ^i0?,2  " tPctHhinvi^oi,*  + "101,3  (6-27a) 

and 

4 = ">20?,3  - 2Pc2"’2"201,3al14  + "201,4  . (B  27b) 


01,3  = 


0?  + 02 
01  + 03 
a‘(  + o'i 


These  define  relations  between 
channel-pair  error  variances  and 
channel  error  variances, 


P#,  = d77r77"  <==0-6>  a11  cclual>  • 


P*a  = (=0.6,  all  Oj  equal)  . 

If  we  allow  for  possible  correlation  between  the  r.v.  <I>,  mid  <I>2,  their  joint  density 
function  is  [18] 


P<lq  «|>2(PlV?2) 


2^00!  002  >fi  - 


X exp 


<„,)  L\“ 


2 2 
_f2_  + rfA 
00  J 00,  002  WJ 


The  suburray  1-t.o-subarray  2 correlation  coefficient  J)0  is  defined  by 
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£[<«*,  *l>a] 


°l.°a  Orf,  u, 


*l°*t 


(5-29) 


The  probability  that  I <^>i I < ir  and  |^2I  = ® ?r,  i.e.,  the  probability  that  array  ambi- 
guities will  be  correctly  resolved,  is 


J7T  r>an 

-v  '-an 

Equation  (6-30)  becomes  just 

fir  fan 

/><[*, (<Pl)  d*T  P<|»2(^2)d^2  . 

TT  rVTT 


(5-30) 


(5-31) 


lor  Manila  = °- 

In  general,  pfll|U  ^ 0,  and  to  calculute  the  exact  pc.()VWul)  requires  (a)  determination 
of  f>ahan  and  (b)  integration  of  the  bivariate  density  function  of  Eq.  (5-28)  between  the 
appropriate  limits. 

The  array-to-array  correlation  coefficient  is  readily  found,  if  Eq.  (6-29)  is  used  with 
Eq.  (6-27a,  b).  The  expected  value  of  (‘I’j CI>«2 ) — is  given  by 

tf[<l>i.<l>2l  = #K«iAiPi,3  - »«lA«Pi)2)(n2A^ii4  - w2Aipli3)] 

= «iM2E[(Av>i  - A^aKAv?!  - Aip4)] 


“mjMjjE’KAipt  “ Av»2)(A^  - Av?4)] 
-n,  m2/i[(Av?i  - A^3)zl 


+mlm2E[(Av>i  - A^HA^j  - Av>3)]  . 


(6-32) 


I 

The  taking  of  the  expected  values  of  the  vurious  cross-products  of  channel-pair  orror  yields 

| E[<lyl>2]  - nin2o'x  --  mln2o'i  - n,  tn2(o'f  + o$J  + mtn2oz  . (6-33) 

Equation  (5-33)  expresses  E[<I>1‘I>2]  in  channel-error  form.  An  equivalent  expression,  in 
channel-puir  orror  form  is 


E[‘I’(  ‘1’21  " «in2Pl,3;l,40l,8c,l,4  ~ mlrt2 Pi, 2 ; 1,4°1,4 


'nlm2Pl,:»;l,;»°f,3  + m1n,2Pl,2;l,aol,20l13  » 


(6-34) 


where 


* 
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PlJll.ft 


£f(Awi  - Av?;)(A <p\  - Av?fe)] 


°U°U 

By  definition,  Pi ,351,3  = 1.  Now,  if  all  channel  phase  errors  ai  are  equal,  this  im- 


By  definition,  p,,3.13  = 1.  Now,  if  all  channel  phase  errors  0;  are  equal,  this  im- 
plies that  of  = 0.5  ofj  * 6.5  ofjji -PR  and  yiat  p<j; l k = 0.5,  j * k.  For  this  situation, 
both  forms  of  £[(I»1<l’'2l  become 


£[‘I>i<!>2]  = orcH-PRf(«in2  " min2  +mlm2>  “ 2n!m2]  , all  of  - OcH-PR  • 

(5-35) 


Equation  (5-35)  can  also  be  expressed  as 


\ och-trK"1!  - «i)(m2  '■  Ma)  - • a11  °f  ~2  °CH-PR  - 


(5-36) 


Finally,  the  use  of  Eq.  (5-36)  with  Eq.  (5-29)  yields 


JL 

P<ii,fl2  = 2 


(mj  - ri})(m2  - n2 ) ~ ntm2 


(m?  - + nf  - m2n2  + n!)1/2 


, all  of  ■=  — o£h-PR  . 


(5-37) 


Equation  (5-37)  is  solely  u function  of  the  subarray  ratio  integers  for  of  - 0.5  ofm-PR* 
Furthermore,  as  with  uny  correlation  coefficient,  -1  < Pai,oa  < 1> 


Array-to-array  correlation  coefficients  for  the  other  four-element  array  configurations 
of  Fig.  4-1  and  Table  44  are  readily  derived.  Tuble  5-1  lists  the  correlation  coefficients 
for  all  five  array  configurations. 


It  is  convenient  to  define  two  standardized  variables  with  standard  analytical  tech- 
niques for  integrating  bivariate  density  functions,  as 


<l>n 


‘i’l  _ l,'2 

* ^ 1 tJ  ~ % 


Equation  (5-30)  is  thus  transformed  into 

4/j  /%fi 


Pc-ovorall  ” f ^ dx 


L 2(1 -p2) 


(x2  - 2 pxy  +ya) 


dy 


f * f* 

J-h  J-u 


g{x,  y,  p)  dy  , 


(5-38) 


where 
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e 5-1—  Subarray  to  Subarray  Correlation  Coefficients  for  the  Four-element  Array  Configurations  of  Fig. 


* All  Cj  = t ffcH-PR  > • = I to  4.  t After  subarray  1 is  resolved  and  is  determined  (see  Fig.  4-1). 
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P Pax  ,au  • 

Equation  (6-38)  is  the  standardized  form  for  the  integral  of  the  bivariate  Gaussian  density 
function  [19,24,25], 

An  alternate  form  for  Pc.ovcrun  which  uses  the  well-known  L{h,  k,  p)  and  /!(*)  func- 
tions [19a] 


JVov«rall  3 2L(A.  A,  p)  + 2L(h,  k,  -p)  + A(h)  + A(k)  - 1 , (5-39) 


whore 


L(h,  k,  p) 
A(h) 
A(k) 


oo 

g(x,  y,  p)dx, 


■*2/2  dt  - Probability  of  being  correct,  subarray  1, 


~t2l2  dt  - Probability  of  being  correct,  subarray  2 (resolved 
by  l). 


Many  expressions  for  calculating  L(h,  It,  p)  are  available;  one  which  is  suitable  for 
/;,  h > 1 and  Ipl  < 0.95  as  is  the  case  for  arrays  of  interest  in  this  report,  is  based  on  the 
series  expansion  [19b] 


L(h,  k,  p)  = Q{h)Q(k)  + 

fi"0 


(ii  + l)! 


(5-40) 


where 

Z(n){x)  = -.jL  Z{xh 
dxn 

Z(n  + *)(x)  = -*/(«+')(*)  - (n  + l)Z(")(x). 
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A special  value  of  L(h,  k,  p)  is 

L(h,  ft.  0)  = Q(h)Q(k) 

= | |1-  A(h)Ul  ~Mk)]  . (5-41) 


This  indicates  statistical  independence  between  ambiguities  in  subarrays  1 and  2,  for 
p~  0.  Thus, 


Pc-ovorall  " [l  ~ A(h)][l  ~ A(k)]  + A(h)  + A(k)  - 1 
- A(h)A(k). 


(5-42) 


Obviously,  the  overall  probability  of  ambiguity,  p0.ovt,ri,ii  is  1 - Pc-overull  from  Eq. 
(5-39),  and  is  given  by 


JVovcrall  “ 2 - [2L(/i,  ft,  p)  + 2/4h,  ft,  -p)  + A (ft)  + /!(/<)] 


(5-43) 


where 


h =■• 


% ’ 


Qtlf 


°<h  ' 
p = ’ 

and  the  functions  L(h,  ft,  p),  A(ft),  A(ft)  are  as  defined  by  Eq.  (5-39), 

The  following  points  may  be  noted  about  Eq.  (5-43). 

• /Vovorall 18  independent  of  the  arithmetic  sign  of  p,  since  the  sum  I Ah,  k,  p)  + 
L(h,  It,  -p)  is  independent  of  the  sign  of  p, 

• Pa-overall  *s  greatest  (other  parameters  being  equal)  when  p is  zero,  since  L(h.  k,  p)  + 
Hli,  K p)  > 2 L(h,  It,  0). 


• The  formula  for  pa  is  general  and  can  bo  used  to  calculate  tho  performance  of 
any  four-element  array,  where  ft,  It,  and  p are  properly  defined. 


5.5  Tabulated  pa  and  Array  Spucings  for  Optimum-Realizable  Four-Element  Cascaded 
End-Phase  Arrays 


In  this  section,  the  probability  of  ambiguity,  p0,  the  array  spacings  corresponding  to 
the  optimum  realizable  arrays  for  4 < 2 < 42,  and  selected  lengths  from  42  to  100  half- 
wuvelongths  will  be  given  for  the  cascaded  end-phase  configuration.  It  was  shown  in 


t IT-  — 1. 
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Sec.  4.5  that  the  performance  of  the  cascaded  end-phase  configuration  exceeds  that  of  all 
other  four-element  array  configurations. 

It  will  be  convenient  to  use  as  the  point  of  departure  in  obtaining  these  optimum 
spaclngs  the  fact  that  if  the  overall  array  length  .<!  is  a perfect  square,  then  the  individual 
subarray  numerators  are  >/)T (Sec.  5.1).  An  example  is  £ = 16;  hero  m t = m2  = 4.  (It 
will  be  seen  that  expression  of  m2  ns  1!  without  removing  a common  factor  a , will  in 
general  lead  to  a consistent  procedure  for  classifying  arrays  in  which  £ is  not  a perfect 
square.) 

After  considering  the  so-called  square  arrays  as  those  in  which  £ is  equal  to  i2, 
i = 2,  3,  4,  . . . , it  is  natural  to  examine  arrays  in  which  £ is  equal  to  the  product  i(i  + 1). 
An  example  is  £ = 20.  Here  i ~ 4,  t + 1 = 5.  It  may  also  be  noted  that  20  is  the  geo- 
metric mean  of  16  and  26.  That  is,  “mid-square”  arrays  are  those  whose  overall  array 
lengths  are  “midway,”  in  a geometric  sense,  between  adjacent  “square”  arrays. 

A final  part  of  this  array  classification  methodology  will  group  into  one  remaining 
category  those  arrays  whose  overall  lengths  do  not  fall  into  the  other  two  categories. 

This  array  classification  may  bo  summarized  as  follows: 

Class  I:  £ - (2,/  = 2,  3,  4,  ... 

Class  II:  £ = i(i  + 1),  i - 2,  3,  4,  ... 

Class  III:  £ --  any  length  not  falling  into  Classes  I or  II. 

The  synthesis  procedure  is  based  on  the  development  in  Sec.  5.3,  and  is  mechunized 

as  computer  program  CLSIII  (see  Appendix  B for  collected  programs),  In  the  remainder 
of  this  section,  un  oxumplo  of  the  use  of  CLSIII  for  ouch  of  the  array  classes  will  be 
given.  Then  pu  for  vurious  array  lengths  vs  tf(’n.pR  will  be  tabulated. 

Class  f 

The  first  examplo  synthesizes  the  optimum  realizable  array  for  £ = 25.  The  proce- 
dure la  best  described  by  tin  examination  of  a sample  printout  from  CLSIII.  Table  6-2a 
shows  the  CLSIII  printout  for  £ = 25,  when  un  mt  value  of  6 (known  by  inspection,  be- 
cuuse  26  is  a perfect  square)  and  a trial  value  of  nx  = 1 is  used.  Candidate  values  of  n1 
can  be  obtained  from  Tuble  3-4,  which  provides  the  allowable  p\q  values  for  use  in  syn- 
thesizing threc'-olement  urrays.  It  will  be  remembered  from  Sec,  4.0  that  a cascaded  end- 
phase  four-element  array  can  be  thought  of  as  two  three-element  arrays. 

CLSIII  yields  for  each  trial  value  of  n2  " Jmi  (see  Eq.  6-19)  the  pan. meters 

• Ratio  m2:n2 

• Spacings:  p,  <7,  and  r In  half-wavelengths 

• Subarruy  t to  suburray  2 correlation  coefficient 
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Tabic  5-2a— Program  CLSII1  Printout  for  C = 26,  >nl  - 6;  Trial  n,  = 1 
cl'; 1 1 1 


S I'lMA-CMAMMV.L-PAIR  » 20.78A6  ELECTRICAL  DEIKEKO 
<J'IEDALL  ARRAY  LEM3TM  - 25  MALE -'ft*' IF.  LEMA  T Hp, 

Ml  ■ 5 

the  tp i al  value  gf  hi  ir.  71 


Ml  1 

M 1 **  + M2 I M2 

ALnHA 

,»LEN1TH  p**0**+r' 

C0.C0-FF. 

H-AMD 

"EDOM? 

(POT ) 

5 t 

1 **■< 

25  t 5 

5 * 

25 

l a 20 

0.261905 

1 1 .A  106 

99  999 

5 l 

1 *** 

25  1 10 

5 „ 

25 

2 ?.  15 

0. 175219 

10.2966 

99999 

5 1 

1 *** 

25  t 15 

5 , 

25 

•J  12  10 

7 .50939S-2 

10.2966 

99999 

5 i 

1 •** 

25  l 20 

5 , 

25 

A 16  5 

-2.38095S-2 

1 1 . A 1 0 6 

99999 

• Overall  compound  probability  of  ambiguity  (assuming  subarray /ambiguities 
are  Independent)  in  percent 

• Test  of  whether  the  resultant  synthesized  array  is  redundant,  i.o.,  whether 
p,  q,  and  r have  common  factors, 

It  is  stressed  that  CLS1II  provides  a synthesis  under  the  assumption  that  subarray  ambigu- 
ities are  independent.  The  correlation  coefficient  is  displayed  as  an  aid  to  the  user  to 
show  that  umong  arrays  with  the  same  independent  compound  p0,  the  array  with  the 
largest  absolute-value  correlation  coefficient  will  huve  the  lowest  overall  pa . The  effect 
of  subarray-to-subarray  p is  secondary,  os  pointed  out  in  Sec,  6.4.  To  establish  optimality 
of  a particular  array  configuration,  an  exact  calculation  requires  consideration  of  tho  array 
correlation  coefficient.  Table  5-2a  shows,  that  of  all  arrays  of  length  26  based  on  an 
m^'.ni  « 6:1  ratio,  the  array  mt:nt  © m2 : n 2,  a;  p r-  6:1  @ 25:10,  6;  +0.176219  is  the 
best.  That  is,  although  both  the  arrays  with  p,  q,  and  r spacings  of  2,  8,  16  or  3,  12,  10 
have  the  same  Independent  p0,  consideration  of  the  correlation  coefficient  in  a exact 
analysis  will  show  die  2,  8,  16  array  to  be  better. 

Of  course,  the  choice  of  = 6:1  is  unwise,  because  Table  3-4  shows  that  the  op- 
timum for  a three-element  array  when  = 6 is  either  r*  j =2,  or  3.  Tables  5-2b  and  c 
result  when  CLSIII  is  rerun  for  n j =2  and  3 respectively.  There  are  four  candidate  arrays 
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5:2  ©25:10 
5:2©  25:15 
6:3  © 25:10 
6:3®  25:15, 


all  exhibiting  the  same  pu  of  9.17%  (due  to  och-PR  “ 20.78  electrical  degrees,  as  forced 
by  the  criterion  of  Eq.  (5.28)).  Of  these  four  arrays  however,  the  array  defined  by 
5:3  © 25:15  has  the  largest  absolute  value  of  correlation  coefficient,  p = "0.289373. 
Consequently,  this  array  has  the  lowest  overall  pa  for  given  och-PR  arrays  length  2b. 
(Exact  calculations  of  pa  will  be  deferred  until  approximate  synthesis  procedures  for  the 
other  array  classes  have  been  described.) 


It  is  easy  to  show  that  in  Class  I arrays  the  optimum  n(  must  be  equal  to  or  greater 
than  mxl2.  If  Eq.  (5-37)  is  taken  as  a starting  point,  then  expression  of  nt  asm<(l  + 6()/2, 
where  S,  is  a deviation  of  nt  from  mtf 2,  the  array-to-array  correlation  coefficient  can  be 
put  into  the  form 


^ 1 + 35j  + §2  ~ 5 | 62 

P"'"  ■ (3  + 6?)1,2<:3*sa)1'2’ 


(5-44) 


Tublo  5-2b-- Program  CLS11I  Printout  for  K - 25,  mt  = 5;  Trial  nx  = 2 

550  DATA  25,5 


CL? I I I 


5 IRMA-CHANNEL-PAIR  * 20.784$  ELECTRICAL  DECREES 
PWRiJALL  ARRAY  LENOTH  * 25  HALr-VAVEUCNUTHS 
Ml  « 5 


THE  TRIAL  VALUE  0F  III  IS  72 
Ml  |NI***MSIN2 


ALPHA 

, LSNOTH 

00 . C0EKF  » 

P-AMB 

CPCT) 

REDIIN  7 

5 l 
5 » 

2 * ** 

i 

25 

t 5 
2 3 

20 

5.0Q626E-2 

10.2966 

99999 

5 i 
5 , 

2 

25 

25 

i 10 
4 6 

15 

-2.6315  BE -2 

9. 16857 

99999 

5 I 
5 , 

2 * * * 

25 

25 

l 15 
6 9 

10 

-0. 105263 

9.16857 

99999 

5 1 
5 , 

2 *»* 

25 

25 

j 20 
8 12 

5 

-0. 175219 

10.2966 

99999 
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Table  5- 2c— Program  CLSIII  Printout  for  V = 26,  m j - 6;  Trial  n | = 3 

550  DATA  25,5 
CLS1  I 1 


SIGMA-CHANNEL-PAIR  ■ 20.7846  ELECTRICAL  DEGREES 
0WERALL  ARRAY  LENGTH  «■  25  HALF-VAVELENOTHS 

Ml  » 5 

THE  TRIAL  VALUE  0F  Nl  IS  73 


Ml  lNl**»M2lN2 

1 ALPHA 

, LENGTH 

p***U***R 

C0 .C0EFF. 

P- AMD 

REDUN? 

CPCT5 

5 

t 

3 

25  l 

5 

| 

5 

4 

25 

3 2 

20 

-0. 175219 

10.2966 

99999 

5 

I 

3 *** 

25  i 

10 

5 

$ 

as 

6 a 

15 

-0 ,236542 

9.16057 

99999 

r 

5 

1 

3 * + * 

25  i 

15 

5 

# 

25 

9 6 

10 

-0 , 289474 

9. 16057 

99999 

s 

1 

3 + ** 

25  i 

20 

! 5 

* 

25 

12  a 

5 

-0 *325407 

10.2966 

99999 

I 

I 

I 

The  6 1 do  not  now  depend  on  the  size  of  m(  and  n/,  but  only  on  their  rutio.  Hence,  oven 
though  nt2  and  n2  might  bo  expressed  in  CLSIII  as  26  and  16  respectively,  if  the  common 
factor  a « 6 is  suppressed,  m2  and  n2  ure  functionully  equivalent  to  5 und  8 respectively . 

It  is  obvious  that  subarrays  in  which  m,  - 6 must  utilize  n(  = 2 or  3 in  order  to  minimize 
their  respective  subarruy  ambiguities.  But,  the  values  of  5(  thut  maximize  the  correlation 
coefficient  of  Eq.  (6-44)  must  have  the  same  algebraic  sign.  Hence,  the  must  be  posi- 
tive, which  means  thut  the  /q  must  bo  equal,  and  greater  than  m,/2. 

As  a second  example  of  a Class  I synthesis,  consider  V - 16.  Tuble  6-3  shows  the  re- 
sult of  running  program  CLSIII  for  K = 16.  It  is  soon  thut  for  m1:n1©w2:m2=  4:3 ©16:8, 
i the  urruy  is  redundant,  us  indicated  by  the  figure  2 in  the  RKDUN?  column.  This  is  the 

fuctor  by  which  p ~ 6,  q « 2,  and  r = 8 are  redundant.  Thus,  the  optimum  array  for 
y - 16  is  4:3 ©16:12,  with  urrny-to-urmy  correlation  coefficient  - 0.423077  ~ 11/26. 

One  further  observation  in  respect  to  Class  l urruys  is  that  the  optimum  subarray 
rutios  automatically  ensure  that  the  spucings  p—q—r  will  bo  maximized.  The  urruy 
4:3  ©16:12  has  minimum  spacing  q - 3,  whereas  the  urray  4:3  ©16:4  has  minimum 
spacing  p - 1.  Thus,  for  V - 16,  the  optimum  array  can  be  operated  over  u 3:1  bandwidth 
before  the  shortest  spacing  approaches  one  hulf-wavelongth,  whereas  in  the  other  two 
Hn'uyB,  the  shortest  spacing  is  already  at  one  half-wavelength.  For  separations  less  thun 
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Table  B-3— Program  CLSII1  Printout  for  C - 16,  mt  = 4;  Trial  nj  = 3, 1 

550  DATA  16,4 
CLSI1I 


SIQMA-CHANNEL-PAIR  « 25.9808  ELECTRICAL  DEGREES 
OVERALL  ARRAY  LENGTH  - 16  HALE-WAVE LENGTHS 
Ml  * 4 

THE  TRIAL  VALUE  0F  N1  IS  73 


HI iNl ***M2lN2 

ALPHA, LENGTH  p***Q*+*R  C0.C0SFF. 


P-AMB 
( PCT  > 


REDUN7 


4 i 3 ***  16  i 4 
4 , 16  31 


4 I 3 ***  16  « 8 
4 , 16  62 


4 i 3 *+*  16  I 12 
4 , 16  93 


IS  -0 • 046 134 

8 -0.40032 

4 -0.4230T7 


10.6339  99999 

9.76767  2 

10.6339  99999 


550  DATA  16,4 

clsii i 


mOMA-CHANNEL-PAIR  • 25.9808  ELECTRICAL  DEGREES 
OVERALL  ARRAY  LEN0TH  ■ 16  HALF -WAVELENGTHS 
M 1 ■ 4 

THE  TRIAL  VALUE  0F  N1  IS  71 


Ml  l N l ***M2 i N2 
ALPHA,  LEHOTH  p.M-*Q***R 


41  1 *+*  16  ( 4 
4,16  l 3 12 


4 I 1 ***  16  l 8 

4 , 16  2 6 8 


4 t 1 **  * 16  ! 12 

4,16  394 


C0.C05FF.  P-AMD 
(PCT) 

0.192308  10.6339 
8.0064IE-2  9.76767 
-3.84615E-2  10.6339 


HEDIJN7 


99999 


99999 


5 

'! 
■ i 


. < 
-) 


i 


► 
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one  half-wavelength,  mutual  impedance  effects  between  antenna  elements  may  affect  the 
phase  tracking,  thus  increasing  ambiguities, 

Class  II 

The  example  for  a Class  II  synthesis  is  for  B = 20.  The  integer  factors  of  20  are  4 
and  6 (=4  + 1).  In  the  absence  of  array -to -ur ray  correlation  It  might  be  speculated  that 
there  are  eight  “optimum"  arrays  formed  by  combinations  of  miOii  ©m2;n2  using  nu- 
merators 4 or  5 and  denominators  3 or  1 (with  m,  = 4)  and  3 or  2 (with  n\j  ~ 6).  As  in 
the  case  of  fi  ~ 16,  there  may  be  factors  common  to  m\  and  n2  (even  after  the  common 
factor  has  been  suppressed),  in  which  case  the  program  CLSIII  will  indicate  the  redundant 
arrays,  as  before. 

Table  6-4u  through  d show  the  results  of  CLSIII  for  B » 20.  For  each  listing,  the 
arrow  indicates  a candidate  array.  Examination  of  all  four  listings  shows  that  the 
array  4:3  ©20:12  has  the  largest  urray-to-array  correlation  coefficient,  and  hence,  will 
exhibit  the  lowest  pa.  The  arrays  5:2  © 20:15  and  6:3  © 20:15  have  p,  q,  r Hpacings  of 
6,  9,  5 und  9,  6,  5,  respectively,  in  eontrust  to  the  optimum  4:3©  20:12  array  which  has 
spacings  9,  3,  8. 

Tublo  5-4a“Program  CLSIII  Printout  for  fi  =>  20,  = 4;  Trial  n t - 3 

S30  DATA  20,4 

CLSIII 


S10MA-C1IANNEL-PAIH  ■ S3. 2.17V  ELECTFUCAL  DE0REES 
0WKRALL  ARRAY  LENGTH  - 20  HALF-VAVELENQTHS 
MI  « A 

THE  TRIAL  VALUE  0F  Nl  IS  73 


Ml tNl*+*M2iN2 


ALPHA, LENGTH 

C0.C0EFF. 

P-AMD 

REDUNT 

cpct> 

4 

t 

3 *** 

20 

1 A 

A 

» 

20 

3 l 

16 

-0.332875 

1 l .9703 

99999 

4 

X 

3 « ** 

20 

i a 

A 

* 

20 

6 2 

12 

-0,301771 

10.403 

2 

A 

1 

3 *+* 

20 

l IS 

A 

> 

i>0 

9 3 

8 

-0 • A l 3505 

10,403 

99999 

A 

t 

3 *** 

20 

l 16 

A 

J 

20 

12  A 

A 

-0,423659 

l 1 .9703 

A 
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Table  5-4b— Program  CLSIII  Printout  for  i!  = 20,  mj  = 4;  Trial  nj 

550  DATA  20. A 
CLS  I 1 1 


SIGMA-CHANNEL-PAIR  ■ 23.2379  ELECTRICAL  DECREES 
OVERALL  ARRAY  LENQTH  - 20  HALF -WAVELENGTHS 
Ml  - A 

THE  TRIAL  VALUE  0F  N1  IS  71 


Ml »N1 ***M2| N 
ALPHA. LENGTH 

2 

p***Q***H 

C0.C0EFF. 

P-AMB 
<PCT ) 

RE DUN 7 

A t 
A . 

1 «** 

20 

20 

1 A 
1 3 

16 

0.21183 

1 1.9703 

99999 

A I 
A . 

1 *** 

20 

20 

I 8 
2 6 

12 

0. 127257 

10.A83 

2 

A t 
A . 

1 A* 

20 

20 

s 12 
3 9 

8 

3. 181A2E-2 

1 0 • A 8 3 

99999 

A I 
A . 

1 *** 

20 

20 

3 16 
A 12 

A 

-6.05228E-2 

1 1.9703 

A 

Table  5-4e— Program  CLSIII  Printout  for  8 - 20,  ml  = B;  Trial  nx 

550  DATA  20. 5 


CLSIII 


SIGMA-CHANNEL-PAIR  » 23.2379  ELECTRICAL  DEGREES 
0VERALL  ARRAY  LENGTH  - 20  HALF-WAVELENGTHS 
Ml  » 5 

THE  TRIAL  VALVE  OF  N1  IS  72 


Ml tNI***M2tN2 
ALPHA. LENGTH 

p*+*Q***R 

C0.C0EFF. 

P-AMB 
<PCT  > 

RED1JN7 

5 t 
5 . 

2 *** 

20 

20  i 

s 

2 3 

15 

3.18 1 A2E-2 

10.A83 

99999 

5 l 
5 , 

2 *** 

20 

20  t 

10 

A 6 

10 

-6.62266E-2 

9.89A53 

2 

5 3 
5 . 

2 *** 

20 

20  l 

15 

6 9 

5 

-0.159071 

I0.A83 

99999 
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Table  5-4d— - Program  CLSIII  Printout  for  6 = 20,  m j = 5;  Trial  n j = 3 
550  DATA  20.5 


CLSIII 


SIGMA-CHANNEL-PAIH  ■ 23.2379  ELECTRICAL  DEGREES 
OVERALL  ARRAY  LENGTH  - 20  HALF -WAVE LENGTHS 
Ml  - 5 

THE  TRIAL  VALUE  0F  N1  IS  73 


Ml  lNl***M2tN2 
ALPHA, LENGTH 

P*+*Q*#*R 

C0.C0EFF. 

P-AMB 
CPCT  ) 

REDUN7 

5 » 
5 , 

3 +** 

20 

20  1 

5 

3 2 

15 

-0.190885 

10.483 

99999 

5 i 
5 j 

3 *** 
20 

20  1 

10 

6 4 

10 

-0.264906 

9.89453 

2 

5 t 
5 , 

3 *** 

20 

20  l 

15 

9 6 

5 

-0.318142 

10.483 

99999 

Class  III 

Over  any  span  of  array  length,  there  are  many  more  overall  lengths  In  Class  III  than 
in  Class  I and  Class  II.  To  avoid  unnecessary  iteration  through  CLSIII,  it  is  desirable  to 
know  which  mi  are  good  starting  choices  for  particular  lengths  K.  Consideration  of  the 
optimum  m t for  Class  I and  the  candidate  m1  for  Class  II  arrays  should  lead  one  to  the 
belief  that 

• If  ^ lies  between  i2  and  i(i  + 1),  choose  mj  « i 

• If  fi  lies  between  i(i  - 1)  and  i2,  choose  - i. 

The  following  are  examples  of  this  procedure. 

• For  « = 17; /2  = 16,  i((  + 1)  = 20,  choose  = 4 

• For  ? = 21;  /(« - 1)  “ 20, 12  - 25,  choose  mt  *-  5. 

Figure  5-1  gives  a geometric  interpretation  to  this  procedure.  The  trial  m}  for  Class 
III  are  given  with  question  marks  following  the  trial  integer,  anticipating  a later  discussion 
of  lengths  where  this  intuitive  trial  m j fails  (these  cases  are  “pathological”  in  the  sense 
that  they  are  associated  with  mpS]  pairs  that  are  not  as  near  the  Mj  = m\t  2 criterion 
as  their  neighbors  (see  Table  3-4)).  In  any  event,  as  will  be  seen,  this  geometric-selection 
criterion  for  ml  in  Class  III  arrays  is  natural,  and  stems  directly  from  the  fact  that  the 
cascaded  end-phase  configuration  is  based  on  a common-factor  concept  that  most  simply 
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Fig.  B-l—  m\  for  Claw;  I and  11;  candidate  m i for  Class  III:  3 < K < 42 


expresses  the  subarray  parameters  in  terms  of  the  overall  length.  This  is  in  contrast  to 
array  representations  that  have  been  used  by  others. 

Tables  5-5a  and  b shows  the  result  of  the  application  of  CLSIII  to  the  example 
fi  ---  21.  The  optimum  array  is  6:3©  21:10,  6; -0.258501. 

Exact  pa.  Including  Correlation 

Two  additional  computer  programs  are  given  in  Appendix  B,  named  AMBIG1  and 
AMBIG2.  AMBIG1  is  configured  to  solve  the  relations  beginning  in  Sec.  5.4  at  Eq.  (6-39) 
for  arbitrary  mlf  m2,  n2,  <v,  and  arbitrary  channel-pair  phase  errors  2.  °i,3>  and 
a l 4.  AMBIG2  is  configured  to  provide  pu  for  arrays  under  the  assumption  that  all 
channel-pair  phase  error  variances  are  identical.  Thus,  AMBIG1  has  the  nature  of  an  ex- 
perimenter’s tool,  giving  the  deoigner,  for  example,  the  ability  to  examine  the  effect  of 
putting  higher-quality  components  in  one  channel-pair.  AMBIG2  is  more  useful  for  tab- 
ulating pa  of  various  arrays  over  a range  of  channel-pair  phase  error  distribution  one-sigma 
values. 


Table  6-6  provides  pa  vs  channel-pair  phase  error  for  several  arrays  of  if  = 16,  20, 
and  25.  AMBIG2  was  used  to  calculate  the  pa.  Figure  5-2  is  a plot  of  pa  vs  channel-pair 
phase  error  for  l!  = 16,  and  Fig.  6-3  is  the  corresponding  plot  for  2 = 25.  In  both  cases, 
because  of  the  proximity  of  pa  vs  channel-pair  phase  error  for  certain  arrays,  only  two  curves 
are  shown.  Table  5-6  makes  these  relationships  clear. 
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Comparison  of  Approximate  and  Exact  Synthesis  Procedures 

It  was  emphasized  in  the  beginning  of  this  subsection  that  an  approximate  synthesis 
procedure  was  sought,  in  part  to  eliminate  the  necessity  of  calculating  a bivariate  density 
function  for  every  array  configuration  analyzed.  An  additional  reason  was  to  keep  the 
theory  behind  the  basic  synthesis  of  an  array  from  being  obscured.  Table  5-7  shows  the 
approximate  performance  of  the  procedure. 

This  table  is  a list,  over  the  range  3 < t!  < 42,  of  the  optimum  as  selected  by 
three  different  criteria. 

• Geometric-mean  criteria— Fig.  5-1 

• CLSIII  program 

• AMBIG2  program. 

A check  appearing  in  the  “comments”  column  means  that  all  three  methods  indicated  the 
same  (and  correct)  mj.  in  the  case  of  Class  II  arrays,  the  inability  of  the  two  simpler 
criteria  to  indicate  the  proper  cannot  be  faulted,  as  they  do  indicate  the  proper  mj  — 
i or  / + 1— within  the  constraint  of  not  considering  array-to-array  correlation. 


Table  6-5u— Program  CLSIII  Printout  for  S = 21,  m\  = 5;  Trial  «j  = 3 

550  DATA  21,5 
CLSIII 


S10MA-CHANNEL-PAIR  - 22.6779  ELECTRICAL  DE0REES 
OVERALL  ARRAY  LENGTH  « 21  HALF -WAVELENGTHS 
Ml  ■ 5 

THE  TRIAL  VALUE  0F  N1  IS  73 


Ml  lNl+*+M2tN2 
ALPHA, LENGTH 

P***Q***R 

CO.C0EFF. 

P-AMB 

<PCT) 

RKDUN7 

5 i 

5 , 

3 * ** 

21 

21  l 

5 

3 2 

16 

-0.187155 

10.2619 

99999 

5 i 
i , 

3 *+* 

21 

21  1 

10 

6 A 

1 1 

-0.253501 

9.57561 

99999 

5 » 
5 , 

3 +** 

21 

21  t 

15 

9 6 

6 

-0.312255 

10.0414 

3 

5 1 
5 , 

3 *** 
21 

21  i 

20 

12  3 

1 

-0 . 3 A 1 022 

1 1 .3056 

99999 

! 


/ 


t 


fl 
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Table  5-5b— Program  CLSI1I  Printout  for  8 = 21,  m5  ~ 5;  'rrial  rij  - 2 

550  DATA  21 <5 
CLSIII 


S JGMA-CHANNEL-PA JR  - 22.6779  ELECTRICAL  DECREES 
0VERALL  ARRAY  LENGTH  - 21  HALF- WAVE LENGTHS 
Ml  « 5 


THE  TRIAL  VALUE  0F  N1 

IS 

?2 

M l *N  1 ***M2 IN2 
ALPHA/ LENGTH  P+*#Q***R 

C0.C0EFF. 

P-AMB 

(PCT) 

REDUN7 

5 t 2 *** 
5/21 

21 

t 5 
2 3 

16 

3 • 62235E-2 

10.2619 

99999 

5 t 2 *** 
5 / 21 

21 

l 10 
4 6 

1 1 

-S.67442E-2 

9.57561 

99999 

5 1 2 *** 
5 / 21 

21 

l 15 
6 9 

6 

-0.146944 

10.0414 

3 

5 i 2 *** 
5 / 21 

21 

i 20 

8 ie 

1 

-0.21803 

1 1 .0056 

99999 

The  situation  is  different  for  array  lengths  between  31  und  42.  The  reason  can  be 
deduced  from  Table  3-4.  For  - 6,  the  only  allowable  are  either  1 or  B.  Now,  the 
ideal-unrealizable  quadratic  form  for  mi  = 6 is  (6 >6-  6-3  + S'S)^2  - 5.1962.  By  con- 
trast, the  two  realizable  quadratic  forms  are  (6*6-  6*1  + 1»1)W2  = (6*6--  6*5  + 6,5)1/2  = 
5.6678.  For  mx  - 6,  the  ideal-unrealizable  quadratic  form  is  (6*6-  5*2.5  + 2,6,2.5)1/2  =■ 
4.3301.  The  realizable  quadratic  forms  form!  = Bare  (5*  5 - 5*3  + 3*  3)*/ 2 = (5*6 -6*2 
+ 2*2)1/2  - 4.3589.  For  m i - 7,  the  ideal-unrealizuble  quadratic  form  is  (7-7  - 7*3.5 
+ 3.6 • 3.6)1/15  = 6.0622.  The  realizable  quadratic  forms  are  (7*7-  7*4  + 4*4)1/2  = (7*7 
-*7‘3  + a*3)1/2  = 6.0828. 

The  calculations  above  show  that  the  subarray  1 ambiguity  variable— (see  Eq.  (5-1))— 
for  mj  = 6 (realizable)  is  proportionally  farther  from  its  ideal-unrealizable  variable  than 
are  the  realizable  variables  for  m!  = 5 or  ttij  - 7 from  1 heir  ideal-unrealizuble  counter- 
parts. The  practical  impact  of  this  deficiency  for  amr  using  =*  8 is  that  as  the  array 
length  migrates  further  in  either  direction  from  8 = 36,  eventually  an  array  utilizing 
mi  = 6 or  ffl[  = 7 will  perform  better  than  one  using  = 6,  even  though  the  "goomotric- 
mean  criteiia  of  Fig.  5-1  are  met. 

The  forogoing  discussion  illustrated  one  way  in  which  the  approximate  analysis  leads 
to  incorrect  conclusions  on  the  optimum-realizable  array  for  a given  length.  Another  way 
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Probability  of  Ambiguity  (pc)  vs  ^ch-pr  ^or  Several  Class  I and  Class  II  Four-Element,  Cascaded  End-Phase  Arrays 
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is  expressed  in  percent. 


Criterion  A:  mi  selected  by  “geometric  mean.” 

Criterion  B:  rrtj  selected  by  Program  CLSUL 

AMBIG2:  verified  by  exact  analysis  program,  at  a Och-PR  more  likely  to  be  used  in  practical  applications. 
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in  which  the  approximate  procedure  loads  to  the  wrong  conclusions  is  as  follows.  The 
approximate  procedure  is  based  on  the  attempt  to  force  the  overall  pu  toward  a value  of 
8.9%— by  forcing  the  chunnel-pair  phase  error  to  assume  a value  based  on  overall  array 
length— without  regard  to  correlation  of  ambiguities  between  the  suborrays.  As  an  ex- 
ample of  the  effect  of  correlated  ambiguities  on  array  design,  consider  i!  = 37.  On  the 
basis  of  the  ovorull  urray  length-square  root  criterion,  one  would  probably  evaluate  candi- 
date arrays  based  on  Wj  = 6 and  = 7 by  program  CLSIII.  This  approximate  proce- 
dure indicates  the  following  two  arrays, 


6:5  ©37:18,6 
7:4  ©37:21,7 


pa  - 10.4165% 
p„  = 10.3100% 


% Program  CLSIII, 
%J  °CH-PR  “ *7.0 


(assuming  independence  of  subarray  ambiguities) 

When  program  AMBIG1  is  applied  to  the  same  array  parameters,  and  to  the  same  channel- 
pair  phase  error  distributions,  the  results  are 


6:5  ©37:18,6;  0.46517  pa  = 9,79: 

7:4 ©37:21,7;  0.265752  pa  = 10.18’ 

(including  effects  of  subarray  correlation) 


,7912%'! 

.1879%] 


Program  AMBIGl, 
®CH'PR  = 17.0848° 


For  selected  values  of  chunnel-pair  phase  error  that  are  more  likely  to  be  used  in  u system 
design,  Progrum  AMB1G2  yields 


6:6 ©37:18,6,  0,46517 
7:4@37:21, 7:  0.256752 


= 0.1946%,  1.1632%  1 
>j„  » 0.3171%,  1.4693% J 


Progrum  AMBIG2, 

«OII-l*R  ~ 12°> 

respectively 


(including  effects  of  subarruy  correlation) 


The  above  tabulations  show  that  on  the  basis  of  the  approximate  procedure  employ- 
ing progrum  CLSIII  one  might  be  tempted  to  choose  the  urruy  based  on  nij  - 7 as  the 
best  performer.  However,  more  precise  evaluation  of  the  arrays  with  the  aid  of  AMBIGl 
and  AMB1G2  (once  the  simpler  procedure  exemplified  by  CLSIII  is  used  to  identify  can- 
didate urruys),  shows  the  optimum  array  to  be  based  on  - 6. 


Tabulated  p„  vn  OcH  I’li 

Table  5-8  gives  pa  vs  ('cii-pk  f°r  Ihc  optimum-realizable  four-olomont  cuscuded  end- 
phase  array  over  un  1!  range  from  4 to  42  half-wavelengths.  The  arrays  have  optimum  p, 
q,  und  r spacings  such  tb  it  p„  of  0.1  to  15%  for  the  0(jh-i>k  specified  results. 

It  can  he  seen  from  the  table  that  the  optimum-rculizuble  urray  for  array  lengths 

24,  30,  and  36  actually  have  u pa  greater  thun  tho  urruys  one  half-wavelength  longer,  or 

25,  31,  und  37  half-wavelengths,  respectively.  This  fuel  has  apparently  not  been  reported 
previously  in  the  literature  on  multielement  interferometers.  The  reason  for  this  behavior 
is  thut  lengths  24,  30,  and  36  ure  highly  composite  numbers,  i.o., 
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Table  5-8— pa  vs  Array  Length— Optimum-realizable  Four-element  Arrays— 4 < £ < 42  (pa  expressed  in  percent) 
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24  = 2*2*2*3, 
30  = 2‘3*6, 

36  = 2’2*3*3. 


There  are  fewer  degrees  of  freedom  in  synthesizing  these  arrays,  than,  for  example, 
in  C = 26  = 2-13.  In  each  of  the  lengths  24,  30,  or  36,  m2:n2  is  either  much  greater,  or 
much  less  than  2:1.  This  is  in  marked  contrast  to  the  example  6 = 26,  where  the  m2:n2 
ratio  is  26:15,  which  is  extremely  close  to  the  ideal-unrealizable  ratio  of  2:1. 

As  an  aid  to  the  system  designer  in  making  tradeoffs,  the  information  in  Table  5-8 
is  graphed  in  Figs.  5-4a  and  5-4b. 


A figure  of  merit  for  optimum-realizable  arrays  may  be  defined  as  the  ratio  of  pa 
ideal-unrealizable,  given  Och-PR  35  Per  Eq.  (5-23),  to  the  pa  optimum-realizable,  given  the 
same  OqH.PR,  evaluated  in  Program  AMBIG1.  That  is, 


Figure  of  Merit  = f 


Po-ID  I gCH-PR>  Eq.  (5-23) 
Po-RE  I °CH-PR>  Eq.  (6-23) 


0.088281 

Pa-RE  I aCH-PR*  Eq.  (5-23)* 


(5-45) 


Figure  5-5  shows  the  figure  of  merit,  P vs  array  length  for  4 < C < 42.  Several  relation- 
ships may  be  deduced  from  the  graph. 


• With  the  exception  of  the  ideal-unrealizable  array  for  6 = 4,  Class  I (fi  = i2)  arrays 
whose  m1  are  even  do  not  have  as  great  a F as  arrays  whose  are  odd 


• Class  II  (6  = i(i  + 1))  do  not  have  as  great  a P as  the  adjacent  Class  I arrays— 
exceptions  are  6 = 16,  36 


• Arrays  whose  lengths  are  prime  have  P greater  than  the  mean  P of  0.898— the 
only  exceptions  are  6 = 7, 13 


Table  5-9  gives  pa  vs  oCH.PR  for  the  optimum-realizable  four-element  cascaded  end- 
phase  arrays,  Classes  I and  II  over  an  6 range  from  4 to  100.  The  tabulated  pa  are  graphed 
in  Fig.  6-6.  Figure  5-7  shows  the  figure  of  merit  for  these  arrays.  It  should  be  noted  that 
as  the  array  length  exceeds  23  half-wavelengths,  the  pa  exceeds  0.1%,  for  0cH-PR  = 12°. 

If  OcH-PR  *8  reduced  to  10°,  arrays  up  to  length  31  are  realized  before  a pa  of  0.1%  is 
exceeded.  Arrays  longer  than  this  require  some  form  of  calibration  (see  Sec.  3.1)  If  the 
larger  pa  associated  with  uncalibrated  arrays  longer  than  this  are  unsuitable  for  a particular 
application. 


5.6  Concluding  Remarks  on  Four-element  Array  Synthesis 

The  synthesis  concepts  presented  in  this  section  for  four-element  arrays  enable  the 
ESM  system  designer  to  rapidly  synthesize  arrays  of  any  length.  These  techniques  are 
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ARRAY  LENGTH,  / 

Fig.  6-B  — Fltfurn  of  merit,  F v»  array  length  for  four-element  arrays 


based  on  the  concept  of  the  ideal-unrealizable  array.  The  two-step  procedures  given,  in 
conjunction  with  the  computer-aided  techniques  exemplified  in  the  programs  CLSIII, 
AMBIG1  and  AMBIG2,  provide  the  designer  with  (a)  an  approximate  synthesis  leading  to 
several  candidate  arrays,  and  (b)  exact  analyses  to  fix  the  parameters  of  the  optimum- 
realizable  array  once  the  approximate  synthesis  has  been  performed. 

Although  them  is  a small  amount  of  trial  and  error  in  the  procedures,  this  should 
be  viewed  in  the  context  that  a brute-force  analysis  for  arrays  of  length  2 = 16  leads  to 
well  over  75  arrays,  all  of  whose  pa  have  to  be  evaluated  before  the  optimum-realizable 
array  spacings  can  be  specified. 

By  means  of  the  concepts  presented  in  this  report,  it  has  been  possible  to  answer  in 
the  affirmative  the  speculation  of  Hanson  [21]  on  the  existence  of  optimum  four-element 
arrays. 

A final  observation  (not  stressed  in  the  development)  is  that  one  need  not  be  re- 
stricted to  implementing  the  optimum-realizable  array  for  a given  length  if  one  or  more 
of  the  spacings  are  too  small  with  respect  to  the  overall  frequency  range  of  operation 
desired.  The  computer  programs  in  Appendix  B allow  the  designer  latitude  to  choose  be- 
tween array  spacings  that  will  mininitee  the  overall  pa  vs  those  which  are  close  to  the 
optimum-realizable,  but  which  will  maximize  the  minimum  interelement  spacing  in  the 
array. 
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Table  5-9—  pa  vs  Array  Length— Optimum-realizable  Four-element  Arrays— 4 < £ < 100,  Classes  I and  II 
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ARRAY  LENGTH,/ 

Fig.  5-7— Figure  of  merit,  F v*  array  length  for  four-element  array*, 
■ Class  I and  II,  4 < 12  < 100 


6.0  SUMMARY  AND  CONCLUSIONS 

This  report  presented  a theory  of  three-  and  four-element  phase-only  interferometers 
for  application  to  high  resolution,  low  probability  of  ambiguity  direction  finding. 

Although  three-element  interferometers  have  been  widely  used,  the  theoretical  funda- 
mentals of  these  arrays  apparently  have  not  been  widely  published  in  a form  accessible  to 
the  system  designer.  With  respect  to  four-element  arrays,  it  was  possible  to  establish  with 
use  of  some  geometric  aids,  in  conjunction  with  the  concepts  of  the  ideal-unrealizable 
array  and  the  aubarray-to-subarray  correlation  coefficient,  (a)  the  cascaded  end-phase  array 
as  the  optimum  configuration  for  four-element  three-integer  set  arrays,  and  (b)  the  optimum 
(i.o.,  the  lowest  probability  of  ambiguity  subject  to  zero-mean  channel-pair  errors  with  equal 
standard  deviation  in  all  channel-pairs)  four-element  array  spacings  for  arbitrary  overall  array 
length. 

Work  is  in  progress  to  extend  the  results  reported  hore  to  arrays  of  more  than  four 
elements,  and  to  define  the  improvement  in  accuracy  of  estimated  angle  when  phase  infor- 
mation from  all  the  apertures,  rather  than  from  only  the  farthest-spaced  pair  of  apertures, 
is  used. 
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It  is  hoped  that  the  theory  and  computer-aided  design  procedures  given  here  will 
stimulate  both  further  research  into,  as  well  as  wider  usage  of,  phase-only  interferometer 
arrays  in  those  applications  requiring  good  angular  resolution  over  wide  fields  of  view. 
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Appendix  A 


EXPANSION  OF  THE  NORMAL  PROBABILITY  INTEGRAL  (BETWEEN  SYMMETRIC 
LIMITS)  AROUND  THE  NORMALIZED  ARGUMENT  2.000 


This  Appendix  provides  an  expansion  of  the  normal  probability  integral  (between 
symmetric  limits)  around  a normalized  argument  of  2.000.  The  expansion  is  the  basis  of 
the  approximate  four-element,  array  synthesis  computer-aided  procedure  (see  Appendix  B— 
Program  CLSIII)  used  in  Sec.  B.S  of  the  main  body  of  this  report. 

Three  functions  from  Ref.  A1  related  to  normal  error  functions  appropriate  to  the 
expansion  desired  are 

NBS  No.  26.2.2 


P(x)  = 


\/2n 


in- 


-t2 

exp  — dt , 


(Al) 


NBS  No.  26.2.4 


Mx)  = £ exp  ~ dt  « 2 P(x)  - 1 , 


(A2) 


NBS  No.  26.2.9 


d „ ( x - m\  1 1 

s p \o—)  - s ■ -jr,  “p 


(x  - m)2 
2o 2 


1 r,  (x  - m\ 
- t * 


where 


Z{X)  ~ °XP  2 ' 

A Taylor’s  series  expansion  for  Eq.  (A2)  in  the  vicinity  of  argument  x is 
A(x  ± A)  A(x)  ± AA'(x)  + “ A"(x)  ± ~ A'"[x) . 


Now  we  have 
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^ '■  ^ i*  • b-«v  ■ 
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*'<*>  - & f,<«>  - £ '<-*) 
= 2Z(*), 

A"(x)  = 2Z'(*)  = ~2xZ(x) , 
and 

A"'(a)  = 2Z(x)[x2  - 1]  . 

The  use  of  Eqs.  (Ad,  6,  and  7)  in  Eq.  (A4)  yields 

A(x  ± A)  * A(x)  ± 2AZ(*)  - xZ(x)  ± ~ 


^ xZ(x)  ±2$  (x2-  1 )Z(x) 


■Y  *4*  A 

The  transformation  * + A -+  — ~ — results  in 


* + A 


A\  /AY  /*\  . 1 /AY  /** 


For  .v/cr  = 2,000,  Eq,  (A9)  simplifies  to 

* (2  ♦ £)  * A<2>  + [2  (£)  - 2 (fj  + (|f]  ' (A10) 

From  Ref.  Al,  exact  values  of  A( 2)  and  Z(2)  = — — = e-2  are 

Y 2?r 

A(2)  = 0.96449  97361 
Z(2)  - 0,06399  09666. 

Over  a range  of  - 0.4  to  +0.4  relative  to  a mean  value  of  x - 2.000,  the  approximate  value 
of  A(*  + A)  compared  to  the  exact  value  from  Ref.  A2  is  given  in  the  table  below. 


(a-  + A)/o 


LE3 


A(Arg.) 
A101 


LEI 


A(Arg.) 
Ref.  A2 


‘C 


0.890674 

0.890401 

1.73E-04 

0.910929 

0.910869 

6.99E-06 

0.928162 

0.928139 

1.27E-06 

0.942668 

0.942567 

7.97E-07 

0.9645 

0.9645 

0. 

0.964272 

0.964271 

8.64E-07 

0.972209 

0.972193 

1.56E-06 

0.978684 

0.978562 

8.19E-05 

0.983871 

0.983606 

2.66E-04 
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Appendix  B 

COMPUTER  PROGRAM  LISTINGS  AND  COMMENTS  ON  THEIR  USE 


This  Appendix  provides  listings  of  three  computer  programs  for  analyzing  and  syn- 
thesizing four-element  cascaded  end-phase  arrays,  along  with  examples  of  their  use.  The 
programs  are  written  in  BASIC  language,  and  have  run  satisfactorily  on  a time-shared 
system  utilizing  a Digital  Equipment  Corporation  System  10  at  the  Naval  Research 
Laboratory. 

The  three  programs  are  called  CLSI1I,  AMR1G1,  and  AMBIG2.  Descriptions  of  the 
programs  are  given  below. 


CLSIII 

Figure  B1  is  a listing  of  program  CLSIII.  The  user  specifies  the  desired  overall  array 
length,  fi  and  a trial  subarray  ratio  1 numerator,  mi  as  data  input  via  line  660.  The  pro- 
gram automatically  selects  a channel-pair  phase  error,  c»ch-PR  bused  un  overall  array  length 
i!  (according  to  Eq.  (5-23)  of  Sec.  5.3  of  the  body  of  this  report)  that  would  result  in  u 
probability  of  umbiguity  of  approximately  8.89%  for  the  ideal-unrealizable  array  of  the 
given  length.  The  program  will  then  query  the  user  for  his  trial  subarray  ratio  1.  denom- 
inator, «]_. 

Tho  program  then  prints  out  the  parameters  of  all  possible  array  configurations  for 
the  particular  set  of  K,  mv  and  n i chosen,  listing 

• ©m2:n2  ratios 

• Factor  a,  and  length  K 

• Array  spacings  p,  q,  and  r 

• Array-to-arruy  correlation  coefficient,  pa  a 

• Approximate  pa,  in  percent,  assuming  independence  of  subarray-to-suburray 
ambiguities 

• Indication  of  redundancy  in  array  spacings:  "99999”  indicates  "no  redun- 
dancy”; a small  integer  indicates  that  the  spacings  p,  q,  and  r have  this  factor  in  common. 

An  oxumplo  of  u CLSIII  printout  is  given  in  Fig.  B2  for  K - 21,  mj  = 6,  and  = 3. 
The  printout  shows  that  the  urruy  whose  spacings  aro  p c 6,  q ■-  4,  und  r - 11  is  the  best- 
performing urray  of  length  21,  glvon  = 5 and  n1  ~ 3.  Actually,  as  pointed  in  tho  main 
body  of  the  report,  this  is  tho  optimum  urray  for  length  21.  The  user  has  tho  freedom  to 
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oust  1 1 


10  DIM  X< 10>,Y( 10) /3< 101 

20  DIM  D<  10)/E<  10>/F<  lO>iP<  1O)/0OO),ZUO) 

30  READ  L/ M 1 

40  LET  K9  ■ SORTOwL) 

50  LET  S9  « 1B0/K9 

60  PRINT"R I3MA-CHANMEL-PAIR  -"(SO)1'  ELECTRICAL  DECREES" 

70  LET  Cl  - EH°< -2)/5QRTC2*3. 1415986) 

00  °RINT"0DERALL  ARRAY  LENOTH  L"HALF-'W»ELENOTHS" 

90  RRINT"M1  »")  M I 

100  LET  J9  ■ I NT ( L/M  I ) 

110  IF  J9*M1-L«0  THEN  130 
120  DC!  TO  1 AO 
130  LET  JO  > J9-I 

1 AO  PRINT 

150  PR l NT"THE  TRIAL  UAL'IE  OF  N1  IS  "J 
160  INPUT  Ml 
170  PRINT 
180  PRINT 

190  PRINT"  Ml  INI ***M2tN2" 

200  PR  I NT" ALPHA") "i " J "LENRTl!"j " P+*#Q* **R"> "C 0 . C0ISFF . " , “ P-AMU" 
BIO  ORINT  " ","  »i»  ( PCT  ) 41 

aao  RRINT 
230  PRINT 

240  F0n  J 1 TO  J9 
250  LET  A ■ H l *M l -Ml *N 1 +M 1 *N I 
260  LET  n . |.*L-L*J*M  1 +J*t1  l * J+M  1 
270  LET  C ■ (Ml -Nl >*<L-J*MI )-Nl*L 
280  LET  P(J)  « .5*C/SQRT (A*B> 

290  LET  D1  - K9/S0RT  < A > - 2 
300  LET  D2  ■ Ml+K9/S0nT<n)  - 2 

310  LET  31  « .9544997  + CI*(2*Dl  - 2*D1+D1  + DI*C1*D1> 

320  LET  32  « .9544997  ♦ C1*(2*D8  - 2+D2*D2  + D2*D2+D2) 

330  LET  33  » 31*32 

340  LET  0 ( J ) - lOO  + d-33) 

350  LET  D(J)  - N 1 * J 

360  LET  ECJ)  ■ (Ml-Nl >+J 

370  LET  F<  J ) - L - J*M 1 

380  F0R  K ■ J9  + 3 T0  2 STEP  -1 

390  LET  T *■  0(J)/K  - INT(0(J)/K) 

400  LET  II  - E(J)/K  - INT  <E  ( J > / 1<  3 

410  LET  u - F(J)/K  - If|T(F(J>/t<} 

4 20  LET  M « T + tt  + u 

430  IE  *•'  » 0 THEN  470 
440  NEXT  1< 

450  LET  3(0)  « 99999 

460  30  T0  4BO 

470  LET  3 C J ) - K 

480  NEXT  J 

490  FOR  J « I T0  J9 

500  PRINT  Mil 

510  PR  I N'r  Mt)"j")LiD(J)JE(J>)F<J)iP(J>i0(J)i3(J) 

520  PRINT 
530  PRINT 
540  NEXT  J 
550  DATA  25/5 
360  END 


Fig,  HI  — Lilting  of  program  CLS1U 
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550  DATA  21/5 
CLSI1I 


SJGMA-CHANNEL-PAIR  ■ 22.6779  ELECTRICAL  DECREES 
0VERALL  ARRAY  LENGTH  = 21  HALF-VAVELENGTHS 
Ml  ■ 5 

THE  TRIAL  VALUE  0F  N1  IS  73 


Ml iNl***M2lN2 

ALPHA/ LENGTH  P*+*Q***R  C0.C0EFF.  P-AMB  REDUN7 

<PCT  > 


5 r 3 *+*  21  « 5 

5/21  3 2 16  -0.187155  10.2819  99999 


5 s 3 +**  21  l 10 

5/21  6 A 11  -0.258501  9.57561 


99999 


5 l 3 ***  21  l 15 

5 / 21  9 6 6 -0.312255  10.04U 


3 


5 < 3 **+  21  l 20 

5/21  12  8 1 -0.341022  11.8056  99999 

Fig.  B2— Example  of  program  CLSIII  printout 


explore  the  performance  of  other  arrays  of  length  21  by  modifying  his  rn1  and  nj  inputs, 
should  he  wish  to  search,  for  example,  for  arrays  whose  minimum  spacings  exceed  those 
obtained  for  the  optimum-realizable  arrays  (from  pa  considerations)  discussed  in  the  main 
text  of  this  report. 

If  a radically  unsuitable  trial  m3  is  chosen,  say  m3  = 2,  for  P = 21,  the  program  will 
indicate  an  overall  pa  that  is  grossly  in  error,  compared  to  the  exact-independent  pa,  be- 
cause the  range  of  validity  of  the  expansion  for  the  probability  integral  for  subarray  2 
will  have  been  exceeded.  Normally,  trial  m3  will  be  chosen  by  reference  to  Fig.  6-1  in 
Sec.  6.5  of  the  body  of  this  report.  Figure  B3  shows  the  exact  probability  of  ambiguity 
(independence  assumed)  for  a four-element  array  in  terms  of  the  normalized  arguments 
for  the  individual  subarrays,  centered  on  a normalized  argument  of  2.000  for  each  sub- 
array. That  is,  if  each  subarray  has  a normalized  probability  function  argument  of  2.000, 
then,  the  pa  for  this  array  is  8.89%.  Suppose  an  array  is  characterized  by  a normalized 
argument  of  1.9  for  subarray  1 and  2.3  for  subarray  2.  The  exact  pa  for  the  overall 
array  would  be  7.766%;  the  approximate  pa  returned  by  program  CLSIII  would  be  7.767%, 
or  0.008  percentage  points  low.  It  can  be  seen  that  the  extremes  of  error  in  program 
CLSIII  occur  when  both  arguments  are  nearly  the  same  value.  That  is,  if  both  arguments 
are  1.6  (normalized),  the  exact  pa  is  20.719%,  and  the  approximate  pa  is  0.031  percentage 
points  low.  It  can  be  appreciated,  however,  that  in  most  array  syntheses,  when  the  nor- 
malized argument  of  one  subarray  is  less  than  2.000,  the  normalized  argument  of  the 
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3.25 

/ 


2.3* 

-.010 

-.034 

/ 

4.24 

3.7G 

2.2* 

TABULAR  VALUES  IN  THIS  REGION 

-.003 

-.010 

-.027 

ARE  SYMMETRICAL  AROUND  THE  DIAGONAL 
LINE-EQUAL  ARGUMENTS  IN 
AM  Ajly) 

/ 

5.48 

4.87 

4.37 

;,0002 

-.002 

-.006 

-.020 

2.1* 

7.02 

6.25 

6.84 

5,15 

/ 

-.002 

- .006 

- .026 

2.0- 

‘ 8.89 

7.96 

7.20 

6.60 

8.11 

-.0002 

.... 

-.0002 

- .002 

..  .006 

-.025 

11.16 

/ 

10.03 

9.11 

"8.36“ 

7.76 

7.29 

-ata' 

-.001 

- .001 

_-.001 

-.003 

-.006 

-.026 

13.88" 

~ 12.52 

’ 1 1.41 

10.50 

9.77 

~9.18 

8.71 

/ 

-.011 

••  .007 

••  .006 

-.006 

- .006 

-.007 

■ .013 

-.030 

17703 

16.46 

14.14 

13.00 

12.17  ’ 

11.45 

10.87 

10.41 

/ 

• .031  - .021 

-.017 

-.016 

- ,017_ 

• .017 

j:J518 

-.024 

-.041 

20.72  "irt  90 

17^38 

16.07 

16.01 

14.14 

13.44“ 

12.87 

12.42 

1 1 

> 

. . i 

1 

1 1 

1.8  1.7 

1 

1.8 

1 

1.9 

' f 

2.0 

1 ' 

2.1 

1 

2.2 

2.3 

2.4 

ARGUMENT  OF  A|(x!  OR  Ajly) 


NOTES:  Dp.-l-AiUlAjlyl;  1.8£*S2,4,  1.6S 2.4. 

2)  - ,002/8. 36  INDICATES  EXACT  pa  8,38%,  APPROXIMATE  p, 
IS  0.002  PERCENTAGE  POINTS  LESS 


Fig.  B3 -Comparison  of  approximate  and  exact  pa  for  normalized  arguments  A[(4  A2( >') 
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other  subarray  is  generally  larger  than  2.000;  Fig.  B3  shows  that  the  errors  in  the  approxi- 
mation tend  to  cancel  when  this  condition  exists. 

In  summary,  program  CLSI1I  is  suited  to  the  purposes  for  which  it  is  intended;  to 
eliminate  the  need  to  calculate  the  exact  pa  for  each  array  candidate,  and  to  provide  the 
user  with  a readily  applied  overview  of  the  performance  of  various  arrays  on  a relative 
basis. 

AMBIG1 

Figure  B4  is  a listing  of  program  AMBIG1.  This  program  provides  an  exact  calcula- 
tion of  pa  for  a four-element  cascaded  end-phase  array  upon  the  user’s  specifying  the  fol- 
lowing data  inputs: 

• Subarray  1 — (line  1110) 

W|,«i,  and  Pi i2;i ,3* 

• Subarray  2 --(line  1120) 

m2  = n2  jot , p1(3;I(4 

• Miscellaneous  parameters— (line  1130) 

a 1 , 2 » °1,3*  °T,4*  Pl,2;l,4»  a 

An  example  of  an  AMBIG1  printout  for  - 5:3,  = 21:10,  a = 5,  ol  a - oq  3 =• 

«l,4  - 12.5°,  and  p1(2.li3  = p 3;M  = Pt,2;l,4  " +0.5000  is  given  in  Fig.  B5. 

Program  AMBIG1  can  also  be  used  to  calculate  the  pa  of  arrays  in  which  the  Oj  y, 
j ~ 2,  3,  4 ure  not  equal.  For  exumple,  suppose  that 

01,2  - 10“ , 

01.3  - 14°, 

0 1.4  ="  12u. 

Suppose  further  that  the  channel-pair  phase  error  correlation  coefficients  are  known  (by 
measurement  of  the  joint  channel-pair  phase  errors)  to  be 

Pl,2;l,3  “ 0.457143, 

01,3jl,4  = 0.380952, 

0l,3;l  ,4  ~ 0.533333. 


*Pl.*:i,3  'H  0i«’  correlation  coefficient  between  channel-pair  errors  in  the  1 , 2 channel-pair  to  the  1 , 3 channel- 
pair.  The  value  of  this  coefficient  (and  the  two  below  in  lines  1 1 20  and  1 130)  is  usually  set  at  +0.5000  for 
design  purposes. 
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10  READ  Ml ,N1*P1 

SO  READ  M2,N2,P2 

30  READ  S1,S2,S3,P3,Q 

40  DIM  H(3Ij2)iU2),!1(2)j’;<2) 

50  PRXNT"F0tJR“5LEMENT  PHASE  INTERFEROMETER t " 

60  PRINT"C0MP0UND  PROS . -0F-AMBIGUITY  F0R  ARRAY  SYNTHESIZED' 

70  P«INr"BY  CASCADING  TV 0 THREE-ELEMENT  INTERFEROMETERS 

80  PRINT 

90  PRINT  "S IRMA>CH-PR-- 1 , 2 ="<  S 1 1 "ELECTRICAL  DEGREES" 

100  PRINT  "SIGMA, CH-PR-- 1 , 3 S2; "ELECTRICAL  DEGREES" 

110  PRINT  »S IRMA, CH-PR-- 1,4  S3J"ELECTRICAL  DEGREES" 

120  PRINT 
130  PRINT 

140  PRINT"  ARRAY  N0.1  PARAMETERS!" 

150  PRINT"  'LARGE'  =*"JM1,"  'SMALL*  «"iNl 

160  PRINT 

170  PRINT"  ARRAY  N0.2  PARAMETERS  > " 

180  PRINT"  'LARGE'  ■" » M2, " 'SMALL*  -*'JN2 

190  PRINT 

200  PRINT"  ARRAY  N0. 2-ARRAY  N0.1  RESOLUTION  FACT0Rl";(3 

210  PRINT 

220  LET  A1  » S0R<<<Ml*Sl)t2)-2.*Pl#(Ml*Sl)*<Nl*S2)+CCNl*S8: 
230  LET  M *»  5QR( < (M2*$2)«2>-2. *P2*(M2*S2 >*(N2+S3)+(<N2*S3! 
240  LET  P8  * M1*M2*P1*S1*S2  + N I *NR*P2*S2*S3 
250  LET  P9  - ~M1*N2+P3+S1*S3  - N1*M2*S2*S2  + P8 
260  LET  °4  - P9/(A1*A9) 

270  LET  Cl  « 1/(2*3.1415926) 

280  LET  C2  * 1 . /SQR < 1 . -P4*P4 ) 


220  LET  A1  » SQRC  < <M1*S1  ) i 2 ) -2 . *P 1 * (M 1 *S l > * (N 1 *S2 ) + ( <NUS2)t2>  ) 
230  LET  A 9 « S0R( ( <M2*52> ♦ 2 ) -2 . *P2*(M2*S2 )* (N2+S3 )+( <N2* S3 ) t 2 ) ) 
240  LET  P8  * M 1 *M2#P 1 *S 1 *S2  + N 1 *N2*P2*S2*S3 
250  LET  P9  - ~M l *N2+P3*S 1 *S3  - N1*M2*S2*S2  + P8 
260  LET  n4  - P9 /( A1 *A9 ) 

270  LET  Cl  " 1/(2*3.1415926) 

280  LET  C2  * 1 . /SQR ( 1 . -P4*P4 ) 

290  LET  A2  « A9/Q 

300  PRINT*'SIG“Z(ARRAY  N0.1)  -'* J A1  , "DEGREE S" 

310  PRINT 

320  PR  I NT" S IG -V (ARRAY  N0.2)  »") A9, "DEGREES" 

330  PR  I NT"  S IG-V  (ARRAY  N0.2)  «"  I A2,»  "DEGREE5  (RE50LVED  BY  N0.D" 

340  PRINT" " 

350  PR  INT"E  <7.*V  ) -"J  P9,  "DEGREES" 

360  PRINT" " 

370  PRINT"ARRAY-T0-ARRAY  C0RR.  C0EFF’ T • ■»;P4 

380  PRINT" " 

390  PRINT 

400  PRINT" ’* 

410  LET  U1  ■ 180/A1 
420  LET  112  » 180/A2 


410  LET  U1  ■ 180/A1 
420  LET  112  » 180/A2 

430  PR INT"PR0BABI LITY  FUNCTION  PARAMETERS!' 
440  PRINT 

450  PR I NT "ARGUMENT  F<1>  «";Ml 
460  PR  I NT  "ARGUMENT  F(2)  w»jU2 


47  0 LET  IK  l > ■ U1 

480  LET  11(2)  - U2 

490  F0R  J - 1 T0  2 

500  IF  U(J)>5.4513  THEN  630 


510  LET 
520  LET 


LET  ?,3  * 0. 

LET  T - U<J)/(2t .5) 
LET  S - T 

LET  Y2  » UCJ)*U(J)/2 
550  LET  D • 1 
560  LET  D » D + 2 


Ki|(.  B4  — Ll(tUn((<>f  progruni  AMBIOl 

121 


ROBERT  L.  GOODWIN 


570  LET  T =>  T*(2*Y2/D) 

500  LET  S = S + T 
590  IF(T/S-1E-10)>0  THEN  560 
600  LET  7,3  ■ < 2/5Q»T ( 3 . 14  1 5926)  >*S*EXP  < -Y2  > 

610  LET  7 ( J ) « <l.-73>/2 
620  00  T0  640 
630  LET  7(J)  - 2.5E-08 
640  NEXT  J 

650  LET  Al»l-  2*Z ( 1 ) 

660  LET  B 1 * 1 - 2*7(2) 

670  F0R  J « 1 T0  2 

680  IF  II(J)  >5.4513  THEN  950 

690  LET  H(0,J)  - l 

700  LET  H ( 1 j J ) — U(J) 

7 10  LET  A *-tKJ) 

720  LET  B ■ 1 

730  F0R  I ■ 2 T0  20 

740  LET  C * A 

750  LET  A — A+U{J)  -(1-1 )*B 

760  LET  B - C 

770  LET  H(I,J)  « A 

780  NEXT  I 

790  NEXT  J 

800  F0R  J m 1 T0  2 

810  L.ET  E « 1 

820  LET  F * 0 

830  F0R  K * 0 T0  20 

840  LET  X «■  P4*(3-2*J) 

850  LET  D - H<K,1>*H(K,2>*<Xt<K+1)> 

860  LET  E ■ E*(K+l ) 

870  LET  F * F +■  D/E 
880  NEXT  K 

890  LET  M3  ■ + U(2)*U(2) 

900  LET  FI  - F*Ct*EXP(-U3/2> 

910  LET  L< J ) - 7 ( 1 ) *Z ( 2 ) *•  Fl 
920  NEXT  J 

930  PRINT"L<1)  L<  2)  «»)L(2) 

940  00  T0  960 

950  LET  L< 1 ) * L<2)  ■ 0 

960  PRINT-PROS. -C0RR< l ) «"JA1,"  PR0B.-C0RR(2 > -"JBl 

970  PRINt"PR0B. -AMB(  1 ) »>";2>*«ZC1  ),*•  PR0B.-AMB(2>  *"J2*2(2) 

980  PRINT" " 

990  PRINT 

1000  LET  El  » 2.  - (AI+B1+2.*L( 1 )+2.*L(2) ) 

1010  PRINT 

1020  PRINT" " 

1030  LET  E2  ■ 100+E1 

1040  LET  E3  ■ INT(K2*10»4  + .5>/10t4 

1050  PRINT  "PR0E1,  0F  AMBIGUITY  =>"  J E 3)  "PERCENT  ( IMCL.  CORRELATION)" 
1060  PRINT 

1070  LET  D2  * 1 00* ( 1 -A  1 *B 1 ) 

1080  LET  DO  - I NT (02* 1 0 1 4 + .5>/10t4 

1090  PRINT  "PR0B • 0F  AMBIGUITY  J D 3) "PERCENT ( ASSUM I NO  INDEPENDENCE) 

1100  PRINT" " 

1110  DATA  5,3,.  5 
1 120  DATA  21 , 10,  .5 
1130  DATA  12.5, 12.5, 12.5, .5,5 
1 140  END 


Fig.  B4~  Lilting  of  program  AMBIG1  (continued) 
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1 110  DATA  5,3,  .5 

l 120  DATA  21, 10,. 5 

1130  DATA  12.5,  12,5, 12.5,  .5,5 


FOUR-ELEMENT  PHASE  INTERFEROMETER ! 

COMPOUND  PR0B. -0F- AMBIGUITY  FOR  ARRAY  SYNTHESIZED 
BY  CASCADING  TV0  THREE-ELEMENT  INTERFEROMETERS---- 


SIGMA, CH-PR--1 ,2  - 12.5  ELECTRICAL  DEGREES 
SIGMA, CH-PR--1 ,3  - 12.5  ELECTRICAL  DEGREES 
S IGMA, CH-PR-- 1 , A - 12.5  ELECTRICAL  DEGREES 


ARRAY  NO. 1 PARAMETERS! 
'LARGE’  « 5 


• SMALL*  ■ 3 


ARRA"'  NO.  2 PARAMETERS! 
•LARGE*  - 21 


•SMALL'  ■ 10 


ARRAY  NO. 2-ARRAY  H0.1  RESOLUTION  FACTOR!  5 
S 13 -2  < ARRAY  NO . 1 ) - 54,4862  DEGREES 
S IG -V (ARRAY  NO. 2)  - 227.418  DEGREES 

S IG “U (ARRAY  NO. 2)  - 45.4835  DEGREES (RESOLVED  BY  N0.I) 


E < 2 *’>f  ) *-3203.  1 3 


DEGREES 


ARRAY-TO -ARRAY  C0RR.  COEFF'T.  *-0.258501 


PROBABILITY  FUNCTION  PARAMETERS! 


ARGUMENT  F(l)  * 3.30359 
ARGUMENT  F<2>  « 3.95748 
L(l)  - 1 .07035E- 1 0 
PROD. -C0RR< 1 > - 0.999045 
PR0B.-AMBO)  - 9.54531E-4 


L(2>  - 4.10727E-7 
»R0B.-C0RR(2)  - 0.999924 
PR0B.-AMB<2)  - 7 • 57352E-5 


PR0B • 0F  AMBIGUITY  - 0.1029  PERCENT ( INCL . CORRELATION) 
PR0B.  OF  AMBIGUITY  * 0.103  PERCENT (ASSUMING  INDEPENDENCE) 


FiK,  B5  — Example  of  program  AMBIG1  printout  — e<)uul  channel-pair  phaie  error* 
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Those  two  seta  of  values  define  zero-moan  channel  errors,  according  to  Eq.  (3-7)  of  Sec.  3.2 
of  the  main  text  (for  p j 2;t  :i,  and  with  suitable  subscript  changes— the  other  correlation 
coefficients) 


; | 

i ' : r 


I 


I 

1 


a i - 8.0000°, 

a2  - (o'f(2  -o‘f)1/2  = (100-  64)l/2  - 6.0000°, 

°3  = (a?,3  - o'f),/2  = (196--64),/2  - 11.4892°, 

a4  = (»'f,4  -o f)1/2  = (144  -64)1/2  - 8.9443°. 

Figure  B6  shows  the  AMBIG1  printout  for  the  data  inputs 

mt  = 5,  n(  - 3,  Pli2ii,3  0.457143  (line  1110) 

m2  “ 21,  «2  “ P t ,3*1,4  = 0 380952  (line  1120) 

fft|2  » 10°,  oli3  * 14°,  ah4  - 12°,  (line  1130) 
Pl,2jt,4  3 0.5333S3,  nt  - 6. 


AMBIG2 

Figure  B7  is  a listing  of  program  AMBIG2.  This  program  is  similar  to  program 
AMBIG1,  providing  exact  calculations  of  pa  for  a four-elemont  cascaded  end-phase  array. 
The  details  of  the  probability  functions  ure  omitted  in  the  printout,  and  all  channel-pair 
phuse  error  distributions  ore  presumed  equal  to  aCH.PH  in  tho  inPut.  Hence,  in  AMBIG2, 
all  the  channel-pair  phuse  error  correlation  coefficients  ure  forced  to  equal  +0.5000. 

Figure  B8  shows  a printout  for  the  sumo  example  that  was  used  in  the  printout 
given  us  Fig.  B5.  The  data  input  on  line  860  of  AMBIG2  hus  the  form: 

mi  •-  5,  n,  --  3,  m2  = 21,  n2  - 10,  a 5. 

The  exuct  pa  (including  the  effect  of  suburruy-to-suborray  correlation)  is  calculated  over 
the  range  on  Ocn-PR  ^rom  ^ 26  electrical  degrees  in  1-degree  steps  in  Ihis  example. 

The  range  and  step  size  on  cun  be  varied  reudily  by  altering  lines  290  and  300 

in  the  program  us  required. 


I 


* 
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I no  DATA  5,3,  .-tiST  I A3 

1 120  DATA  21, 10, .300952 

1130  DATA  10. ,14., 12. ,.533333, 5 


AMD  10 1 


FOUR-ELEMENT  PHASE  INTERFEROMETER  I 

COMPOUND  PROB. -OF -AMBIGUITY  FOR  ARRAY  SYNTHESIZED 

BY  CASCADING  TV0  THREE -ELEMENT  INTERFEROMETERS  — 

S IGMA, CH-PR-- 1,2  » 10  ELECTRICAL  DEGREES 
S IGMA, CH-PR-- 1 , 3 ■ 14  ELECTRICAL  DEGREES 
S IGMA,  CH-PR— 1,4  ■ 12  ELECTRICAL  DEGREES 


ARRAY  NO . 1 PARAMETERS! 

•LARGE*  ■ 5 ’SMALL'  ■ 3 

ARRAY  NO. 2 PARAMETERS’ 

' LARGS  • •«  21  'SMALL'  - 10 

ARRAY  NO. 2-ARRAY  NO. 1 RESOLUTION  FACTOR!  5 

SIG-ZCARRAY  N0.1)  - 48.4U9  DEGREES 

SIG-N (ARRAY  NO . 2 > - 271.949  DEGREES 

S IG -U (ARRAY  NO. 2)  - 54.3897  DEGREES ( RESOLVED  BY  N0.1) 
E < Z*V ) --6908.  DEGREES 

ARRAY-TO -ARRAY  CORR.  COEFF'T.  »-0. 524671 


PROBABILITY  FUNCTION  PARAMETERS! 

ARGUMENT  FC  D ■ 3.7  17  87 
ARGUMENT  FC2)  • 3.30945 

LCD  "-1.59420E-13  LC2)  ■ 7.6I799E-6 

PROB.-CORRC 1 ) ■ 0.999799  PR0B . -CORR < 2 ) ■ 0.999065 

PROB.-AMBC 1)  ■ 2 . 00920E-4  PR0B.-AMBC2)  * 9.34742E-4 


PROB.  OF  AMBIGUITY  « 0.112  PERCENT C INCL.  CORRELATION) 

PROB.  OF  AMBIGUITY  » 0.1135  PERCENT (ASSUMING  INDEPENDENCE) 


BO --Kx .unpin  of  program  AMBK1I  printout— unoqual  cl.anitol-palr  phatto  arrow 
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AMBITS 


10  READ  MUNI<M2<N2<0 
20  DIM  S <50 ) 

30  DIM  H<31*2>«L<2)#U(2),Z<2) 

40  PRINT"F011R-ELEMENT  PHASE  INTERFEROMETER t " 

50  PRINT"C0MP0UND  PR0B. -0F-AMBIGUITY  FOR  ARRAY  SYNTHESIZED" 

60  PRINT”DY  CASCADING  TV0  THREE-ELEMENT  INTERFEROMETERS " 

70  PRINT 

80  PRINT"  ARRAY  N0.1  PARAMETERS!" 

90  PRINT"  'LARGE'  •"IMI<"  'SMALL'  "''INI 

100  PRINT 

110  PRINT"  ARRAY  N0.2  PARAMETERS!" 

120  PRINT"  'LARGE*  -"JM2,"  'SMALL'  -»!N2 

130  PRINT 

140  PRINT"  ARRAY  NO. 2-ARRAY  N0.1  RESOLUTION  FACT0RI",;Q 
150  PRINT 

160  LET  A 1 ■ SORT  <M  l *M 1 -M l *N 1 +N 1 #N  l > 

170  LET  A 9 - SORT (M2*M2-M2*N2+N2+M2  > 

180  LET  P9  - <M1-M1  > + <M2-N2 > -N 1 +M2 
190  LET  P4  - .5*P9/<A1*A9> 

200  LET  Cl  ■ 1 /< 2*3 . 1415926) 

210  LET  C2  * 1 ./SGR<  1 .-P4*P4> 

220  LET  A2  - A9/G 
230  PRINT 

240  PRINT"ARRAY-T0-ARRAY  COHR.  COEFF'T.  -”JP4 
250  PRINT 

260  PRINT"CP.PH.SIG .">"#1  P-C0Rn.","#2  P-C0RR. "< "C0MP* D. AMD. " 
270  PRINT"  <DEGR)'S"  M<"  'V*  <PCT>" 

280  PRINT 

290  FOR  M - 6 T0  26 
300  LET  S <M ) « 5 +<M-1> 

310  LET  U1  - 180./<S<M)*Al  > 

320  LET  U2  - 1 80  . /< S CM )*A2 > 

330  LET  M< 1 ) « U 1 

340  LET  U < 2 ) - U2 

350  F0R  J - 1 T0  2 

360  IF  U<J)>5.4513  THEN  490 

370  LET  23  ■ 0. 

380  LET  T - tl<J>/<2t.5> 

390  LET  S ■ T 

Kl({.  B7  — Listinu  of  program  AMBIG2 
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400  LET  Y2  - U(d)*U<d>/2 
410  LET  D « 1 
420  LET  D « D + 2 
430  LET  T “ T*<2*Y2/D> 

440  LET  S - S + T 

450  IF<T/S-1E-10>>0  THEN  420 

460  LET  23  ■ <2/SQRT « 3. 1 4 1 5926 > ) *S*EXP ( -Y2 ) 

470  LET  2<J)  - <l.-Z3>/2 

480  30  T0  500 

490  LET  2<J>  « 2.5E-08 

500  NEXT  d 

510  LET  A7  - 1 - 2*Z ( 1 ) 

520  LET  B7  ■ 1 - 2*2 <2 ) 

530  F0R  d - 1 T0  2 

540  IK  U(d>»5.45l3  THEM  800 

550  LET  H<0, J ) ■ 1 

560  LET  H( 1 , J ) — U<d) 

570  LET  A 
580  LET  B * 1 
590  F0R  1 » 2 T0  20 
600  LET  C « A 

610  LET  A -<1-1>*B 

620  LET  b « 0 

630  LET  H(L.i)  ■ A 

640  NEXT  I 

650  NEXT  J 

660  F0R  d ■ 1 T0  2 

670  LET  E ■ 1 . 

680  LET  F - 0. 

690  F0R  K « 0 T0  80 

7 00  LET  X « P4* < 3-2* J ) 

710  LET  D >«  HCK,l>*H(K,2>*(Xt<K+l>> 

720  LET  E « E*<K+ 1 ) 

730  LET  F « F + D/E 

740  NEXT  K 

750  LET  U3  ■ U(1)*’J<1>  + UC2>*UC2) 

7 60  LET  Fl  •«  F*Cl*EXP<-U3/2> 

770  LET  LCJ)  « 2<I)*?,C2)  + Ft 

780  NEXT  d 

790  00  T0  810 

800  LET  LC l ) ■ L<  2 ) - 0 

810  LET  El  - 2.  - CA7+B7+2.*L< 1 >+2.*L<2) ) 
820  LET  E2  • 100*E1 

830  LET  E3  - INT<E2*iO*4  + »5>/l0t4 
840  PRINT  S<M5,A7,S7,E3 
850  NEXT  M 

860  DATA  5,3,21,10*5 
870  END 

Fig.  B7— Lilting  of  program  AMBIG2  (contlnund) 
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290  F0R  M - 6 T0  21 
860  DATA  5,3,21,10,5 


AMBIG2 


F0UR-ELEMENT  PHASE  INTERFER0METER I 

C0MP0UND  PR0B.-0F -AMBIGUITY  F0R  ARRAY  SYNTHESIZED 

BY  CASCADING  TV0  THREE-ELEMENT  INTERFSR0METERS--- 

ARRAY  N0.1  PARAMETERS! 

•LARGE*  ■ 5 •SMALL'  - 3 

ARRAY  N0.2  PARAMETERS! 

•LARGE*  - 2 1 'SMALL*  - 10 

ARRAY  N0. 2-ARRAY  N0 . 1 RESOLUTION  FACTOR!  5 


ARRAY-T0- ARRAY  CORR.  COEFF'T.  «-Q. 258502 


CP. PR*  SIG  • 
<DEGR) 

#1  P-C0RR. 

#2  P-C0RR. 

C0MP*  D.AMB 
( PCT  ) 

10 

0.999964 

0.999999 

0.0037 

1 1 

0.999826 

0.999993 

0.0181 

12 

0.999421 

0.999963 

0.0616 

13 

0.99851 

0.999358 

0.163 

14 

0.996818 

0.99959 

0.3583 

15 

0.994095 

0,999026 

0.6851 

16 

0.990146 

0.99801 

1 . 1766 

17 

0.984864 

0.996385 

1 .8573 

18 

0.978219 

0.994009 

2.7409 

19 

0.97025 1 

0.990775 

3.8303 

20 

0.961053 

0.936617 

5.119 

21 

0.95075 

0.981509 

6.5929 

22 

0.93948T 

0.97546 

8.2329 

23 

0.927414 

0.968508 

10.0165 

24 

0.91468 

0.960715 

11.9201 

25 

0.901423  0.952155 

Fig.  B8  — Exumplu  of  program  AMBIG2  printout 
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